ΑΝ 
ANALYTICAL 
CALCULUS 


BY 
ΕΑ. MAXWELL 


VOLUME IIl 


ΠῚ TOA 5Π|πο180 yeonsjeuy Uy : 7 THMX VW 


CAMBRIDGE UNIVERSITY PRESS 


SIF 


$5.50 
in U.S.A. 


AN ANALYTICAL CALCULUS 
By E. A. MAXWELL 
Fellow of Queens’ College, Cambridge 


This is the third of a series of four volumes 
covering all stages of development of the 
Calculus, from the last years at school to 
degree standard. The books are written for 
students of science and engineering as well 
as for specialist mathematicians, and are de- 
signed to bridge the gap between the works 
used in schools and more advanced studies, 
with their emphasis on rigour. 

The volume treats of the theory of func- 
tions of more than one variable. Partial 
differentiation and multiple integrals are ex- 
plained in detail. Each section has examples; 
at the end of each chapter there are problems 
from school-leaving and open scholarship 
examinations. A final chapter deals with the 
tracing of curves from their cartesian and polar 
equations. This volume is most suitable for 
the first or second year at British universities. 

The fourth volume deals with differential 
equations and functions (such as Spherical 
Harmonics) which arise from them. 


Please see back flap for some press opinions of 
the series. 


£1.60 net 
in U.K. 


THE LIBRARY 


THE HARRIS COLLEGE 


CORPORATION STREET, PRESTON. 


All Books must be Returned to the College Library or 
Renewed not later than the last date shown below. 


Ju. * 


515 MAX 


A/C 033477 
UL 
30107 000 5 471 


ΑΝ 
ANALYTICAL CALCULUS 


FOR SCHOOL AND UNIVERSITY 


by 


E. A. MAXWELL 
Fellow of Queens’ College, Cambridge 


VOLUME Il 


CAMBRIDGE 
AT THE UNIVERSITY PRESS 
1966 


PUBLISHED BY 
THE SYNDICS OF THE CAMBRIDGE UNIVERSITY PRESS 


Bentley House, 200 Euston Road, London, N.W.1 
American Branch: 32 East 57th Street, New York, N.Y. 10022 


First Edition 1954 
Reprinted 1960 
1962 

1966 


= τς πὸ Ἂς; 
~~  »». 
CH Bath: 
~ " £2, ww 


PRESTON 


i 


Printed in Great Britain at the University Printing House, Cambridge 
(Brooke Crutchley, University Printer) 


PREFACE 


CONTENTS 


page vii 


IntRopDuctTOoRY Notre: COORDINATE SYSTEMS 
IN SPACE OF THREE DIMENSIONS 


PAGE 1 


CuarteR XIV: PARTIAL DIFFERENTIATION 


. The use of geometrical lan- 


guage 


2. Continuity 


3. Differentiation; introduc- 


tory example 


. Notation 
. The increment of a func- 


tion of several variables 


. Geometrical interpretation: 


Cartesian coordinates 


. Geometrical interpretation: 


polar coordinates 


. Plane Cartesian and polar 


coordinates 


- Connected ‘independent’ 


variables 


- The chain rule 
. Differentials 


PAGE 


= 
5 


8 
10 


11 


13 


16 


18 


22 
25 
28 


12. 


The differential of a func- 
tion of two dependent 
variables, in the form 


ez Cz 


. Digression on implicit func- 


tions 


. Differentials of functionally 


related variables 


. Small increments 
. Partial differential coeffi- 


cients of higher order 


. Implicit function of a single 


variable 


. The tangent to the curve 


J (2, y)=90 


. Exact differentials 
. Integrating factor 
. Taylor’s theorem for several 


variables 


Revision Exampies X, ‘ScHotarsHIP’ LEVEL 


ΟΗΑΡΤΕΞ XV: MAXIMA AND MINIMA 


1. The general conditions 
2. To distinguish between 


maxima and minima 


74 


77 


3. 


Lagrange’s method of un- 
determined multipliers 


Revistow Exampies XI, University LEVEL 


PAGE 


79 


85 


vi CONTENTS 
CHAPTER XVI: JACOBIANS 


PAGE 
1. Introductory example 90 6. Ratio of areas under trans- 
2. The Jacobian defined 95 formation (Approximate 
3. The chain rule for Jacobians 97 theory ) 
ἡ 7. Ratio of volumes under 
4, The ‘reciprocal’ theorem 98 transformation (Approxi- 
5. Dependent functions 100 mate theory) 


REVISION Examptrs XII, UNIveRsItyY LEVEL 


Cuarptrer XVII: MULTIPLE INTEGRALS 


- Double integrals 

. Notation 

. Triple integrals 

. The evaluation in practice 
of multiple integrals 

5. The double integral as a 

volume 


mw he 


109 6. ‘Elements of area.’ Double 


115 in 
115 7. ‘Elements of volume.’ 
Triple integrals 


118 8. The elements by Jacobians; 
change of variables 
120 


Revision Examptes XIII, University LEVEL 


CuapterR XVIII: THE SKETCHING OF CURVES 


1, Towards a general method 

2. The method for small values 
of z 

3. The method for large values 
of x 


4. To locate the multiple points 
of the curve f(z, y) = 0 


13 ὅ. Some rules for finding the 
asymptotes in simple 
146 τον 
6. The radius of curvature at 
a multiple point 
154 7, The sketching of curves 
from their polar equa- 
168 tions 


REVISION EXAMPLES XIV, ‘ScuoLaRsHiP’ LEVEL 


ANSWERS TO EXAMPLES 


INDEX 


PAGE 


102 


104 
106 


121 


127 


132 


136 


170 


173 


175 
181 


188 


195 


PREFACE 


This volume continues for two or more variables the work under- 
taken for one in the preceding volumes. Though it seems wise to 
relax standards somewhat, the aim is still reasonable rigour; for 


example, the formula Ϊ f(r, θ) rdrd@ for a double integral in 


polar coordinates and the formula | Ϊ Ϊ f(r, 0, d) r? sin θαγαθαῴ for 


a triple integral in spherical polar coordinates are derived by 
methods intended to provide a convincing proof without recourse 
to the usual ‘small elements’ of apparently approximately rele- 
vant values. I hope that those who pass from this book to the 
greater exactitude of analysis will be able to do so with some 
appreciation of the nature of analytical argument, 

I have allowed myself the pleasure of a chapter on the sketching 
of curves, which is rather out of fashion at the moment. There is 
no better field for developing that interplay of intuition and logic 
which is half the thrill of mathematics, and a revival, kept 
within bounds, is to be desired. As far as I know, the treatment 
(though foreshadowed by Salmon) has not been developed in 
such detail in text-books at this level. 

I must again express my thanks to Dr Edmonds and Dr Cassels 
for their valuable help, and to pupils who have assisted me in 
the checking of answers. My gratitude is also due to the staff of 
the Cambridge University Press for their care and attention. 

KE. A.M. 


QuEENS’ COLLEGE, CAMBRIDGE 
March, 1954 


Opportunity has been taken while re-printing to make a few 
changes and corrections. These are usually slight, but where use is 
being made of this and the earlier printings together, attention 
should be drawn to pages 92, 154, 176, 179. 


E. A.M. 
November, 1959 


INTRODUCTORY NOTE 


COORDINATE SYSTEMS IN SPACE OF 
THREE DIMENSIONS 


The object of this note is to explain briefly how some of the common 
coordinate systems are defined in space of three dimensions. For 
subsequent developments, reference must, of course, be made to 
a text-book of analytical solid geometry. 

Let O be a point fixed in space, and Oz a fixed line through it 
(fig. 104). For convenience of reference, we shall allude to Oz as the 


Att 


Fig. 104, 


‘upward vertical’. Through O draw the plane perpendicular to Oz; 
for reference, we call this plane ‘horizontal’. 

In the horizontal plane through O take two perpendicular lines 
Ox, Oy, forming a pair of coordinate axes there in the normal sense 
of plane Cartesian geometry. (In the diagram, the axis Ox appears 
te occupy its usual position in the plane—straight across the page.) 
The three mutually perpendicular lines Ox, Oy, Oz form a set of 
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Cartesian axes in space, and the position of a point P may be deter- 
mined by coordinates relative to them, exactly as in the familiar 
case of two dimensions: 

Through P draw lines PL, PM, PN perpendicular to the planes 
yO0z, zOx, xOy respectively. The values of PL, PM, PN, with appro- 
priate signs, fix the position of P precisely, and the values 

x=LP, y=MP, z=NP 
are called the Carrmstan coorpinates of P for that set of axes. 

It is customary to fix the positive senses along the axes by the 
‘right-handed corkscrew’ rule: a corkscrew, turning from positive 
Oy to positive Oz, is driven along positive Ox; turning from positive 
Oz to positive Ox, along positive Oy; turning from positive Ox to 
positive Oy, along positive Oz. 


Two other systems in common use may be described by com- 
pleting the ‘box’ of which the coordinate planes yOz, zOx, xOy form 
three faces while P is the vertex opposite to O. The lines PL, PM, 
PN appear as edges, and the vertices opposite L, M, N are called 
ζ΄, Μ', N’ respectively. 

Consider first CYLINDRICAL COORDINATES, in which the point N 
in the plane xOy, hitherto identified by its Cartesian coordinates 
τ, ἡ, is identified alternatively by its polar coordinates referred to 
O as pole and Oz as initial line. If we write 

p=ON, ὁ = Lx0N, z= NP, 
then p, ¢, z are called the cylindrical coordinates of P. 

The angle ¢ may also be thought of as the angle described in the 
positive sense (as determined by a right-handed corkscrew) when 
a plane rotates about Oz, from the position zOz into the position 
20}. 

It is customary to take p as positive, though there are times when 
the unrestricted sign proves more convenient. 

Consider next SPHERICAL POLAR COORDINATES, in which the 
position of P is determined as follows: 

(i) it les in a plane (namely, zOP) making an angle ¢ with a 
fixed plane (namely, zOx) after rotation about a fixed axis (namely, 
Oz); 

(ii) it lies, in that plane, on a line through the origin O making 
an angle @ with the fixed axis; 

(iii) it is at a distance r from O. 


COORDINATE SYSTEMS 3 
Thus, in the diagram (fig. 104) 
¢ = Lx0N, @ = 2:0}, r= OP. 


It is customary to take r as positive, though (as for p) the un- 
restricted sign is sometimes more convenient. 
The relations between the sets of coordinates are clear from the 


diagram. Thus zg =pcos, y =psing; 
a=rsin@dcos?, y=rsin@sing, z=rcos0; 
p=(e+y*), r= να γ" 2); 
Φ -- tan (γ|), @ = tan (p/2). 


Fig. 105. 


An alternative form of diagram (fig. 105) may help to make the 
‘spherical’ nature of the coordinates r, 0, ¢ somewhat clearer. 
A sphere, of centre O and radius 7, cuts the plane Oy in the ‘hori- 
zontal’ circle shown and cuts the plane zOP in a circle of which the 
‘upper’ semicircle is indicated. The ‘turning’ which defines the 
angles ¢, @ may be identified easily. 


CHAPTER XIV 
PARTIAL DIFFERENTIATION 


In the earlier volumes of this book, attention was directed mainly 
towards the properties of functions of a single variable; we must 
now examine functions of several. There are, of course, many of 
these properties which develop by natural extension; but the 
increased scope afforded by the additional variables introduces 
difficulties of which the functions readily avail themselves. A 
searching study, even of quite elementary features, would quickly 
lead us beyond the limits of the present work, and must be found 
elsewhere. The intention is to present a picture, accurate as far 
as it goes, of how the theory extends, so that the reader will be 
able, on the one hand, to make use of his knowledge, and, on the 
other, to see later what is involved when he comes to a fuller 
examination of the details. 

We should, however, state at once that it seems necessary in one 
or two places to relax the standard of rigour at which we aimed for 
a single independent variable. The language of ‘small elements’, 
‘negligible quantities’ and so on will be adopted more freely when 
clarity in the picture might otherwise be sacrificed. 

One other point deserves mention. It is now reasonable to assume 
a greater maturity of outlook and experience than could be expected 
for the first volumes; in particular, the reader should have some 
familiarity with more advanced topics such as determinants and 
coordinate solid geometry. The subject-matter of other branches of 
mathematics normally studied about the level implied by this 
volume will usually be incorporated without further reference. 


1. The use of geometrical language. Suppose that 
f(x,y, 2, ...) 
is a function of the independent variables 
A Ἐν «εν 


We shall often picture the numbers 2, Y,2Z,... as the coordinates of 
a point referred to ordinary rectangular Cartesian axes. (Since this 
will seldom be done except for two variables x,y or three variables 
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x,¥y,%, ideas involving space of more than three dimensions will 
arise only in their most elementary form.) Thus our analytical ideas 
may appear clothed in geometrical language. 

For example, the function 

(1/xy), 
is defined for all non-zero values of x,y, and we may speak of it as 
‘defined for all points of the (x,y)-plane except on the axes of 
coordinates’; also, the function 
1 

is defined for all values of x,y,z not simultaneously zero, and we 
may speak of it as ‘defined for all points of space except the 
origin’, fd 

Again, we shall meet functions defined only for limited ranges of 
values of the variables. Thus, in a theory of real variables, the 


function Jil—ai—yt—a%) 


is ‘defined for all points in the region enclosed by the sphere 
x*-+-y*+ 2% = 1’, and the function 


v(x? + y? — 42) 


is ‘defined for all points outside the circle x?+y? = 4x’. We shall 
use such language without further explanation when it seems to 
arise naturally. 


2. Continuity. A function f(x) of a single independent variable 
« is continuous at a point x = aif the difference between f(z), f(a) 
remains small for all values of x sufficiently near to «; and this 
informal statement of the principle of continuity suggests at once 
the basis for its extension. For the sake of exposition, consider the 
case when the number of independent variables is three. 


aes f(x,y; 2) 


be a given function of the three independent variables z, y,z, 
defined for all points in a region of space which includes the point 
P(x, 8, y). The function f(z, y, z) is said to be continuous at the point 
(x, β,γ) if, given any positive number e, we can find a positive 
number ἡ (depending on ¢) such that the difference 


| f(x,y, z) —f(«, 2,7) | 
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remains less than ¢ for all points (x, y, z) whose distance from (a, β, y) 
is less than 9; that is, if 


| f(x,y, 2) f(a, B, y) | «ε 
νί( -- αὐ + (y—B)? + (2-y)"} «η. 

A slight modification is sometimes useful. We have interpreted 
“near («, β, γ)᾽ as ‘lying within a sphere of radius η΄; but the sphere 
may, if we wish, be replaced by a cubical box of edges 27, giving the 
alternative inequalities 

[στα] «η, ἰν-- β] «ἡ, |z-y|<7. 
As may be expected, the sum, difference and product of two con- 


tinuous functions can be shown to be continuous; so also is their 
quotient at any point where the denominator does not vanish. 


whenever 


To illustrate the implications of this concept of continuity, let 
us consider a function of two variables x, y. The expression 


xy? 
a? + y? 
is defined for all values of x, y not both zero, and we take the func- 


tion f(x, y) based on it, but defined also to have the value zero when 
z=y=0. Thus 


oi both 
F(%,y) “Bay (x,y not both zero), 


= 0 (x = y = 0). 
The function is obviously continuous away from the origin, and 
we now discuss the continuity for x = y = 0. 
Take a given positive number e, and draw the circle of centre 


the origin and radius ἡ. Let Q(x,y) be any point, other than the 
origin, inside the circle. In terms of polar coordinates, we may write 
x=rcos0, y=rsin8g, 
where O<r<y. 
tnd 2 
Then f(e,y) = TEE Coos’ 
= r* sin? @ cos? 0, 

the cancellation being justified since r>0. But 


r<%y, sin?@ 0520 «1, 
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and so, for all points (x, y) within the circle, 


f(x,y) <7. 
In particular, if we take 7 = Je, then 


| (wy) |<e. 
Hence the function is continuous at the origin, being numerically 
less than ¢ at all points within a distance γε of it. 


As a second example, consider the function defined by the 
relations 


fle,y) =” 


Pty (x, y not both zero) 


=0 (x=y=0). 


To investigate continuity at the origin, take a given positive 
number ε, and draw a circle of centre the origin and radius k. Let 
Q(x,y) be any point, other than the origin, inside the circle. In 
terms of polar coordinates we have, as before, 

r? sin 0 cos@ 


F(z,y) a 2 


r 
= sin@ cos @. 


But it is Not now possible to ensure that this value shall be less 
than ¢ for all points within the circle; for the function sin Θ᾽ cos@, 
or ᾧ βίη 20, takes all values between — 4 and 3 as 0 varies, so that, 
for given small ε, there will certainly be points inside any circle 
(however small its radius) for which f(x,y) exceeds ε. For example, 
f(x,y) = 4 at all points except the origin on the line y = x, however 
close to the origin these points may be. In other words, the number 
ῃ of the definition cannot be obtained, and the function is discon- 
tinuous. 
It may be remarked that the formula 


f(x,y) = sin @ cos θ 


shows that the function is, so to speak, ‘continuous for approach to 
the origin along either of the axes’, its value for points of either 
axis being zero; but for approach along any other line there is a 
discontinuity, namely, a break from the value sin @ cos 0, at points 
other than the origin, to zero at the origin. Thus it is continuous as 
& function of « when y is fixed at zero, and vice versa; but not 
otherwise. 
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EXAMPLES I 


Investigate the continuity at the origin of the following functions, 
given that, in each case, the value of the function at the origin 
itself is zero. 


χϑῃϑ Ι 
l. a + ἡ 2. a? + yy?” 
3 fae ‘ee 
. t+ yt” . xt +-y*° 


3. Differentiation; introductory example. We have already 
studied in considerable detail the variation of functions of a single 
variable. For more than one variable the picture is naturally more 
complicated, for the values of the function are then affected by 
changes in each of the variables, and the whole variation arises as 
a combination of several effects, each making its own contribution 
to the total. 

Some idea of what is involved may be gathered by considering 
how the volume of a given quantity of gas depends on its tem- 
perature and pressure jointly. It is well known that (under suitable 
conditions) the volume v, the temperature ¢, and the pressure p are 
connected by the formula 

pu = Rt, 


where # is a numerical constant. Thus 


If t, p undergo small changes to values ¢+ δέ, p+ dp, the volume 
varies to a value v + dv, where 
R(t + dt) 
dv = --- - -, 
υτ Ῥ τ ὃρ 


and so, by direct subtraction, 
sya R(pdt—tdp) 
p(p + dp) 


Our concern is, of course, with small variations, and we assume 
that dt, dp are so small that their squares and their product are 
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negligible in comparison with p. By the binomial theorem we have 
the approximate relation 


prop p Pp 
“3-9. 
p p 
R op 
that δυ = — (» δὶ -- ἐδ (1-2). 
so tha pl? Pp) p 


Neglecting squares and products of small quantities we obtain the 
formula of variation 


R 
dv sone γον 


approximately. 

The significance of this formula is that it gives an expression for 
dv which is linear in the increments δέ, dp of the independent 
variables. 

There are two special cases which deserve attention: 

(i) Variation at constant pressure. In an experiment under 
conditions of constant pressure, the value of the increment dp is 
actually zero, and so the formula reduces to the relation 


dv = — dt, 
; | dv R 
leading to the relation ἐν 


This is, of course, the differential coefficient, with respect to t, of the 
function v= Ri/p, calculated on the assumption that p is constant. 

(ii) Variation at constant temperature. For constant temperature, 
the value of δὲ is zero, so that (approximately) 


Rt 


dv = a 
or ou Εἰ 
dp 


and this is again the differential coefficient, with respect to p, of 
the function v= Ré/p, calculated this time on the assumption that 
tis constant. 

2 M III 
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4. Notation. Suppose that 


w=f(x, y,2z) 
is a function of three independent variables x, y, z. On differentiating 
w with respect to x, on the assumption that y, z are treated as constants 
during the process, we obtain a differential coefficient which it is 
customary to write in the form 
ow 
Ox * 

Ow ow 
ay’ Ὁ; denote differentiations with respect 
to y (on the assumption that z, x are treated as constants) and with 
respect to z (on the assumption that x, y are treated as constants). 

These coefficients are called the PARTIAL DIFFERENTIAL CO- 
EFFICIENTS of w with respect to x, y, z respectively. Other notations 
in common use are hee 
and 4 Wy Wy 
For special purposes, we shall occasionally use the symbols f,, f/, f. 

The partial differential coefficient of f(z,y,z) with respect to x 
is the limit, if it exists, 

flim fe thy, ἘΠ 
h->0 


Similarly the symbols — 


y, 5 remaining constant; with similar formulae for ἤν; f,. 


ILLUSTRATION 1. If w = xy%2?, 


OO nay OW ogee © otis 
then a - οἰ ν2᾽, τ Υ - ye 
InLusTRATION 2. If z= οὐ: ρηδψ, 
Oz ; 0z 
— = ax ... 2 ax 
then ay "- oem sin by, ay be cos by. 


EXAMPLES II 


Evaluate ἕω ὍΣ for each of the following functions: 


Ox’ Cy 
Lk. ate’. 2. 23 /y?. 3. xsiny. 


4. coszy. 5. e* cos 2y. 6. e* +e, 
7. 8, 8. (1+2)%e-, 9. y Ju. 


INCREMENT ll 


Evaluate — cng ἴω for each of the following functions: 


0x’ Cy’ Oz 
10. ay4/z*. 11, at+y% +27, 
12. e*siny cosz. 13. sinz siny sinz, 


14, (1+2)% εν, 15. wtan— yz, 

5. The increment of a function of several variables. We 
now apply the ideas illustrated in §3 to functions in general. For 
conciseness we restrict the statement to functions of two variables, 
but the extension may easily be completed when required. 


a z=f(x,y) 
be a given function of two independent variables x, y. Suppose 
that x, y receive small increments to the values x+ dx, y+dy and 
that z consequently assumes the value z+ dz, where 


χ δ = f(~+dx,y+ dy), 
so that dz = f(u+dx, y+ dy) —f(x, y). 
It follows, by inserting and cancelling the term f(x, y + dy), that 


bz = {f (a+ du, y + dy) — f(x,y + dy)} + (Πα, y + by) f(x, y)}- 
Now the two functions f(x + dx, y + dy), f(x, y+ dy) differ only in 
respect of x, the second variable retaining its value y+ dy in each. 
Considered together, they are effectively functions of the single 
variablez, as if the variable y were reduced temporarily to the status 
of a parameter. We therefore apply the mean-value theorem for 
a single variable (Vol. I, p. 61) and obtain the relation 


f(a+ dx, y+ dy)—f(x, y+ dy) = δι ζ(α - 0, dx, y + dy), 
where 6, lies between 0, 1, and where the suffix 2 is inserted at the 
symbol of differentiation f’ to imply that f is differentiated with 
respect to x only, the second variable meanwhile remaining con- 
stant (at value y+ dy). 
In the same way, and with similar notation, 


f(u,y + dy) -- [(,Ψ) = dy fy(x,y + θ, δῳ), 
where @, lies between 0, 1, and where the suffix y is inserted at the 
symbol of differentiation f’ to imply that f is differentiated with 
respect to y only, the first variable meanwhile remaining constant 
(at value x). Thus, in all, 
dz = Oxf,(u +0, dx, y + dy) + dyf, (x, y +O, dy). 


2-2 
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[Note. For a function w of three variables x, y, z, we should obtain 
similarly the relation 


dw = dxf) (x+0, dx, y+ dy,z+ δ) 


+ OY γα, ν + Og 0y, 2+ dz) 
+ dzf (x, y,z+0, 62), 


and so on.] 
For example, if z=/(2z,y)=2%y?, 
then file, y) = 32?y?, 


ἔα, Ψ) = 2a*y, 
so that fi(x+ 0, dx, y + dy) = 3(x+ 6, dx)? (y + dy), 
Fy(@, y+ θ, δ) = 2x*(y +, dy). 
What we have proved is that, if 
dz = (x + dx)? (y + dy)? — xy?, 
then there exist numbers 0,, ,, each between 0, 1, such that 
dz = 3(x+ 0, dx)" (y+ dy)? δα + 223(y + 0, dy) dy. 

This is as far as we can go with exactitude, but an approximate 
discussion, when dx, dy are regarded as very small, is of importance 
for two reasons: first, that it is a case which often occurs in applica- 
tions; secondly, that the approximation sets the stage for an 
extension of the idea of differentials already used (Vol. I, p. 42) for 


functions of a single variable. 
We now make the assumption that the two partial differentia 


ΡΟΝ filey), file,y) 


are continuous functions; we do not want our approximations to be 
disturbed by ‘jumps’ in these values. Under this assumption, the 
difference between the products 


Oxf (x+0,dx,y+dy), dy fy (x.y + G2 dy) 
and Oafi(x.y), Syfi(x,y) 
will be very small indeed. We write 
fi(ut+ 0, dx, y + dy) —fi (x,y) = €, 
ἤνία,ν + 92 5y) —f(@,y) = 6.» 
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where €,, 69 both tend to zero as dz, dy tend to zero simultaneously. 
be dz = {fi (x,y) +6} du + {fi (x, y) + €9} Oy. 

The APPROXIMATION then assumes the form 
δὲ = fila,y) δὲ +f i(o,y) ὃν, 
or, in more suggestive notation, 


Oz ΟΖ 

= δα + By Oy. 

An important feature of this approximation is that the partial 
Oz 


ey 


dz = 


differential coefficients a filz.y), —=f)(%,y) are evaluated at 


the point (a, y) itself. 
6. Geometrical interpretation: Cartesian coordinates. 


Let Ox, Oy be the axes for a system of rectangular Cartesian 
coordinates in a horizontal plane. This is illustrated in the diagram 


P (x,y, 2) 


Fig. 106. 


(fig. 106), where the reader may regard himself as looking ‘down’ 
upon axes drawn in the usual position. Draw a straight line Oz 
vertically upwards. 

Given a point P in space, let the vertical line through it meet the 
plane zOy in Q, and draw QR perpendicular to Ox. Denote by 
x, y, 2 the lengths OR, RQ, QP respectively; then the triplet x, y, z 
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may be used as coordinates for the point P in space, just as the pair 
x, y is used for a point in a plane. If P is the point (x,y,z), then 
Q is the point (x, y, 0) in the horizontal plane. 

The coordinate z gives the height of P referred to xOy as zero 
level. In particular, if x, y, z are connected by a relation 


a f (x,y), 
where we assume f(x,y) to be a single-valued function defined for 


each pair of values of x, y, then, as x, y (and consequently z) vary, 
the point Q@ moves about the plane zOy, while P describes the 


Fig. 107. 


surface whose height at any point is equal to the corresponding 
value of the function. We say that this surface represents the function 
f(z,y). 

For instance, it is an easy example on the theorem of Pythagoras 
to show that the function 

z= + /(1-—2?—y’) 

is represented by the hemisphere of centre O and unit radius lying 
above the plane xOy. 2 

We now assume, for convenience of language, that Ox is due 
east and Oy due north. We shall also regard the surface z = f(x, y) 
as a ‘hill’ and P as the position of a ‘climber’ moving about on it. 

Suppose that the climber is at the point P (fig. 107) defined by 
the values x, y of the easterly and northerly coordinates, and that 
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he wishes to climb to the point P’ defined by 7+ dz, y+ dy. The 
whole crux of the difference between functions of one variable and 
functions of two lies in the fact that, whereas for one variable the 
motion along the curve which represents the function is defined all 
the way, for two variables the surface may be traversed by an 
innumerable choice of paths. Moreover, each way of leaving P will 
demand a gradient all of its own. The partial differential coefficients 
correspond to selective choices from the available paths. 

One way of passing from P to P’ which can easily be described in 
mathematical terms is, first, to move the distance dx easterly to 
B, and then to move the distance dy northerly. The climber thus 
describes in succession the two arcs PB, BP’ of the diagram. 

Let us now suppose that P’ is very close to P, so that the arcs 
PB, BP’ are very small. The arc PB may be regarded as almost 
straight, and so the ‘rise’ between P and B is proportional to the 


length dx; say dz (easterly) = α δὰ. 


Similarly, BP’ is almost straight, and so the ‘rise’ between B and 
P’ is proportional to dy; say 
6z(northerly) = fdy. 
If dy is the total ‘rise’ between P, P’, then 
dz = dz (easterly) + dz (northerly) 
= αδα Ὁ βδυ. 

We add that, if the climber had gone first northerly and then 
easterly, following the course PD, DP’ in the diagram, then, for 
distances so small that the paths may be regarded as straight, 
PBP'D is (in normal cases) a parallelogram, and so 

dz = Pdy+adz 
for the same values of α and f. 

Two simple observations now complete the illustration. Geo- 
metrically, a, β' are the gradients of those curves which are the 
sections of the hill in the easterly and northerly directions respec- 
tively. Analytically, we see, by putting dy = 0, that « is the ratio 
ὃζ : dx calculated on the assumption that y is constant; thus 


a 
“δ᾽ 
and, similarly, f= S 
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Hence the partial differential coefficients τ : are identified as 
the gradients of the surface z=f(x,y) 
in the x- and y-directions respectively. 


EXAMPLES III 
1. Show that the function 
z= 1—,(2e—2?-—y?) 
is represented by an ‘inverted’ hemisphere. 

Prove also that the gradients in the 2- and y-directions at the 
point (7, y,z) are in the ratio (x—1):y, and that these gradients are 
equal only for points on a certain diametral plane. 

2. Prove that the gradient at the point (x,y,z) of the surface 

z2=1—2z?— 4y?, 
for motion in the plane y=mz, has the value 
— 2x (1+4m?)/,/(1+m2). 


7. Geometrical interpretation: polar coordinates. Re- 
ferring again to§ 6, the position of the point Q in the horizontal plane 
may be described alternatively by means of polar coordinates r, 0. 
When thisis done, the height z appears as a function of r, 0 in the form 


z=g(r,0). 

The path from P to P’ then falls naturally into the two parts PB 
‘radially’ and BP’ ‘round a circle’ (fig. 108). The gradients are 
quite different from those described in § 6. 

If Q, Q’ in the plane have coordinates r, 0 and r + dr, 0 +66, where 
dr, 60 are regarded as small, the lengths of QA, AQ’ do not differ 
greatly from dr, 700 respectively. Also PB, BP’ are almost straight 
lines. 

As before, we indicate the ‘rises’ by the notation 

dz(radial) =‘rise’ from P to B, 
ὃζ (circular) =‘rise’ from B to P’. 
Denoting by α΄, β' the gradients up PB, BP’, we have the relations 
dz(radial) =a’QA =«a'dr, 
6z (circular) = f’AQ’ = f’rd8, 
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so that the total ‘rise’ dz is given by the formula 
dz = α' dr+ fr dé. 


Since = ; = are the limiting values of ὃὲ --- δὺ (@ constant) and 
ὃζ + 00 (r constant) respectively, we have the formulae 
C2 , 
gees 
Oz ; 
ap = Pr 


Fig. 108, 


Hence the values of the gradients are 


Oz 
5;" radially, 
VS gos ; 
and 00" circularly’. 
ILLUSTRATION 3. Suppose that the height at the point (2, y) is z, 
where z = ax? + by? 
The gradient for motion in the x-direction is 
Oz 
a = 2ax ; 
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and the gradient for motion in the y-direction is 
Oz 
aye 
For ‘radial’ and ‘circular’ directions, however, we transform 
to polar coordinates, so that 
z = r*(a cos? @ + bsin* 0). 


The gradient for motion radially is 


preteen 


ai = (ax? + by’); 


and the gradient for motion in the ‘circular’ direction is 


loz =r(—2acos@ sin 8 + 2bsin @ cos 0) 


r 00 
|. Se b) xy 
εἰσ στ, 
These four gradients are usually quite distinct. 


EXAMPLES IV 


Oz 02 OZ loa 
Evaluate — ax’ Dy? Br’ 700 for the surfaces: 
1, 2 = Qzy. 2. 2= 23-7, 


8. Plane Cartesian and polar coordinates. Suppose that, in 
a given plane, P is a point which (for convenience) we take in the 
first quadrant (fig. 109). The Cartesian coordinates x, y and the 
polar coordinates r, Θ are connected by the relations 
r= γα" Ἐν), 9 = tan (y/x), 
x = rcos0, y =rsind. 
Thus r is a function of the two variables x, y, where 


r= ων"). 
Hence a πάν 1: 
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Similarly 6 = tan (y/z), 
00 -Ἅ 
so that πεν aye 
eS 
ὃν x+y 


The partial differential coefficients of x, y are simpler, namely 


and 


Fig. 109. 


[It is important to observe that relations such as of 1 = 


Ox or 
do not hold.] 

These partial differentiations may be illustrated graphically. 
Consider, for example, variation with respect to x, in which y 
remains constant (fig. 110). 

If P moves to a near point P’ in such a way that y is constant, the 
line PP’ is parallel to the z-axis. Let the circle with centre O and 
radius OP cut OP’ in U. Then, if δα, dr, 60 are the increments in 
«, Ὁ, 0, we have the relations 


éx= PP’, ér=UP’, ὃθ τοι -- ΖΡΌΡ, 


the latter being negative since @ increases in the counter-clockwise 
sense, 
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The arc PU cuts the radius OP’ at right angles, and, if P’ is very 
close to P, we may regard the are PU as almost straight. Then 
PUP’ is a right-angled triangle, in which 2UP’P = 0+ 60. Hence 

. = cos UP’P = cos (+00) 


approximately. In the limit, as dz>0, the quotient dr + δὰ tends 


to the value = and 60 tends to zero. Hence 
ἐν = cos? = = : 
Ox r 
agreeing with the calculated result. 


Fig. 110. 


Again, the length of the arc UP is γί — 60), or —r 60. Also we have 
the approximate relation 


er. 5 
pp’ = sin UP'P, 
foo OC, 
or “3 sin (+60), 
6d | ὁ 
or oa — sin (7 δύ). 
In the limit, as da—> 0, we have the relation 
00 i. 
= = --- ϑηθ 
Ox r 
mae 
= — 


again agreeing with the calculated result, 
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To illustrate partial differentiation with respect to r, let us con- 
sider the coefficient = (fig. 111). This requires a variation from P 


to . in which 6 remains constant, so that P’ lies along the radius 
OP. Itis then an easy matter in elementary trigonometry to derive, 
for increments dz, dr, the relation 


so that, as dr 0, oe = θ τῷ τ 


Fig. 111. 


Finally, consider, for partial differentiation with respect to 0, 
the coefficient si . This requires a variation from P to P’ in which 
r remains constant, so that P’ lies onthe circle of centre O and radius 
OP (fig. 112). Draw PM, Ρ' Μ' perpendicular to the axis OX, and 
PU perpendicular to P’M’. Then the increments 60, dz are given 
by the relations 60 = ZPO P’, 


—dx = ΜΗ, 


the minus sign arising since dx is negative. The length of PP’ is 
ré0, and the radius OP is at right angles to the are PP’. 

For the interpretation, we suppose that P’ is very near to P, 
80 that PP’ is almost a straight line. Then 


LUPP' = 4n-—ZUPO = jn -8. 
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UP ; 
But pp? = 598 UPL, 
δου ἡ. 
so that τὸ = Sin 0 
approximately. Hence in the limit, as 690, 
=-—rsin#, 


as calculated. 


Fig. 112. 


EXAMPLES V 
Illustrate graphically the relations 


or ‘ 0b κα 

i, sure. o By a 
cy _Yy Oy 

z eal 4, ap = 9059 


9. Connected ‘independent’ variables. It often happens that 
a function z= f(x, y) 


is given in terms of two variables x, y, but that those variables are 

in fact related, each being a function of a single variable ¢ (which 

might, in a special case, be ἃ or y). Then z is a function of t, say 
z=Fi(t). 


[For example, the point (x,y) might be restricted to lie on a 
curve, such as the parabola x = al?, y = 2at.] 
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To find an expression for the differential coefficient 
a 
Tel (ἡ), 
ς ᾿ 5 ΒΩ re de 
assuming the existence and continuity of the coefficients κ᾽ By’ and 
the existence of — Te? oe 


Suppose that ¢ receives an increment wu, as a result of which 2, ἡ 
become +h, y+k. Then 


F(t+u)—F() _ feth, ytk)—f(x,y) 
0 Φ, erage = Cee 
_ f(at+h, yt+k)—feyth}+ S(eyt+h —fle,y)} 
U 
_ Mf (e+ Oh, y +k) + λυ, y + Oak) 
U > 


as for the similar argument on p. 11. Hence 


F t+u —Fit h ’ k , 
ea me = yi et tO hiy +k) +~ fy(e.y +k). 
N ow proceed to the limit, letting uw tend to zero. Since 4 = 


38 WG owe ~ ae’ 


also h, k both tend to zero with wu. Hence, by the usual theorems on 

sums and products of limits (Vol. I, p. 25) and the continuity of 

ὃ 

a of , the relation becomes 
a oy 


F d 
= = Flute n+ 2a, 9), 


dF δ ἄς Y dy 
dt dudt ὃν dt’ 
This relation is also written in the form 


dz 0dzdx 4 0 dy 
dt dxdt + oy at? 


where it is understood that z is expressed in terms of t on the left 
and of x, y on the right. 


or, on rewriting, 
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ILLUsTRATION 4. Suppose that z is the square of the distance of 


a point P(x, y) on the ellipse 
a2? 
att p32 
from the focus (ae, 0). Then 
z= (w—ae)*+y". 


1 


The variables x, y are not independent, but can be expressed in 
terms of a parameter ὁ in the form 


x=acost, y= bsint. 


Then = = 2(~—ae) = 2a(cost—e), 

oz 

-- = 2y = 2bsint, 

dy νυ 

dx dy 

== —asint, = bcost, 
Hence, by the formula, 

= — 2a7(cos t—e) sint+ 2b? sint cost 


2a7e sin t — 2(a? — b?) sin t cost 
= 2a®esint—2a*e*sint cost. 
For example, the value of z is a maximum or minimum when 


dz , 
= 0; that is, when Gat an G, 
or cost = l/e. 


The latter solution is impossible (since e <1), and so the only turning 
values occur when sin? = 0; that is, when P is at an end of the major 
axis. 

ILLUSTRATION 5. Suppose that 

z=f (x,y), 

where x, y are linear functions of the variable ¢; say, 
zt=a+pt, y=b+dt. 

eee 

dt ox dt cy dt 

Gz. δὲ 


Then 
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EXAMPLES VI 


Evaluate F(t) by the rule given in this section, and verify your 
answer by calculating /’(¢) explicitly and then differentiating it, for 
each of the following examples: 

1. f (zy) = aAy*; ἃ) τα ἐδ, y = ἐξ, 

2. f(x,y) = rev; x=cost, y=sint. 

3. f(x,y) = σα Ἐν; x=cost, y=sint. 

4. f(z,y)=2e%; z=Ff, y = 2t. 


10. The chain rule. The work of the preceding paragraph may 
be generalized. Suppose that the two* variables x, y are each 
functions of the three* variables ξ, 7, €; say χ(ξ, 7, €), y(&, 7, €). Then 
a given function w(x, y) of the variables x, y can be expressed in 
terms of £, y, € as a function ρ(ξ, 7, ζ). We obtain an expression for 
the partial differential coefficient 

op 
0g” 
in terms of the partial differential coefficients of w, x, y. 

When ἔξ, 7, € receive increments δέ, dy, δζ, the variables x, y 

become 2+ dx, y+ dy respectively, and w(x, y) becomes 
w(x + dx, y + dy). 
The increment in w (or p) is (compare p. 11) 


w(a + dx, y + dy) — w(x, y) 
= δχιυ (+, dx, y + dy) + dyw)(x,y+O,dy). (1) 


Moreover, 2, y are themselves functions of £, ἡ, ζ, so that the 
corresponding increments are 


da = a(§ + 6, 9 +09, €+ δζ) —2(E, 7, Ὁ 
= δξας(ξ- 68,9 +69, 6+ δ) 
+ dy x,(£,9 + $209, 6+ δῷ) 
+ δζας(ξ,η, ζ-Ὲ $396), (2) 


with a similar expression in y. 
* The numbers, two, three, are only chosen for illustration and may be any 
Positive integers. 
3 M III 
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On substituting from equation (2), and the corresponding equa- 
tion with y, into equation (1), we obtain as the coefficient of δξ in 
the increment of w the expression 

w(x +0, dx, y + dy) x(E +4, 68, ἡ +09, 6+) 

+ w(x, ¥ 1 0, dy) με(ξ + Vy ὃξ, ῆ + δη, c+ δζ), 


where the number y, replaces the ¢, of equation (2). 
Now restrict our attention to those cases in which the partial 
differential coefficients are continuous. Then (p. 13) 


0 
εὐς(α Ἐθι dx, y+ dy) — w,(2, Y)+ey ae i 


Ox 
w) (x,y +O.dy) = wi(x,y)+e, = oe ΠΝ 
wi + φιδξ,η Ὁ 90,6480) = ENO +a = ὅξελι, 
E+ Yr BE, 0+ δὴ, C+ δ) = ψέξ, η, Ὁ) Ἐλς = Et Ay 


where A,, A, tend to zero as df, dy, δ tend to zero; so that dz, dy also 
tend to zero, and hence ¢,, 65 also. Thus the increment in w (or, to 
be more exact, in p, which is the same as τὸ but with £, 7, ζ as the 
independent variables) is do, where 


Ow Ca Ow oy 
to = (sea tay ae) 
Ow On Ow 2) 
+ δὴ 
(τ δῃ ὃν Oy 
(wet Ox Ow ὅν at) δέ 
Ox att oy of) > 
plus terms which vanish with δ, dy, df. 
Hence we have the relations, forming the chain rule, 
op _ ow Ou Ow oy 


o& ou a Oy Oe? 
ap _ ewe , dwey 
Oy Ox δὴ Gy On’ 
Op _ Ow ox 4 Ow oy 
ee ἂς OC" dy BC" 
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More generally: 

If wis a function of the m variables x1, Xp, ...,C,_ which are themselves 
functions of the n variables £,, &o, «+s &ns eA. eg w when expressed in 
terms of £1, b, .--, 18 @ function p, then 

ὃ. @ [ep ox; |. 
x An. $= 1,2, ..., ἢ). 
δὲ, ,Ξιϑα, 06, ( ) 

ΤΙ ΒΊΒΑΤΙΟΝ 6. Let w=zxy?2’, 


where *=cos@singd, y=sin@sing, z=cos¢. 
Then 
Ow dwdu Cwoy  dwoz 


δ ~ dx 00” by 00" bz 00 
= —y25 sin sin ὁ + 2xyz* cos 0 sin ᾧ + 3xy*2*.0 
= —sin?@ sin? ¢ cos? ὁ + 2 cos? @ sin θ᾽ sin® ¢ cos* d 
= (2 -- 8 51η3 0) sin @ sin® ὁ cos* d; 
Ow dwox prey 5 Ow 02 
δῴ dx od ἘΝ 36° Oz δῴ 
= y*z3 cos 0 cos ᾧ + 2xyz* sin 0 cos ᾧ — 3xy"z* sin ᾧ 
= sin?@ cos @ sin? ¢ cos'¢d + 2sin®@ cos@ sin? ὁ cos* ὁ 
—3sin?@ cos 0 sin‘ ¢ cos? ᾧ 
= 3sin?0 cos@ sin? ¢ cos® ¢(cos* ὁ —sin* ¢). 


EXAMPLES VII 


1. Given that w= αύξ, 
where a=cos@singd, y=sin@sing, z= cos¢, 
Cw ow 
evaluate θ᾽ δ᾽ 
2. Given that w= u?—v*, 
where U=e+yY+2Z, v= τῇ, 
Ow Ow dw 
evaluate ὄχ’ By? ae" 
3. Given that w = e“ cos», 
where U=2L+Y, VS=U-Y, 
Cw ow 
evaluate 55 ay 


3-2 
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4, Given that w = 2 log yz, 
where r=u?+v%, y=ut+v, z=Uu-d, 
Cw Cw 
evaluate oy? Oo’ 


11. Differentials. We now transfer our attention from the 
approximate theory of small increments to a more exact theory 
based on the concept of differentials. 

Suppose that z=f(x,y) 


is a function of the two independent variables x, y, and that these 
variables receive increments dx, dy not necessarily small. Since we 


know what is meant by the partial differential coefficients if = 
(see ὃ 4), we may form the function 
Oz ΟΖ 


which, for given increments, has a definite value. It is called the 
DIFFERENTIAL of z for the increments dx, dy, and is denoted by the 
symbol dz, so that, as a matter of definition, 


0z dz 
dz = Oe δὰ + dy oy. 
In particular, we require meanings for the differentials of 2, y 
themselves. If we regard x as a function of the two independent 


variables x, y, then One On 
se 
so that dx = 1.dx+0.dy 
= δ. 
Similarly dy = dy. 


Thus the differentials of the two functions 2, y for the increments 
dx, dy are just dz, dy respectively. Hence for given increments 
dx, dy the differentials dx, dy, dz are connected by the (exact) 
relation 

Oz dz 
dz = 55 dx + δὴ 

Returning for a moment to the geometrical illustration of §6 
(pp. 13-16), we can give an intuitive interpretation of the differ- 
entials, though a precise statement would be hard at present. The 


—_— 


dy. 
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surface (fig. 107) z = f(x, y) meets the vertical lines BA, P’Q’ at the 
points B, P’; suppose that the ‘tangent plane’ at P meets them in 
B, P’. The sections of the tangent plane in the easterly and northerly 
directions really are straight lines, and their gradients are exactly 
dee 
On’ Oy 
calculated at P. Hence the rise from P to P’ is precisely 


ΟΖ ΟΖ 

or, by definition, the differential dz. Thus the differential dz is the 
increment in z when the surface z = f(x,y) is replaced by the tangent 
plane to it at P. For small increments, the two values are approxi- 
mately equal. 


12. The differential of a function of two DEPENDENT 


variables, in the form dz = = dx+ ὅν ἀν. 
With the notation of ὃ 9 (p. 22), multiply the relation proved 


there, namely dz dzdu Odzdy 


di " ὃς dt" dy dt’ 
by the differential dt. Then 


dz Οχ ἄχ, dz dy 
au = mata 

Now z (on the left), x, y are functions of the single variable ¢, 

and we know (Vol. I, p. 43) that their differentials dz, dx, dy are 


defined* by the relations 


dim ie ae Se, dy =~! 


dt dt 
re = ; τ [ = are differential coefficients with respect to t. 


Hence the relation assumes the simple form 
oz Oz 
dz = 555 + dy dy, 


giving the connexion between the differentials dx, dy, dz. 


dt, 


whe 


_* It should be clearly understood that we are appealing directly to the defini- 
tions. We are not morely cancelling dt, though the notation dz/dt for a differential 
coefficient may create an optical illusion to that effect. 
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The significance of this result will be appreciated on referring to 
811 (p. 28), where precisely this formula appears under the hypo- 
thesis that the variables x, y are independent. The work of the two 
paragraphs thus unites to obtain the important theorem that, ἐ zis 
a function of the two variables x, y, then the relation 

dz = = dz + = dy 
is true both when the variables x, y on the right are independent, and 
when they are dependent. 


13. Digression onimplicit functions. A detailed study of the 
properties of implicit functions is difficult, and cannot be attempted 
here. The present paragraph seeks to illustrate, mainly by example, 
the ideas which the reader should have at the back of his mind in 
the work which follows. 

Consider the relation zy*z'-—1=0. 

It can be used in three ways: | 
(i) to express x as a function of y, z in the form 
= ye; 
(ii) to express y as a function of z, x in the form 
y = 2-tat; 

(iii) to express z as a function of x, y in the form 

z= aby, 

These three functional relations are implicit in the given equation, 
and we have, in fact, been able to make them explicit by direct 
solution. 

There are, of course, many implicit relations for which direct 
solution is (in practice, if not in theory) quite impossible; for 
example e*# —sin (4 Ἐν -Ὁ 2) = 0. 

In any event, there are often advantages in keeping the equation 
in its implicit form, without direct solution, and so the corre- 
sponding rules for partial differentiation must be formulated, 

Returning to the example zy*z5 = 1, we obtain the differential 
relations corresponding to the solutions (i), (ii), (iii) in the forms 

dx = — 3y~4z- dy — ὅν z* dz, 

dy = —$2-$4-4dz—12-32-tdz, 

dz = —}ta-tytda—fa-ty8dy, 
respectively. 
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Multiply the first equation on the left by τὶ and on the right by 

its equivalent y*z5; then 
xtdxe = —3y dy —5z1dz. 

Multiply the second equation on the left by 3y-! and on the right 

by its equivalent 3ziz?; then 
3y—tdy = —5z1dz—a" dz. 

Multiply the third equation on the left by 5z~1 and on the right by 

its equivalent aty*; then 
5zdz = —a1dx—3y—dy. 

It is immediately obvious that the three differential relations 
just obtained are identical. Indeed, they could have been obtained 
at once by writing down the differential du of the function 

u = ry®z—1 
intheform du =y'zidx+ 3xy2dy + d5xy%ztdz, 
and noting that du = 0 since τὸ has the constant value zero. 

This example illustrates a technique which we shall adopt 
without further comment. Suppose that, say, three variables 
ΖΦ, y, z are connected implicitly by a functional relation 

μία, y,z) = 0. 
Then we obtain the same relation between the differentials dz, dy, dz 
whether we find it directly by equating dw to zero or by first 
expressing ἃ as a function of y, z (or y of z, x; or z of 5, y). 


14, Differentials of functionally related variables. Suppose 
that three variables x, y, z are connected by some relation, so that 
either may be expressed as a function of the two others. Regarding 
«as a function of y, z, we have 


Ox Ox 
regarding y as a function of z, 2, we have 
oy, Oy. 
dy = δ, dz + δ dx; 
regarding z as a function of x, y, we have 
Oz Oz 
dz = ap ta, 


_ (As we have just explained, these three forms of the relation are 
ἴηι fact identical, the ratios dw: dy: dz being the same in each.) 
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An alternative notation for the partial differential coefficients 
is sometimes found useful as giving greater precision. We write, 
for example, a 
a) 
OY} 
to denote the result of differentiating x partially with respect to y 
when the other independent variable is z. Ἤτοι such notation, the 


three relations given above are 
dx = 5) ἀνε 5) dz, 
ramen st) de 5) ae, 
ἄς = me) det 5) ἂν 


[Thus if the implicit relation is 
at+y2+z22 = 1, 


02: y oO Ζ 
ia | holes Melee 
oy x 


The three ways of looking at the functional relation, and con- 


sequent three ways of looking at the differential relation, all result 


in a linear equation among the differentials dx, dy, dz; and it is often 
convenient to express this linear relation with symmetric notation 
udaz+vdy+wdz = 0, 
(where the symbols wu, v, w denote certain functions of the variables 
x, y, 5) leaving unspecified at this stage which of the variables are 
to be regarded as independent and which as dependent. 
With this notation, the partial differential coefficients are 


Vv 
=) oe 5) - ω 
85), τ Ὧν ἀν 
=) = os Ὥ Age 
δὰ 5 δ᾽ ν᾽ 
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This leads to another observation. Consider, for example, the 

coefficient Oe 

i 

It is calculated by differentiation on the assumption that z (the 

other independent variable) is constant, and this assumption may 

be stated in the equivalent form that dz is zero. But then the 
differential relation udz+vdy +wdz = 0 becomes simply 


udx+vdy = 0, 
and it follows, since 

xl ow 

= 
that 5) Ξε ἡ -- αν. 

ὃν)ς 


In other words, the partial differential coefficient ΕΞ =} , where zis kept 


constant, is equal to the ratio dx = dy of the ΟΝ dx, dy when 
dz is equated to zero. 
These ideas may be extended to any number of variables, and 
the resulting theorems are important in practice: 
(i) If zis a function of the variables 2, 5» ...,7,, then 
_ O% Oz Oz 


whether the variables 2,, 4» ...,%,, are independent or not. 

(ii) If there is a functional relation between n variables 
Ys Yos +++) Ym, then their differentials are correspondingly connected 
by a linear relation of the form 


U, dy, + Ugdygt ...+Undyn = 9, 
or Su,dy, = 0, 
1 


where w,, Uo, ...9 U, are functions of the ἡ) variables. 
Further, if we regard, say, y, as a function of the remaining 
variables ¥5, Ys, ...» 3... then the partial differential coefficient 


Cys 
CYe 
(calculated with 3, Y4; ---»Y, all constant) is equal to the ratio 


dy, --- ἀν, 
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calculated with dy, = dy, = ... = dy, = 0; so that 


ILLUSTRATION 7. Three variables x, y, z are such that each may be 
regarded as a function of the other two. To prove that 


(2) -» 
mor), @)-— 


where, for example, ( means that xis expressed as a function of y, z, 
Ε] 


oy 
and then z kept constant in the differentiation. 

A relation between 2, y, z implies a linear relation between their 
differentials. Thus functions a, ὃ, c of the variables exist such that 


adx+bdy+cdz = 0. 


Hence (=) - =dx+dy when dz= 
eo b/a, 
and (2) =dy+dz when dz=0 
= —a/b, 
can 6)8}-: 
In the same way, (32). = —c/b, 
(Ξ) --ale 
(,--1 
won Θ)6)8)-- 


ILLUSTRATION 8. Four variables u,t,p,v have the property that each 
one of them can be expressed as a function of any two of the others. 


To prove that Ou ιν δὲ 
(si), - (ai), * Ge), (5), 
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where, for example, (=) means that u is expressed as a function 


of t, p, and then p kept conned in the differentiation. 

A relation between three variables implies a linear relation 
between their differentials. Since u, ἐ, p are related, functions 
a, b, c of the variables exist such that 


adu+bdt+cdp = 0. 


Thus (=) =du+dt when dp=0 
re δα. 

There is also a relation between w, ἐ, v, and functions A, B, C of the 

variables therefore exist such that 


Adu+Bdi+Cdv = 0. 


Thus (=) =du+dt when dv=0 
* = —BIA, 
Ou 
and (2) =du+dv when di=0 
cv), 


=—C/A. 
This accounts for the first three partial differential coefficients in 


the relation which we have to prove. There remains (=) , for 
p 
which we need the relation between dv, dt, dp. This is found by 


eliminating du between the two differential relations, giving 


bdt+cdp _ oo ἘΦ σαν. 
α Α 
or aC dv+(aB—bA)dt—cAdp = 0. 
Hence (=), =dv+dt when dp=0 
_ bA- aB 0A B 
αὖ ad Ὁ 
(i), _ (=i), 
ot 


ὦ), ὦ), 
mito), ~ GB), +05), 03), 
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EXAMPLES VIII 


1. Prove that, if wu, y are each functions of the two variables 
x, t, and if wu is expressed (by elimination of t) as a function of 2, y, 


a Nase 


Illustrate this result by considering the cases 
(i) w=2+f, y=at; 
(ii) w=e*cost, y= e*sint. 


2. Each of the variables w, v, τὸ is a function of the two variables 
a, y. By elimination of zx, y, the variable τὸ is expressed as a function 
of v, w. Prove that 


Cudv ducv 
Cu\ Ca Oy oy Ox 
5a). Owdv Cwdv 
Ox Ou Ou o£ 


Illustrate this result by considering the cases 
(i) w=22+y*, v=2*-y*, w= 2zy. 
Gi) w= ety’, σεν, w= —ayty. 
15. Small increments. The example which follows is 
typical of many: 
ILLUSTRATION 9. The sides ὃ, ὁ and the angle A of a triangle ABC 
PIT Latin, oa Gia, Aw oi. 
A number of children (not knowing these figures) make measurements, 
with a maximum error of 2% in length and 3% in angle. Τὸ find the 


greatest consequent error in their estimates of the area of the triangle. 
The area z is given by the formula 


z= fbesin A. 


When dealing with products, it is often convenient to begin by 
taking logarithms, so we write 


u=logz = log$+logb+loge+logsin A. 
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Suppose that the errors are db, dc, dA, leading to an error dz in z or 
du in logz. Then we have the approximate formulae 


du 
ou = τῷ oe 
Ou Ou Ou 
δι = ah b+ = δον τ OA, 
so that ἌΝ ἜΗΝ 
Zz b Cc sin A 


Now the coefficients 1/5, 1/c, cos A/sin A are all positive (an im- 
portant point) and so the greatest error arises when db, dc, dA are 
all positive and at their greatest values. But we then have 
db = 26/100, de = 2c/100, 6A = 34,100, so that 


1). oi (2) 1 (2%) 9054 (34 
τ -$\i00) *c\100) ὁ sind τοῦ 
— £4 3A cot 
~ 100 100 ᾿ 


Now the angle A is measured (as in all work for calculus) in radians, 
so that A = ἐπ; also cot A = 3. Hence 


1, 4 πι8 _ 8+7,3 
2 δ = 300+ 200 = 200 ’ 
_ (8+7 73)z 
so that ieee” “ἀκ ; 


The greatest percentage error in the estimate of z is thus 
ἐ(8 Ἐπ 3) % = 63 % 
approximately. The corresponding error in actual value is 
sane, δὲ 
4$x4x5x sin 30 XT90? 
or 0-34 sq.in. 


EXAMPLES IX 
Repeat the work of this illustration (with appropriate modifica- 
tions) for the following sets of measurements: 
1, O=22, c=3, A= 48°. 
2 6=4, c=5, A= 150°. 
3%. δι. 3, c=3, A = 138°. 
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16. Partial differential coefficients of higher order. If 


2 Oe z=f (a ? y), 
then ae . ay are also functions of the two variables x, y, and possess 
their own differential coefficients. We use notation such as 


ye 
Ca \du} 0a?’ 
Ὁ (d%z\ oz 
ra (at) Sa 
ὃ (=) = 032 0 (τς ὃ 
du \dy) ϑαδν᾽ Oy \dx) 

ὃ (dz\ _ 0% 

ay (og) "δ, 


and so on. 
For example, if z = (ax+by), 
Cz a . 
then a? 3a(ax + by)*, 
oz 
a. = 36 ax+b ᾿ 
oy ( y) 
C2 2 
On = 6a oo Ss: 
28. 
ων" Ὁ = ογδε; 
ε ΞΞ 6b7(ax + by). 


Before going further, we give a proof (under the simplest con- 
ditions only) that the order in which partial differentiations are 
performed is irrelevant. In other words, we prove that, in ‘ normal’ 
cases, 


Oz 0% 
Oudy δια" 
[The proof which follows may be omitted at a first reading.] 
Let z=f(x, y) 


be a given function of the two variables x, y, and consider the 
expression 


E=f(c+h,yt+h)—f(e+h,y)—f(a,y+k)+f(x,y). 
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First we express Z in the form 
H=(f(et+h,yt+k) —f(a+h, ψ)}--ἰ ιν +k) —f(%,y)}; 
and then write u(x) =f(x,y+k)—f(x,y), 
where y is regarded as constant. Thus 
E=u(x+h)—u(z), 
and we may apply the mean-value theorem (Vol. I, p. 61) for the 
function u(x) of the single variable x, Hence 
B=hu'(%+6,h), 
where @, lies between 0, 1, the differentiation of wu being performed 
on the assumption that y is constant. Thus, by definition of wu, 
H=h{f,(x+O,h,y+k)—fi(c+Oh,y)}. 
υ(ψ) =f,(e+ Oh, y), 
where «+ 06,h is now regarded as constant. Then 
H=h{v(y +k)—vy)}, 
and we may apply the mean-value theorem for the function v(y) of 
the single variable y. Hence 


E=hkv'(y + 6,k), 
where @, lies between 0, 1, the differentiation of v being performed 
on the assumption that x is constant. Thus, by definition of », 
EB=hkf},(2+0,h,y+0,k), 
the differentiations of f(x,y) being performed first with respect to 
« and then with respect to y. 
In the same way, we could have written Z in the form 
H={f(ceth,y+k)—f(u,y+k)}—{fet+h,y) —f(x,y)} 
and so obtained the formula 
H=khfy,(u+Osh,y+0,k), 
where @,, 0, lie between 0, 1 and the differentiations of f(x,y) are 
performed first with respect to y and then with respect to x, 

Hence fyalt-+ Osh, y +O.h) = foy(v+O,h,y + 04h). 

Now let h,&-0 independently. On the assumption that the 
two second-order partial differential coefficients are both con- 
tinuous, we reach, as the limiting form of the equation, the formula 

Oi. CF 
Cyou χουν" 


Now write 


ILLUSTRATION: CHANGE OF VARIABLE [5 
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Finally, we remark that, once this result is established, we can 

easily prove more general results such as 
082 08 082 


Ge. Gt 
Hence εἰ μα 


Differentiate with respect to x, y respectively. Then 


διδὸν Oxdydx Oy dx?” Oz 02  O% ad 
and so on, ἐπ On Yanoy ~ P+ur 
EXAMPLES X = d—¢'y/z; 
a δ» Oe Fe ὋΝ Oe . ΘΑ͂» ὃν φ 
Evaluate Su? Dndy? dy®? datoy’ dwoyt for each of the following πὰ o 
functions: = ¢’ 
1, ary. 2. οὖιβ, 3. e*siny. δ). 
ἃ. wlogy. δι᾿ πον. αἰ ἀκ Multiply these equations by x, y and add. Then 
7. a/y. 8. ytan a. 9. z*y*e. ery Say oe Ἐν sat tee ty = = xg, 
InLusTRATION 10. Τὸ prove that, if west aie ἂν ὶ 
or, using the previous result, 
z = wply/z) + Wly/2), f Gia ae 
yee Oe) Soe 2° — + 2ey~—— + y*—, = 0. 
then x? 3 + 2ay ay a+ ay = = 0, dat “ραν 7 oy 
where d(y/x), γγ(ν 4) are functions of the variable y/x. ILLUSTRATION 11. (Change of variable.) Το prove that, if 
Let ¢’, y’ be written to denote ¢’(y/x), w'(y/x), the differential 85. δ; 
coefficients of ῴ, y with respect to the variable y/x. (In other words, δι" Bys? 


. . : ὃ ἀφ ; i 
if ὁ is written for y/x, then ¢’ is -τς with y/x written for ¢ after then the function z must be of the form 


differentiation.) z2=f(z@+y)+g(x—-y), 
Then 55 2 tblylz y=? se wl x) where f, g are arbitrary functions of their arguments. 
= (Wyla)} = ae = 2—y. 
Oz 0208 Oz δὴ 
= ¢'/x, seine Gx OE dx" Oy Oe 
with similar results for y. Thus, by the definition of z, na 
a ( ὧν vy ~ BE" Oy” ‘ 
ie Pe ἢ @ hit tt 
oz 2 ag’ ψ' and Ξ 
ao ar ae oy a a ὃη oy 
z , Ξ 
ΞΞ g’ ὩΣ. 7 og on” 
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It follows that the ‘operation’ 
0 


Ox 
on a function of x, y is equivalent to the ‘operation’ 
δ ὁ 
aE" Gy 
on the same function expressed in terms of £, 7; and that 
0 
oy 
᾿ ᾿ ϑ. ὃ 
is equivalent to a ay 
2 
Thus ite = τς ) (by definition) 
= Ὁ (5; +] (replacing =) 
ox \c& cy Oat 
δ. O\ (ez. & --- ὦ 
(eta) (era) (ρθε 
ὃ [dz Oz\ ΟἿΟΣ cz 
= 25 seq) + (see) 
he oz 42.0% 022 Pi 
OE “CEOn Oy 


and, by similar argument, 
ite δα , Be δε 
ὃυ" δὲ) ck 


Hence the given relation 
i ΕΒ ς 
ox® ον 
Oz 
is equivalent to ara = 0. 
is equiv aE on 
ὃ (oz 
Th (= )- 0, 
ig δὃξ δὴ 
so that a is independent of £ and therefore a function of ἢ only; say 
Oz 
- = U . 
Sa et 
Integrating, = | μ(η) ἄη- Α, 
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where A is constant with respect to ἡ, but may involve &. 
Writing [ὦ dy = (7) 
and A = f(§), 
we have the relation z=f(f)+9(%), 
or, intermsofz,y, z=f(x+y)+g(x—y). 


ILLUSTRATION 12. (Change of variable by the use of operators.) 


Given that κανῷ or ae ὅν 
1 (2 δῴ 
to express +y (as +55) 


with u, v as independent variables. 

This example is typical of many, and is given in detail. (The 
beginner may also like to solve it by straightforward partial 
differentiations.) We assume throughout that the context makes 
it clear when ¢ is being regarded as a function of x, y and when it is 
being regarded as a function of ει, v. 

By the chain rule of §10 (p. 25) 


2p θυ ϑῤδν 
Ox Cudcx dvox 


= δ (20) +2 (ay); 


= ei ἢ Ἢ 
Hence vey - 2a 42) + dary 
ee ae ee 
= wot 2uck, 
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We may regard the symbol od ο 
ὃς ὃ ΠΣ 3 =) 
ἀπ FY 
Ox Cy Od 
as an OPERATOR in the sense that, when it ‘operates’ on the function Ἐν (x ays) 
@ (expressed in terms of x, y), the result is ad 
ret 4 yt +(e) 
What we have just proved is that the operators 4 (uss + =) , 
ee er ὟΣ ἐν 
δ oy? "Oa Toy’ fnall ye Od 9 7p ἔδῥῴ OF. δῴ 
acting on any function of x, y, are equivalent to the operators si % Eat “δον idl ey? en Yay" 
Qu + 208, 20° — au Hence we have the equation 
8 sony 28  αϑό, 06, a6 


el “ 

xg ὁ 2ay- ἜΣ la +35, 
«ΟΦ OO 20° ὃ Og 
=4(u "aya + Ma ay t oa to ΝΥ 3) 

Now consider similarly the relation 


(ἀρ 


operating on that function when expressed in terms of τι, Ὁ, 

[The next steps may be abbreviated with experience, but are 
given in full for clarity.] 

In particular, 


(Pat¥e, a) ret) 
= (20g +202) (ano! + 20 51, 


ou 
nee the functions ΤΗΝ ἀκ 2 P19 ἐῷ d th ut να; os "2 
si toe? ae toy oot ee Operas The left-hand side is 
ὃ ὃ 0 "2 , ; οῴ 
(-ξ Ἐν =) , (205 +2v 2) are respectively equivalent. Y=, (y =) 
Consider the left-hand side. 0 ( δῴ 
[This step should be noted carefully, as it is often calculated Fay (y =<) 
wrongly. ] ὃ ( δῴ 
This is equal to ~ = (#56) 
(«= +z) (-) + (εξ + =| ( ἂς ὃ ( δῴ 
dx” cy) \" cx Ga Gy) \Yey)? +25 (2X), 
δαὶ dx ἡ wt? 24 
ὃ -) ed d¢ 
+Yy— — ee ὦ 
155 ("35 -+(vevan* se) 
0 Hp δῴ 
ἜΣ 19=— - — 
sa (Vey ‘ Cs sty * i) 
can 2 " δὲῴ 
ἐν Vay) Ὡ 
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Hence ib Bbde Bb 2y 
8 _ ay 18 5 hit ms a0 ~ Μ1 
ν δ." σον ον οὃχ _% δῴ 
og 2d 06 th δῴ = τ ρεῖβ 0) Ἐπὶ (γ cos 8) 
= «(005 su +u aa a? ose). __, ag δῷ 
Oc ὃν 
Adding this to the similar equation obtained above, we have 
2 2 Hence c= ye Qay fae pet aaes yee 
24 42) δ Φ δώ 24 92) (2P , OP Buty” Oy? ὃς Yay 
t+ a(S Spe) τ (ὦ εοὴ (58+ Se) αδῴ 06 δῷ ap δή 
Adding, S40 48 = (οὐ Ἐν} 4 52) 
But δ υ = (a®@—y?)? + (2ay)* = (5 Ἐν), ey 
need a+ y® (esa) (51: 54)" a" dy 


δ, ad OS 1δῴ 10% 


so that 
ILLUSTRATION 13. The expression i cate! in polar coordinates. Gy? Ort Or γὴ δῦ 
The expression : . ILLUSTRATION 14. Homogeneous functions. A function of three 
ζ + ὃ as variables f(a, y, 2) is said to be homogeneous of degree n if it possesses 
ey the property 


flux, uy, uz) - Unf (x, 5» z) 
for some number n. Typical examples are 


a? δι 22 


is of great importance in mathematical physics, and it is important 
to be able to transform it to polar coordinates in accordance with 
the usual formulae 


ᾧ τι γοοβθ, y=rsind, e-e4e (n = 1), 
We use the method of the preceding illustration, now expressing eetwie (n= 0), 
the argument with more normal brevity. l 

By the chain rule of ὃ 10 (p. 25), 4 dys 328 (n = -- 8). 
og as og ox δῴ oy Another exampleis j2+jy+/z (n= }), 
Or Ox Or ' oy or 
where it is assumed (for real functions) that x, y, z are positive, 
δι. er Pee ae and that w also is positive. 
Cx oy We prove the following theorem: 
ba. J fiv.y,2) 

so that -- = 

εν Yay" ts a function homogeneous of degree n in x, y, z, then 

og 56, | δῆ δες 

Hence γ᾿; = (r3 4)- νι a tvz) (5 δ; * Sy)" tat Ya tga 
bs ne a τος aan od a ot od Suppose that x, y, Ζ are temporarily fixed. Then 


= Bae . 
any Tg FY ΤΣ ὩΣ" f(ux, uy, uz) 
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is a function of the single variable u. Write 
μα τὸ ξ, αὐ τ ἡ, uz=C. 


Then, in accordance with the rule for differentiation (compare 
Illustration 5 on p. 24), 


d 
ron mS) 45 δ[(ξ,η, ) dy , O(E,0, ξ) ἀξ 
re: πο du* on dut og = du 
Of(E9,0) ἀξ, 78), F(E, 7, ἢ) 
fe ae eee 
But, by definition, 
f (ux, uy, uz) =u"f (x,y,z), 
so that (x, y,z still being fixed) we also have the relation 


d = n—l 
— flux, uy, wz) τειν γα, y, 2), 
Hence, equating values of 2 f (ux, wy, uz), 


a! ¢) oye, f) οὔτέ, TOS) nur, y,2). 


In particular, this result is true for w = 1; and then &, 9, € are 
respectively equal to x, y z. Hence 


δ ᾽ ᾿- 9 4.) 8 ; 
εἶτ, ϑια ὩΣ ῸΝ 5), Fe yen, @) πὰς τον <a 
or, more briefly, ἐν at Ya a ten = nf. 


This result may be extended. By differentiating with respect to 
x, y, z respectively, we obtain the relations 


Pe γε Ῥω, of 
"ant t On στον * *teon 
*f 


of Of af 
“aa dy "δ ay" “dyde 
o*f o*f of of 


τ area +4 dyde t Gat δ: = 


"ἧς 
ng. 
m. 
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Multiply by z, y, z in turn and add, first transferring the terms in 
δ F FD vo the ight hand side, Then 


du’ Oy’ Oz 
eity sates Sat eae ΠΣ + ay 58 
~ n= nf yt ied 
= n(n—1)f. 


17. Implicit functions of a single variable. A function y of 
a variable x is often defined, not directly in the form y = f(x), but 
implicitly by means of an equation 


f(x,y) = 0 


The differential coefficients — κι of y with respect to x may 


y dy 
dx’ dx* 
then be calculated as follows: , 

Since x, y are both functions οὗ z, we have, by the chain rule, 


the relation df af 2 af dy 


Hence ἊΝ ᾿Ξ - ; 


Moreover, applying the same treatment to “ : " respectively, 
4.) 2, a 
dix \ox Ox Oy daz’ 
τί =) = OF OF dy 
da λῶν ΠΝ oy? daz’ 
80 that, on differentiating the relation 
ὦ ὦ ἂν 
0 - ant oy de 
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with respect to 2, 


( 65} 2) (<5 Pcs a dy | of dy 


dx Gedy da) ' \Oxdy  dytdx) dx dydx®  -° 


053 ΟΣ" dy o*f (ἀγλ3 ef dy 
= ΩΣ Soy et alae +s oe" : 
(<) 4-2 of FOF +(2 ) ae 
a με Pe Oy) Ox®  dOxdydudy \ex) Cy? 


(a) 
oy 
ILLUSTRATION 15. 70 find expressions for dy vy when 


dz’ ἀπ 
2*+y%—1= 0. 
Differentiating, we have 


so that ἘΞ τῷ --- τας 


Differentiating again, 


2 
or po ᾿ 


Hence a Se 


EXAMPLES XI 
dy dy. 


Find expressions for qx’ dat in each of the following cases: 


1, y?—4ax+a? = 0. 2. ax*+ βυ3--Ἰ τ Ο. 
3. 2+y° = a3, 4. xy®*+y2? = 1, 


18. The tangent to the curve f(*,y)=0. Let (p,q) be 
a given point lying on the curve whose equation is 


71(α, y) ἊΝ 


TANGENT TO THE CURVE f(z, y)=0 δ] 


The gradient at the point is 2 , where 


of 
dy oa 
ἀττῷ 
oy 


(p. 49), the differential coefficients being evaluated when 2 = Ὁ 
y = q. Hence the equation of the tangent is 


= 
y-q=- Fe Pp), 
oy 
“ of a 
or [2 τ fg 9.3] ὦ 4) 


the square brackets being inserted to imply evaluation with 
t=P,y = ἡ 


ILLUSTRATION 16. 70 find the ἜΤΟΣ of the tangent to the ellipse 


at the point (p,q). 
Of 2% of 
We have ax = - > ς κα > 
Taf _ 29 
Ox (p,q -ὦ, i (p,q) “ δὲ 


Hence the equation of the tangent is 
2 
=P (e—p) +74 (y-4) = 0, 


or, since p?/a* + q*/b? = 1, rete =], 


19. Exact differentials. We are now familiar with the formula 


dz = © a + ΙΝ dy for the differential of a function z of two variables 


x, y; but the converse problem often arises and requires an answer: 
To determine whether, when two functions f(x, y), g(x,y) are given, 
there exists a function z of x, y whose differential is given by the relation 


dz = f(x,y) ἀπ -ἰ σία, ν) dy. 
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If z does exist, then, necessarily, 


= = f(y); Ξ = σία, ν), 


and these two partial differential coefficients satisfy (for normal 
conditions) the relation (p. 38) 


Peyese) 
dy \ex) — da \dy}” 
Hence @ NECESSARY condition for the existence of the function z is that 


ὃ 
5, ὕλιο,Ψ}} = ξοίσί,ν}}. 
The condition is also SUFFICIENT: 
If two functions f=f(x,y), g= σία, y) are such that 
of og 


oy Ox’ 
then there exists a function z of x, y such that 
dz = fdx+gdy. 


[This proof may be omitted at a first reading.] 
Take the function f(x, y) and integrate it on the assumption that 
y is to be kept constant; this gives a function 


F(z,y)= Ϊ » oui” y) dx. 


In the same way, integration of g(x,y) on the assumption that z is 
to be kept constant gives the function 


G(x,y) = | g(x, y) dy. 
(x const.) 

The functions /’, G satisfy the relations (Vol. I, pp. 86-7) 

ae 

ΠΝ 

ey = 9; 
so that, in virtue of the given relation = = , we have 

or ag 


Oyea —dxdy’ 
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If we assume that F is a function of ‘normal’ type for which 
or ὁ 


‘bet , then 

ay da δυῶν er ag 
δρῦν  Sxey’ 
e(F-G)_ ὦ 

or exoy =v. 


Hence (compare pp. 42-3) the difference F —G is of the form 
F-—G=u(x)—v(y), 
where u(x), v(y) are functions of x, y only. It follows that 
F(x, y) + v(y) = G(x, y) + u(x), 
and we take either of these expressions as the function z whose 
differential we require, namely, 
2=F (x,y) +v(y)=G(x,y)+u(2). 
With this choice of z, we have 


oz OF 
δ: on 
te 20 
Cy oy ©’ 
so that fdx+gdy = δ ἄν + = dy 
= dz. 


Derrmition. When the condition 


= ie, y)} = τ tote, y)} 


is satisfied, the expression 
S(x,y) da +(x, y) dy 
ts called an EXACT DIFFERENTIAL. 
ILLUSTRATION 17, The steps of the proof may be grasped more 
clearly by reference to a particular example, say 
3u*(a + y) dx + (x3 + 3y) dy. 
Integrating the first part with respect to x, on the assumption 
that y is constant, we have 
F (x,y) = jx* + xy, 
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and integrating the second part with respect to y, on the assumption 
that x is constant, we have 


G(x, y) = αν + $y’. 


What, in essence, we have to do is to add to F(x, y) a function 


of y, and to G(x, y) a function of x, so as to yield (if possible) the 
same function z. This we do by taking 


a= fat + ay + dy". 
Then dz = (3.3 + 3a°y) dx + (a + 3y) dy, 
as required. 
The point of the proof given above is that, since 
0 0 
sq (@ + Sy) = 5 (B28 + B2%y), 


functions u(x), v(y) must exist such that 
P(x, y) + u(y) = G(x, y) + u(z). 


When they have been determined, either side of this identity gives 
a function z. 


20. Integrating factor. Although there is not, in general, 
a function z whose differential dz is equal to the given expression 


fda σάν, 
it nevertheless happens in several important instances that a 
function 
“Ξε μ(α, ψ) 
exists with the property that the product 
u(fdx σάν) 
is an exact differential. The function μία, ν) is then called an 


INTEGRATING FACTOR of the given expression. 
The condition given in the preceding paragraph becomes 


By ME) = δα HO)» 
or wie = wesoce, 

om ch (δ ὃ 
" vast a5 ~ * (Gy 2a) 


INTEGRATING FACTOR 55 


This equation, which we may expect to be of help in deriving μ, 
is unfortunately somewhat awkward to handle; but in practice 
what we often want is a solution in which μὲ is a function of x only, 
in which case the working becomes more manageable: 

If u is a function of x only, then 


τ =0 
ia ὑπ is the ordinary differential coefficient oe Then 
of ὃ 
so that on = MC 
x. 2 


If it happens that the coefficient of da is a function of x only, the 
function (x) can be found; for then 


The Illustration which follows shows how an integrating factor 
may be found—though an experienced mathematician will usually 
‘spot’ the solution in such a direct example without going through 
the routine process. 


ILLUSTRATION 18. To find the current y at time t in an electric 
circuit with self-inductance L and resistance R under an electro- 
motive force E cos pt, where L, R, HL, p are constants. 

It is known that y satisfies the differential equation 

dy 
di +Ry = Ecos pt, 
or, on multiplying by the differential dt, 
Ldy + Rydt = E cos ptdt. 

We seek to find whether the left-hand side is an exact differential. 

If so, 


ὃ 0 
δὲ (1) ἊΣ ay HY) 


or 0= &, 
which is not true in general. 


56 PARTIAL DIFFERENTIATION 


We next try to find an integrating factor u, by means of which 


the equation will become 
pL dy +pRydt = μ ἢ cos ptdt. 


The right-hand side can be integrated (perhaps with difficulty, of 
course) if w is a function of ἐ only. Let us therefore try to make the 
left-hand side an exact differential on that assumption. We must 


have 2 
ay (ML) = = (UR) 


or, since μὲ is a function of ¢ only, 


d 
L+ = pR, 
du R 
or aon zu 
Hence log x = Rt/L, 
or fim οἴ, 


no constant being introduced since any one particular value of uw 
will suffice. 
The equation is therefore 


Le®l dy + Rye dt = He® cos ptdt, 
or 1, ἀ(ν 6 1 = Εὶ e®L cos ptdt, 


or yerL = (BL) | em cos pidi + C, 


where C is an arbitrary constant. Evaluating the integral by parts 


in the usual way, we easily reach the solution 


RUL EB eRtL ; 
ye Fa Lapis oom + Lp ain pt} + Ο. 


The current y at time ¢ is therefore given by the formula 
E 
y= Biz Lap Bos pt + Lpain pt} + Oe“, 


where C is a constant whose value, for a given problem, depends on 
the initial conditions. 
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21. Taylor’s theorem for several variables. In Vol. IT, p. 44, 
we discussed Taylor’s theorem for functions of a single variable. 
In enunciating the extension we have to use two points of notation 
which we shall explain almost immediately: 


TAYLOR’S THEOREM: To prove that, if f(x, y, z) is a function of the 
three variables x, y, z, then a number 0, lying between 0, 1, can be 
found such that 


f(ath, y+j,z+k) 
= f(x,y, 2) 


+ (rn +ig tks] ΑΝ 


1 δι. πον 0λ3 
ἘΣ ΣΝ 5 f(x,y, 2) 


1 948 g\*-1 
toma lhetiggt*a) S(x,y, 2) 
l ” Se o\* ; 
+i (ig tig the) fla+Oh,y +0j,2+0k). 


It is assumed that the function and all its partial differential co- 
efficients exist and are continuous for the ranges of values to be 
considered. 


[For the ‘shape’ of this expression, note the analogy with the 
binomial theorem. ] 
Noration. By the expression 
δ Ὁ 0\? 
(μι τ), +kg) f(x,y, 2) 
we mean the result of operating p times in succession on the function 


2 Ὁ oe o 
f(x, y,z) by the operator ha +) By - κε: for example, when p = 1,2 


we have respectively af ΩΝ ΕΝ τ 
Tey" 
weds peat ἘΣ a + Bh oe + Ohi oe a 
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In the term 


9.8 .2\" ᾿ 
μος, τι) f(+Oh, y +0j,2-+k), 


it is understood that, arrmr the ἢ operations, the variables 2, y, z 
are replaced by «+ 0h, y+0j, z+ 6k. 

Suppose that x, y, z have definitely assigned values. Form the 
function F(t) of the variable ¢ by means of the relation 


F(t)=f(x+th,y+tj,z+tk). 
If, for convenience, we write 
w=at+th, v=yt+t, w=2z+ik, 
then (x, y, z being temporarily constant) 


δι du, f dv, of dw 
dt Oudt ᾿ dwdt ᾿ dw dt 


f oF, oF 
(Compare also p, 48.) Thus the operation 
4 
dt 
on the function F(é) is equivalent to the operation 
Ee hs ὃ 
on the function f(u,v,w). Applying the same rule to nf , which is 
a function of ¢, and to its alternative expression ἢ = + Jay ὦ + rl J 


which is a function of u, v, w, we obtain the equivalence 


d? 0 0 0λ3 
da (actiggthze) » 
- and so on. 


Now, by Maclaurin’s theorem (Vol. IT, p. 49), @ exists (0<0< 1) 
such that 


F(t) = F(0)+tF’(0)+...+ 


as re (0) +5 he (0t), 
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and so, using the relations just obtained, 


f (u,v, w) 
= f(u, υ, W)o 
+i(ho tig ἢ f(U, 0, W)o 


ἐπ- 0 . Ms) Qo \n-1 


(n—1)! 
ἐπ 0 0 
ἘΠ ΠΣ ΎΣ, ἢ f (u,v, w)e, 


where, on the right-hand side, the suffix 0 indicates replacement of ὁ 
by zero in f and its differential coefficients, while the suffix θὲ 
indicates replacement by θὲ, AFTER DIFFERENTIATION in all cases; 
thus 


0 
ΠΥ ἢ f(u, v, w)o 


ὃ .9 0 
means heey. 2) Ἐν fey?) + he ft.y2), 


rr. é 
and (ia tigthe sg) Su, W) a 


means («2 tigntks =) f(e+hbt, y+j0t, z+ kt), 


the values ἡ τ θέ, y+jOt, z+kOt being inserted in the partial 
differential coefficients arrmr differentiation by the current vari- 
ables x, y, 2. 
Putting ὁ = 1, we obtain the theorem. 
ILLusTRATION 19. Z'aylor’s theorem for the function xyz. If 
f(x,y, 2) =ayz2, 


δ᾽, σι ἢ 


and (r= g tigtks =) f(x, y, 2) = 2jka + 2khy + 2hjz. 
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Hence 
(i) a number @ exists such that 


(2 +h) (y+j) (z+k) 
= TYzZz 
+ hyz+jzat key 
+ ${2jk(a + Oh) + 2kh(y + Oj) + 2hj(z+ Ok)}, 
so that, on expanding and cancelling like terms from each side, 


hjk = 3hjk0, 
and θ-- ἰ; 
(ii) a number ¢ exists such that 


(x+h) (y+j)(2+hk) 
= xyz + h(y + 47) (2+ bk) +j(2+ Pk) (e+ Gh) + k(x + dh) (y + 47), 
so that, on expanding and cancelling like terms from each side, 
jka+khy +hjz+hjk = 26( jka + khy + hjz) + 3¢*hjk. 
Taking the general case, in which h, 7, k are not zero, divide by hjk; 


then pegs 
3 a - ΤΌ τ) - i = 
ϑφ5- (29 -- 1) (+445) 1=0. 


It is an interesting exercise in elementary algebra, which we leave 
to the reader, to prove that, if the expression 


ned 
ΚΕ ΤῈ 


does not lie between —1, -- 2, then there is precisely one value of ¢ 
between 0, 1; but that, if the expression does lie between —1, —2, 
then there are two such values of ¢. In either case, the theorem is 
verified. 


REVISION EXAMPLES X 


‘Scholarship’ Level 
1. Prove that, if zis a function of u, v, where 
U= XY, V=H+Y 
(so that z may also be regarded as a function of x, y), then 
Oz 40% oe ΟἿΣ 


Bat = δῶ μου tye 
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2. If the independent variables x, y in the function f(x, y) are 
changed to ξ, ἡ, where 
f=x+y, 7=%-y, 
so that the function f(x, y) becomes g(£, 7), prove that 
eA a 
Ox ὃν = a8 07” 
Deduce, or prove otherwise, that, if the independent variables 
x, y are changed to uw, v, where 
a+y=log(u+v), x—y = log(u—»), 
so that the function f(z, y) becomes A(u, v), then 
of af Oh ch 
say 0 (Fafa 
3. (i) Prove that, in general, 
ὃ | ( 0 =) 
x—+y—)(e—--y—] f(a, 
( δ. Yq) \ Sy Yaa) SO) 
2 
is not equal to (“3 -- y?) τὴν wy (S- ἘΝ ‘ 
(ii) Prove that, if f is a function of ὦ and 


dam (a*+y°)tan1Z, 


sin et Ze rnd 

and (= tye) (eZ - ἢ = 2(2?+y?) [2 τὰ. 
4, If c= Uw, Y=ut+y, 

prove that caer ss ον ὃν 


ann πον tote 


If z is a function of z, 


022 
νῦν =” 
prove that 
Oz 86% ε΄ ΟΖ 
(u—v) (us5 +055) - (u — 08) oe (w+) (=-F) = Ὁ. 


[A permissible method of proving the second part of the question 
is to assume that z is of the form f(x) +9(y).] 
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5. If ξ = ay, Q=2U+Y, 

find Me when 7 is constant, and oe ὃν when é is constant. 
5 og on” Oy 
03 σἂν x ] 
x ance) = | ey)? ey 
6. If u = e*(xcosy—ysiny), 
O72 O*u 
prove that δ ἐπ = 0, 
O2e% O2e% 
(Gert apg) = lee ey. 


If w=e(xcosy—ysiny), v= e*(xsiny+ycosy), 
eu _ e*{(x+1) cosy—ysiny} 
᾿Ξ re er meee 


prove that ἘΣ Εν ; 


where 2 is considered as a function of the independent variables w, v. 


7. Variables x, y are defined in terms of u, v by the equations 


x=ccoshucosv, y=csinhusinv., 
Write down the values of 
δ, Ἐν ὅν, ὃν 
Ou’ dw’ dv’ dv’ 
and prove that 


ὃ 

νὰν delet y)ainhe COS , 
a2 αὐ" ΤῸ Ά 2 2 
x; (x+y) +aa(t+y) = 4c?(cosh* uw —cos*v), 


8. The function u(x, y) is defined by the formula 
u = e*(xcosy—ysiny). 
Find a function v(x, y) satisfying the conditions 
ὃν 8. ὃν _ aw 
On ν᾽ ὃν τ 
υ(0, 0) = 


Show that u+iv is expressible in the form oa: Fre), with z = x+ty, 


and find the function F. 


REVISION EXAMPLES X 
9. (i) Prove that the equation 
Of Of ΟΝ 2f 
Gat ΤῊΣ ταν yee 9 
— σεν 
18 satisfied by f = e+ y?+ 22" 
(ii) Given that 
2=e"—“"cos2uv, y= e’-*sin Quy, 
Ox ὃν oy Ox 
prove that tees Oe eae 
Prove that, if fis a function of x,y, then 


of - 2(ua- ey) E+ Buy +02), 


2 = — 2uy +02) 2 + 2(ua —vy) 
Hence prove that (2) +) +(Z) 
τ ὦ 


is the same for all functions f. 


oy 


10. Given that 
a=ccoshucosv, y=csinhu sin», 
where c is constant, find expressions for 
ὃν 
Cw’ dv’ 
Given that f is a function of u, v, prove that 


2st ce sinh (u— -in (εἰς, 


oy)’ 
of 095 0 03 
ate f = c?(cosh? wu — cos? v) (sat ae =a)" 
: U 
11, Given that z= eae y= ea 


and that f is a function of x, y, show that 


ΤῊΣ 1) 
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12. (i) Given that 
w=2+y, v= 2+y', 
prove that, if x is considered as a function of wu, v, 


On y Oat ] 
ὃι ὃὥχ(ω--ψ) Ov ϑ8χί(α--ν) 
Prove that Oe O02 
Ou ov l 
dy dy Gay(x—y)” 
du ov 


(ii) The variables x, y, are connected by the implicit relation 
x? + y?— 2aycosé = sin? 0. 
Prove that 0 = + (sin-+z—sin—'y), and deduce that 


sae Ee 
daoy 
13. If ne = SAE — 0080+ cosh g, 
Oc Oy ὃ oy 


prove that 3036’ δ δ᾽ 


The function u(x, y) transforms by means of the above relations 
into υ(θ, φ). Prove that δυλϑ (ϑυλϑ 
(Ὁ 0). (se) + (as) 
Oa cy)  (cos@+cosh ¢)?" 
14. The identical relation 
f(u?—2?, uw? —y?, w®—z*) = 0 
defines u as a function of x,y,z. Prove that 
Poy lew tow ἢ 
zon yoy 202 οἷ 
15. Prove that, if 
x=rsinOcos¢?, y=rsin@sind, z=rcosé, 
and if r is a given function of 0,¢, so that z may be regarded as a 
function of x,y, then 
or. ; a ὃν. ; 
(ssn 0+ rsin 90080) = + =r sin 4 +rsin* 0.089 


— 7p Sind cos @ 008 ὁ = 0. 
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16. Prove that, if 
z=uUut+v+w, y= vw+wutuw, z= uw, 


and if z is a function of x,y so that w is a function of u,v, then 


ow Ow 
ou J ov stu) - 
vw—p—qv+w) wu—p—qwt+u) w—p—g(ut+r)’ 
where Pili 1-5. 


17. Show that, if x,y are functions of ξ,η, and if 
2 


then, for any differentiable function U(z,y), 
oU dUdy dU cy 


en = GE δὴ δὴ 08 
ἢ) ὃὕδε OU Or 
Ay ~~ OE ant Oy OE 
Prove also that 27] 80) 417 517 
eu eu 1 a Pae—Magh 1 a [Maz στη 
“tts π᾿ ΠῚ πὶ 3” 


18. A function f(#,y), when expressed in terms of the new 
variables u,v defined by the equations 
x=t(ut+v), y=w, 
becomes g(u,v). Prove that 
ΓΕ 422 L215) 
νυ 4ιλοχξ “ydacy οὐ yoy)” 
19. Prove that, if f(x,y) is a function of z*+y? only, it satisfies 
the identical relation 
δε. (Oy 
By changing to polar coordinates, or otherwise, prove con- 
versely that, if f(z, y) satisfies this relation identically, then it must 
be a function of x? + y* only. 
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20. Prove that, if x,y,z are functions of two variables u,v given 
by the relations 


ῷ τὸ [(,0}, y=g(u,v), z= h(u, υ), 
then (suffixes denoting partial differentiation) 


(fude ~ fu) = My Jo— hous 
, fia fo Suu9o— AFuvIuIo + fovFs 
(fuIe—So9u) δα =| Iu Go θωωιθὺ-- 2GuvIuJv t+ IuvGu 1. 
Ping hy ἤμυθυ-- ἢ, JuFy t+ yw Ge 
21. The function y of two variables 0, z is defined implicitly by 
the equation y = 0+2d(y). 
If u is a function of y and F(u) a function of u, prove that 


sa \F() gel = a (Fw) so) 


ou ou uw a Ou 
and that Ox = PY) 59’ 2.3 ‘a ΕΣ 


and, generally, that μὴν 
gn Cu 
aan = ὅττι WZ | 


22. Prove that, if 
x=rcos0, y=rsin#, 


and if ¢ is any function of r, 8, then 


op δῴ o¢ sin 0 
ba Op eo πτ' 
and obtain a corresponding expression for a 
Prove that, if g=rsinn8, 
ed 92 
then tata = 0, 


23. Prove that, if 
| U = f(z/y), U, = 7"U, 
oU 0 
ΞΡ Εις θυ 
where r? = 7?+y?, then oma ty 0, 
2 2 8 


oat + Oya * 78 
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Hence, or otherwise, prove that, if the equation 


is satisfied by τ = V(x, y), where V(x, y) is a homogeneous function 
of degree n, then the equation is satisfied also by τ = r-*"V (a, y). 


24. Four variables u,, U2, Us, U, are such that any three are 
independent and each may be expressed as a function of the other 
three. If u,. denotes a , the rate of change of u, with respect to uw, 

2 
when 1,5, u, are kept fixed, and if εἰς; (¢+7) has a similar meaning, 


show that 
Ugg Uzi Uyg = -- ἢ}, Uya Ugg = 414 1055. 
25. The variables 2, y, z satisfy the two equations 


f(a, y, 2) = 0, g(%,Y, z) = 0. 


By eliminating z between these equations, it is possible to obtain 
a relation connecting x, y which defines y as a function of x. Show 


that 
dy aa . 9.- 9. 
da τ ν9ε -- 2Gy 
26. Given that x, y are functions of u, v defined by 
F(x, ψ 0, v) ἊΣ 0, g(x, ¥,U, v) = 0, 


find = in terms of partial derivatives of f, g with respect to 2, y, u, v. 


If x*+y?—25uv = 0, 
uz+vy—1=0, 
prove that = = τ τ when u = v = 1, and give the reason for the 
ambiguity in sign. 


27. If f(x,y) is a function of x, y, where x, y are functions of t 
defined by the relations 


u(z,t)=0, v(y,t) = 0, 


and if f(x,y), when expressed as a function of ἐ, takes the form 
g(t), prove that 


dt 


dg Of Ovdu of dudv) Jou ov 
- τὸ δεδι δα oy" 
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28. (i) If U = xyz, 
where x, y, z are connected by the relations 
χει τα, x+y+z=56 (a,b constants), 


du 


prove that ot (α --- ψ) (ὦ --- 2). 


(ii) If £, ἡ are functions of x, y such that 
E=e*cosy, ἢ = e*siny, 
and if x, y are functions of r, Θ such that 
x=e'cos#, y=e'sin#§, 


where r is a function of 0, prove that 


dr 
de 9 8a y+ 9) 
dy . adr 


1+, tan (y+0) 


29. The equation z= F(x, y) 
is obtained by eliminating u between the equations 


υ = f(u, “), = g(u, 2). 
afar af ay af ὃ 


Pro ee eee ae eee 
—— Bu On Oude δόδω 
of oF ὃ 
Ou Oy δὲ 


30. Prove that, if wis a function of x, y which is transformed into 
a function of r, Θ by the relations 


2? r cos 40, y= rt sin 40, 


Ou = ind 
Cu Cu ou 

2 πὸ" a 
(ii) pia " δ’ 


a ae oe ae 
(ui) a(rsat +s +a) = Bat ha 
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31. Giventhat w=rcos0, ᾧ τὸ γβίηθ, 
06 08 080 
find On’ oy’ ox Oy’ 
and prove that 0970 (2η-- 1)} 


sin (2n0 —inz). 


ΣΙΩΝ gen 
32. If, y are connected with r, # by the relations 
χ 0050 --ἶἰἣν βίη θ᾽ = xsind+ycosd =7, 
88 ὃ 00 OF 
find ax’ Ox’ oy’ oy 
in terms of r, @. 
Show that, if V is a function of x, y, 
eV eV 1(feV 10V 10 
Sait Bye 3 [δα ror TOI" 
33. The variables 2, y are changed to 7,0, where 


z=rsec0, y=rtand. 
Show that 
O24 O2u d%% 1du cos?Ad% sin#@ cos@ du 


it yt OF τὸς τὸ OCC 
34. aes that w,v are functions of the two variables z,y, 
express (sat +) (uv)in terms of the partial differential coefficients 


of uw and v. 
If w is a function, οὐ Ὁ third-order differential coefficients, satis- 


Ou : 
fying the relation Stes τ = 0, and if 


w = (a27+y?+an+by+c)u, 

m qs 0g \2 ; 

ow that (ata) w=, 

35. Four variables ὦ, y, 7, θ are connected by the two relations 

2=rcos?, y=rsind. 

A symbol like (=), indicates that r is expressed as a function of x, y 
and that y is kept constant during the differentiation of r with 
respect to x. Verify the relations 


0 (F),@),-> ὦ G),-@), 


aw (ἢ, -- ἢ, = (@),- "(a 
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Illustrate (iii), (iv) geometrically, taking (x, y), (7, @) as Cartesian 
and polar coordinates of a point in a plane. 


36. If A is the area of a triangle ABC, calculate θΔ when the 


θα 
independent variables are 


(i) a,6,C; (ii) a, B,C; (iii) a,b,c. 
Illustrate your results geometrically. 


37. If the circumradius R, and the area A, of a triangle ABC 


are regarded as functions of ὃ, c, A, prove that 


oR oR OR 

Neda τὰ 

ORCA  eROA 

db Gc” Gc Ob 

38. If (x,y) are the Cartesian coordinates of a point and (r, 0) 


its polar coordinates, define what you mean by “4 ᾿ς 


= }Rsin A. 


Find their values, and illustrate by a diagram, 
apo ell 
θαδὸν λόχον 
39. If ὦ τῷ [(0}, where 
ὥ τε γοοβθ, ᾧ -ὸΖ γβίηθ, 
Ou Ou df idf 


Prove that 


prove that aa +5 > ag ue 
40. If z= 2'f(y/x), 
prove that en ας = 12, 
nh 
«ΟΣ οἷ 


waa t ἐπ ἐν τὰ Ξε r(r—1)z. 


41. If u = f(ax* + 2hay + ὃν), 
Ὁ = g(ax? + 2hay + by*), 


0 ( a ὃ ( ov 
prove that By (u =) =a (u =) . 
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42. A function f of two independent variables r,@ is trans- 
formed into a function g of variables w,s by means of the relations 


roos@=I1/u, tan? - 8. 


Prove that rf a | F = us + (1 +62), 


du’ 06 ds 
Prove also that, if f satisfies the equation 
of 
008 0F + rein 0% = = 0, 
then Y (149%) g(u), 


where ¢ is some function of w. 


43. If x,y are the coordinates of any point and 


C= rcos 86, ¥y¥= rsin 0, 


fol ive δὲ (5) ᾿ (=) : 
06 r const. Ox y const. 


Verify the second result geometrically. 


020 5530 
Prove that aa + By? = 0, 


44, The variables x, y, z are connected by the relation 
x+y? +22 3xyz = 0, 
and P(x, y, 2) Ξε σ᾽ γῆς. 
Determine the value of = oF at (1,1, 1), 


(i) when the nero <i variables are x, y, 
(ii) when they are 2, z, 
and explain geometrically the difference between the meanings of 


0 
τ in the two cases. 


45, Each of wu, v, w is a function of the three variables z, y, z. 
Whenever 2, y receive small increments which satisfy the equation 
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dy—zdx = 0, the corresponding increments du, dv satisfy the 
equation dv—wdu = 0. Prove that 


wae ταί we) = 0 
Ὁ θυ Cy = Gy : 
ov = Ou 0. 

a Oz 


Find the most general functions τι, v, w for which u is a function 
of x only, v of y only, and w of z only. 
46. If d(x, y) is a differentiable function of x, y, and if 
sin u f(v) 
~cosu+eoshy’ %~ cosu+cosho’ 
show that, for a certain form of the function f(v), 


not tux = sin u cosh v2 + cosw sinh oS, 


and find the nae of f(v) in this case, 


41. Ifu=2/y, v = xy, prove that 
2 2 
αὐ εν 20F _ (w 52 OF yy 1, 


δ ov" 
22-0, (om Zh Z), 


48. The area A of a triangle is found from measurements of the 
side a and the angles B, C. Prove that the error 6A in the calculated 
value of the area due to small errors da, dB, 6C is given approxi- 


on e ὃδλ δα ὁ ὃΒ ὃ 80 
Δ “a asinB’ δα" 
49. The side 6 of a triangle ABC is calculated from measure- 
ments of the side a and the angles B, C. There may be an error not 
exceeding ἃ (in either direction) in the measurement of a and an 


error not exceeding ἃ in either or both of B, C, where ἢ, « are small. 
Show that, if B+C <4m, the greatest error in ὃ is approximately 


[ +acot B}. 
a 


Find the greatest possible error in ὃ when B+C > 47, and explain 


the difference in the two cases. 
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50. If f(x, y) =axy + bx+cy+d, 
g(x,y) =a'ay+b'x+e'y+d’, 
ὃ (ἢ ply) 
prove that aH = gt ; 


where p(y) is a quadratic function of y (not containing 2). 

Hence or otherwise prove that f(x, y)/9(x, y) is a function of the 
product XY, where X is an appropriate function of x only and Y 
a function of y only. 


51. A function y of x is defined by the equation f(x,y) = 


Express οἰ ταὶ in terms of p, g, 7, 8, t, where 
of of 
P= Ox?’ πὸ Cy? 
of y ae of 


sia sla Oxcy  dycx’ oy?" 
If the same equation is regarded as defining 2 as a function 
of y, prove that dx dy 


52. Show that, if x = pcos¢, y = psing, then 
ΟΡ eV eV 10} 105} 7 
Ga? * Oy? ρὲ p Op ρὲ δῴ"᾽ 
Hence show that ‘ 
tant! ztan-12, (p? +22)-3 tan“ 
are all solutions of the equation 
CV eV &V 


at os te δὰ 
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CHAPTER XV 
MAXIMA AND MINIMA 


1. The general conditions. Suppose that 


u=f(x, y, 2) 

is a function of the three variables x, y, z (not necessarily indepen- 
dent). The function u is said to have a Maximum at a point P, 
for which x =a, y=b, z =c, if its value at P exceeds that at 
neighbouring points, and a mrnmmum if the value is less. As for 
functions of one variable, maxima and minima are local properties. 

The argument for several variables is naturally more complicated 
than for one only, and the treatment which follows is of necessity 
more sketchy. 

We begin by finding necessary conditions, assuming first that the 
three variables x, y, z are independent. 

Suppose that the function 


u=f(x, y, z) 
has a maximum value when 2 = a, y = b, z = c. Its value is then 
greater than that at any near point; in particular, at all those near 
points for which y = ὃ, z = ο. Thus the function 
f(x, 6, ¢) 


of the single variable x is greater for x = a than for any other near 
value, so that this function has a maximum atx = a. Its differential 
coefficient with respect to x (while y, z remain constant at ὦ, ὃ 
respectively) is therefore zero when x = a. Hence 


af 
Aa 


when z = a, y = ὃ, 5 =c. Similarly 


Suite arg 


oy ᾿ ὃ: 
for those values. 


The presence of a minimum may also be treated in the same way. 
Hence a set of necessary conditions for the function 


u=f(x,y,z) (x,y,z independent) 
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to have a maximum or a minimum value for x = a, y = ὃ, z = cis that 


forw=a,y=b,z=Cc. 

Since these conditions must hold at each maximum or minimum, 
they enable us to locate such turning points. But the conditions are 
not sufficient; moreover, they do not, of course, distinguish between 
maxima and minima. Detailed analysis is beyond the scope of 
this book, but common sense will often provide the answer for 
particular problems. 

The three conditions may be gathered into a single differential 
form. For the differential of the function a is 


du= δ ἀνα dy + Zee, 


which vanishes at a maximum or minimum. Hence the differential 
ued das Lay εἶ dz vanishes when u has a maximum or mint- 

mum value, for arbitrary values of the differentials dx, dy, dz. 
Expressed in this form, the condition du = 0 for a maximum 

or a minimum may be applied whether x, y, z are independent or not. 

For if they are not independent, the conditions of dependence 

enable us to reduce the number of variables until those which 

remain are independent, and the condition du = 0 then follows. 
Note. If x, y, z are not independent, the condition 


du = 0 
for a maximum or minimum is still equivalent to 
Ou ou Ou 


dz 
but it does not follow now that = a a vanish separately. For 
the treatment to be adopted, see § 3. 
InLustRATION 1. Τὸ find the greatest distance of a point onthe surface 
ax? + by? + cz? + 2fyz+ 2gza+ 2hay = 1 
(assumed to be an ‘ ellipsoid’) from the plane z = 0. 
(The surface may be pictured as a somewhat distorted sphere 


with its centre at the origin.) 
6-2 


76 MAXIMA AND MINIMA 


The distance of the point (#, y,z) on the surface from the plane 
z = 0 is simply the coordinate z, and, at a point of maximum (or 


minimum) distance the two partial differential coefficients ss ; ω. 
Ox’ dy 

must vanish. 
Differentiate the equation of the surface partially with respect to 


x, and then put = = 0. Thus 


ax+hy+gz = 0. 
Similarly, hz +by+fz = 0. 

These two equations, together with the equation of the surface, 
determine the points of greatest distance; the z coordinates give 
the actual distances. 

In order to find z, we eliminate x, y as follows: 

Multiply the two equations of condition by 2, y respectively, and 
subtract from the equation of the surface; divide the resulting 
equation by z. Thus 1 . 
gx + fy + cz— “in 0. 


Eliminating z:y:1 determinantally between the three linear 
equations, we have ey 


gz 
h ὃ fz = 0, 
1 
5 1 ἐς-- 
or, after division by z, 
ah g eon 6 
h b 7 1 h b ᾿ ΠΝ 
g¢ J pile 2 
Expanding, 
(abe + 2fgh— af? —bg®—ch®) — (ab — ht) = 
ab —h? 
, h τς oe Ae... | Ὁ. 
wae a | aera appa): 


From the shape of the surface, this gives the greatest distance, 
whose value is therefore 
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2.* To distinguish between maxima and minima. Though 
detailed consideration would carry us much too far, the general 
picture may be exhibited of a method to distinguish systematically 
between maxima and minima. 

For, say, three independent variables, suppose that the function 


+= f(x,y; 2) 
has a turning value at the point (a,b,c). It will be a convenience 
of notation to express the conditions in the form 


ou ou Ou 
meet π = 9 a 


If ἕξ, ἡ, ζ are small, the value of the function at the near point 
(a+f,b+7,c+) is flaré.b+7,c+0, 
where, by Taylor’s theorem (p. 57), 

f(a+&,b+,c¢+)—f(a, 6, c) 


= (ἐξ τς, εξ few 


ἐξ (εξ τη tse) tows) 


Hives 


the differential coefficients being evaluated at «=a, y = ὃ, z=C. 
Taking £, 7, ζ to be so small that powers and products of degree 
greater than 2 may be neglected, and remembering that 


au _ du _u_y 


Ga 0b Ce 


we obtain a relation which we may write in the form 
f(a+é,b+9,¢+€) —f(a, 6,¢) 


= (eae + a 63 ἜΣ": ΡΣ a +2 + 2aaapen\. 
where £, ἡ, € are variables, depending on the point selected near to 
(a,b,c), and where the coefficients a ... are constants. 


* This paragraph may be postponed, if desired. 
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For a MAxImMuM this expression must be negative for all sets of 
values of ξ, ἡ, €; for a MmNIMUM it must be positive. This condition 
is both necessary and sufficient. If, however, the expression varies 
in sign, being positive for some values of ξ, ἡ, ¢ and negative for 
others, the function uw has neither maximum nor minimum. 


The simpler cases of genuine maxima or minima can be settled 
by definite algebraic tests: 
(i) For convenience, write the expression with the notation 
B= Af? + By? + Cl? + 2F 46+ 2GCE + 2HEn, 


and form the determinant whose rows are the coefficients of 2, 27, 
oH oH oH 


2¢ in the partial differential coefficien 2’ Oy’ 2 respectively, 
namely, te» 5. Ὁ 
Hz δ, δὶ. 
Gre 
Write down also the ‘leading diagonal’ determinants of 2, 1 rows, 
βοδιῶν, AH|, A. 
Ἢ 8 


Then it is known that a necessary and sufficient condition for E to 
be always positive, for all values of £, ἡ, €, is that these three deter- 
minants should all be positive. 

Consider, for example, the expression 


3G? + 2η5 + 26? + 2nf — 4ζξ — 4ξη. 
The determinants are 
3-2 -2 


οἷ, ae a fen 


3-2 
= 2, 
—2 2 


and 3. 


Since they are all positive, the expression is positive for all , 7, ¢. 
(ii) The test for αὶ to be always negative is that — should be 
always positive. 
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Thus the function f(x,y) of two independent variables has a 
MAXIMUM at the point (a, b) if 


<8 

ssia— (soap) > 9% 

andamenouif 2-0, Zao, 
ct >0, 

sete (Bat) > 


(Note the direction of the last inequality.) 
(iii) The method can be extended, in an obvious way, to any 
number of variables. For the expression 


Ag? + By? + C62 + Dr? + 2F yf + 2466 + 2HEn+ 2UEr+ 2Vyr+ 2Wer, 


the determinants are 


. 8. ΤᾺ. 5. SN: 
ἐδ στα ἔτλη al meni poor 
GFCOW G FO 
a I 


3. Lagrange’s method of undetermined multipliers. The 
method now to be explained may be grasped more readily if we 
begin with an illustration: 

ILLUSTRATION 2.* Τὸ find the lengths of the axes of the conicin which 
the quadrie ax? + by? + cz* + 2fyz+ 2gza+ 2hay = 1 
meets the plane la+my+nz = 0. 

* The reader who has not yet reached the theory of quadrics may simply 
regard the problem as the determination of the maximum and minimum values 


of the faction a*+y%+z*, where the variables z, y, z are subject to the two 
equations quoted in the enunciation. 
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The squares of the lengths of the axes are the maxima and minima 
of the function 


U=2?+y?+2%, 

At such points, du = 0, so that 

xdz+ydy+zdz = 0. 
Since the point (x,y,z) lies on the quadric, the differentials satisfy 
the equation 

(ax + hy + gz) du + (ha + by + fz) dy + (gx +fy +z) dz = 0; 

and, since it lies on the plane, 

ldx+mdy+ndz = 0. 


Multiply these three equations among the differentials by 1,A,z 
respectively, and add. Then choose A, yw so that the coefficients of 
da, dy both vanish; that is, so that 

w+ A(ax+ hy +gz)+ ul = 0, 
y +A(ha + by +fz)+ ym = 0. 
Then [2- λίρα -ἰ [Κ-ἴ ο2)-Ὁ μΉ] ἀξ = 0. 

Now the three variables x, y, z are subject to two conditions, 

namely, the equation of the quadric and the equation of the plane; 


their three differentials are therefore subject to two linear con- 
ditions, so that one differential may be given an arbitrarily assigned 


value. In particular, dz may be assigned arbitrarily—it need ποῦ 


be zero. Hence the coefficient of dz in the equation last written must 


be zero, so that z+A(gu+fy +cz)+pyn = 0. 


To summarize, the five equations 
ax® + by* + cz* + 2fyz + 2gza+2hay = 1, 
le+my+nz = 0, 
x+A(ax+hy+gz)+pl =0, 
y+A(ha+ by+fz)+pum = 0, 
z+A(gu+fy+cz)+un τὸ 


serve to determine the five unknowns consisting of (i) the three 
variables x, y, z, evaluated at a turning point, (ii) the two multipliers 
A, μ. The sacrifice involved in increasing the number of unknowns 
is amply repaid in symmetry. 
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In problems of this type, where the equations of condition are 
homogeneous polynomials in 2, y, z, a good first move is to multiply 
the equations involving A, 4 by «, y, z respectively, and then to add 
them. This gives 


(x? + y? +27) 
+ Aaa? + by? + cz* + 2fyz + 2gza + Zhay) 
+ p(la+ my +nz) = 0, 


so that u+A(1)+ (0) = 0, 

or u+A=0, 

or A=-4. 

Hence x—ulax+hy+gz)+pl = 0, 
or (2—)e+hy+ge—f1 = 0. 
Similasiy, has (b-=) y+ fe—Em = 0, 


Ι | 
gx+fyt+ (c-<) ‘Ts = 0. 


Also le+my+nz = 0. 


Eliminate the ratios x:y:z:—y/u between these four equations. 
Then 


α-- ἃ gi ἢ 
1 
1 
g f 6m = τὰ 


This is, on expansion, a quadratic equation in =, and therefore in wu, 
whose roots determine the maximum and minimum values of τι. 
The lengths of the axes may thus be found. 


The next illustration will help to consolidate what has been done, 
and also to expand certain details. 
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ILLUSTRATION 3. Τὸ find the maxima and minima, for NON-ZERO 
values of x, y, 2, of the function 

U=xy2, 

under the condition e+y+z=1., 

At a turning value, du = 0, so that 

yzdx + zxdy+axydz = 0. 
By the equation of condition, 
dx +dy+dz= 0. 


Multiply these two equations among the differentials by 1, A 
respectively, and add. Then choose A so that the coefficient of dx 


vanishes; that is, so that gz-+A = 0. 


Then [za +A] dy + [zy +A] dz = 0. 


Now the three variables x, y, z are subject to one condition; their 
three differentials are therefore subject to one linear condition, so 
that two differentials may be given any arbitrarily assigned values. 
In particular, dy, dz may be assigned arbitrarily. Hence the 
coeflicients of dy, dz in the equation last written must both be zero, 


so that zx+A= 0, 
xy+A= 0, 

To summarize, the four equations 
e+y+z=1, 
yz+A=0, 
χα-ξ λα Ξε 0, 
xy+A=0, 


determine the four unknowns 2, y, z, A. 
Multiply the last three equations by x, y, z, and add. Then 


3xyz+A(z+y+z) = 0, 
or 38u+A = 0, 
or A =—3u. 
Thus yz = 3u, 

| 22 = Su, 


xy = 3u. 
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Multiply corresponding sides of these three equations; then 
x*y*z? = 274%, 
or 3 = 27ι.. 


Now w is not zero, under the conditions of the problem (2, y, z 
non-zero). Hence eal 


We may extend this illustration to show how to settle whether 
such points are maxima or minima. 
It is easy to prove that there is only one such point, namely, 


S=y=z= t, 
Consider, then, a near point 
S=ttp, γα, z= $+", 
where p, g, r are small. Since x+y+z = 1, we have 


ptqtr=0. 
Now U = xyz 


= ($+) ($+9)($+7) 
= th(ptqtr)+hart+rp +pq)+pqr. 


Suppose that p, 4, r are so small that the product pgr may be 
neglected. Then 


w= ἷγ τ ἐ{ξ( Ἐφ ν7)}-- ἐ(»5Ὲ 45 +2°)} 
=g7-a(P tg +"). 
Hence w has its locally greatest value at (4, 4, 4), and the value ἦγ is 
therefore a maximum, 


We now apply these ideas more generally. 
To locate the maxima or minima of the function 


& =f (X1, Ve, ...»,Ζ6,} 
of the n variables 21, %2,...,%, subject to the m relations (m <n) 
P3(Xy, eg, «.-,%_) = 0, 


P(X, Xo; ping Reg) = 0,. 


P(X; Boy woes 4.) = 0. 
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At a turning value, du = 0, so that 


SRS of 
le +7 a Akg + 00 b= in. dz, = 0. 
By the equations of condition, 


Multiply these m+1 equations among the differentials by 
1, Ay, Ag, ..-,A,, respectively, and add. Then choose the multipliers 
A in such a way that the coefficients of dx,,dz2, ...,dx,, all vanish; 
that is, so that af 


haa site. oes ee 


. et. soph dag ΣΕ Ὁ, 


On oz, 
Then Ee + λεῖα aren 


Now the n variables ,,...,%,, are idee sai m conditions; their 
n differentials are therefore subject to m linear conditions, so that 
n—m differentials may be given arbitrarily assigned values. In 
particular, dz,,..1,0%m19, -..,d%, may be assigned arbitrarily. Hence 
the coefficients of dz,,,,,...,dz, in the equation last written must 
be zero, so that 


i ΕΝ tag et = 0 


In all, then, we have the m+n οὐδε 


D(24, Le, .2-%,) = 90 (§ = 1,...,m), 


fe tage + dag 0 ({=1,...,n) 


REVISION EXAMPLES XI 85 


to determine the m+n unknowns 2},...,%,, A,,...,A,,. This can, 
in general, be effected, and the value of τὸ ascertained. 

Note. The method breaks down if all the determinants formed by 
taking m columns of the ‘matrix’ 


a ae op 
θα, > C25 > 3 OL y > 
Os. Os Ops 
θα, > O25 > ? 2 Bs 
δῴ, δῴ, 0m 
oe,” Ox,” OL,” 


should vanish. The equations of condition are not then independent. 
But the reader is unlikely to be troubled by this case, at any rate for 
the present. 
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University Level 
1. Find the maximum and minimum values of 


r+y—-1 
a? + Qy?+ 2° 


p\* (ἀᾺλῈ (μὰν 
2. Ifz= (?) + (=) +(%) , where p,q,k are constants, prove 
that z has a stationary value when Ὁ, ὦ, y, g are in geometrical 


progression. 
3. The variable z is determined as a function of x, y by the 
(ae a3 + y%+23—6ayz+3 = 0. 


ΟἿΣ ΟἿΣ Ο2 
a2’ dy?’ dady 
that these values of x, y, z make the function 


z+ 2. 


Find the values of — when x = y =z = 1, and prove 


& minimum. 


4, Find the minimum value of 
x+y? + (ax+by+c)*. 
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5. Find all the stationary values of the function 


a Sat Li 
Y 2 a 


of the variables z, y; discuss how their nature depends on the value 


of the parameter a. Consider in particular the points at which 


ὥ τεῦ τὸ +a. 


6. Find the greatest volume of a rectangular parallelepiped 


which can be placed inside an ellipsoid of semi-axes a, ὃ, ὁ with its 
edges parallel to the axes. 
[The edges 2x, 2y, 2z of the box are subject to the condition 


gw? y2 22 
a sR a 1.] 


7. Find the maximum and minimum values, for real values of 


2%, y, 2, of the quantity ἡ + y¥*+2z* subject to the conditions that 
lx+my+nz = 0, 7 
ax* + by? +c2* = 1, 


where a, b, ὁ are positive and 1, m, n are real. 
Verify that the values determined are real and positive. 


8. Find the maximum and minimum values of the function 
U=xe+yi+z23, 
where x, y, z are connected by the relations 


r+y+z=a, 
αϑ y+ 22 = a", 
9. If ax+by+cez=1 (a,b,c positive) 
show that the values of x, y, z for which 
ΒΞ £2 
pam 4. = +f - 
oy s 
is stationary are given by 
| ax® = by? = c2z®, 


Show that this is a true maximum or minimum if xyz> 0, 


10. Use Lagrange’s method of undetermined multipliers to 
show that the triangle of maximum area which can be inscribed in 
a given circle is equilateral. 
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11. The sum of the twelve edges of a rectangular block is k; the 
sum of the areas of the six faces is },k*. Prove that, when the excess 
of the volume of the block over that of a cube, whose edge is equal 
to the least edge of the block, is greatest, the least edge is εξ; 
and find the other edges. 


12. Find the greatest and least distances from the origin of a 
point on the surface 


(or 


where α, b, c are fixed positive numbers and » is a fixed even integer 
greater than 2. 
[The square of the distance is x + y? + z?.] 


13. Prove that, if a, β, y are positive, there exists a shape of 
triangle for which the maximum value of 


sin’ A sin’ B sin’ Οἱ 
is attained; and that for this triangle 


tan? dA = (6 ἘΡβΈΥ) porep i βία Ὁ Ἐ:Ὲ}) poneg — 6: 86: 
ἬΝ τ ΒΗ αβ 


14, Find the stationary values of 
y* + 453 -- 4yz— χα — 2ay 
subject to . 22° + 3y? + 62? = 1. 
15. Show that, if x, y, z are connected by the relations 
x+y?+22 = 1, 
la+my+nz = 0, 
and a, ὃ, c are not all equal, then the extreme values of the function 
V =ax* + by? + c2z* 
are the roots of the equation 


I? m? η3 


16. On the surface given by 
x? + y? +22 — 2u + 2y+62+9 = 0, 


find the stationary points of the function 2+ y*+2?—yz—zx—ay, 
and investigate their nature. 
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17, The variables x, y, z are connected by the relations 


d(x,y,2)=0, (x,y,z) = 0. 


Prove that a necessary condition for f(x, y,z) to have a stationary 
value at (2%, Yo, 20) is that 


f f of 


Ox Oy 
op op 
ὃς by 
of oY oY 
On oy 
should vanish at (0; 20. 20), where f, ¢, wv are assumed to be 
differentiable. 
Prove that, if e+y+e=—xzyz=1, 


then there are three points (%,y,z) at which 2*+y*+2? has a 
stationary value, and show that they are minima. 


WML we 


2 


18. If f(x, y) =x? — 3ay? + 18y, 
where 32*y -- νὃ — 62 = 0, 
prove that the values of x, y which make f(x,y) a maximum or 
minimum are #=y= + /3. 
19, Investigate the character at the points (—5, 4,4), (1, 1,1) 
of the function z? + y* +z*, where 
e+y+z2=3, 2iyie=3. 


20. Determine the minimum value of 
ax* + by* + cz*, 
where a, ὃ, c are positive constants and the variables are restricted 
by the relationship x+y+z = 1. 
Hence, or otherwise, determine the minimum value of 
px? + gy? +12" + 2yz + 22x + 2ry 


subject to the same restriction, where Ὁ, 4,7 are constants, greater 
than unity. 
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21. Show that the minimum value of 
(aac? + by? + 022?) /a*y?22, 
where az? + by? + cz? = 1, with a, ὃ, c positive, is given by 
= aauray Yuet)“ πωτα' 
where τ is the positive root of the equation 
ει — (be +ca+ab)u—2abe = 0. 
22. If r denotes the distance from the origin to a point on the 
curve in which the plane ἴα + my +nz = 0 meets the surface 
(2? + y? + 2)? = a7? + δὲν + 0222, 
prove that the non-zero maximum and minimum values of r? are 
roots of the equation 
2 γι n* 
ant R_At aA 
23. Ifm,n, p are given positive numbers, find positive numbers 


2“, Ψ, 2, Whose sum is a constant A, such that the function 2™y"z 
is a Maximum. 


0. 


CHAPTER XVI 
JACOBIANS 


1, Introductory example. Suppose that two variables 2, y 
are expressed as functions of uw, v in the form 


x= u*+2v, 
y=ut+». 
Then w, v can also be expressed in terms of x, y; for 
x—2y = u?—2u, 
so that (u—1)*? = #—2y+1, 
or w=1+,(~—2y+1), 
and so v=y—1F /(w—2y+1). 


The given relations express x, y uniquely as functions of u, v; but 
there are two sets of relations to express τι, v as functions of z, y, 


namely, EITHER = 1+ ἡ 2941), 
v=y—1—j(e—2y+1), 

oR, alternatively, [ = 1 -- (ὦ -- 2. - 1), 
v=y—1+/(e—2y+1). 

If, as is natural, we now make the assumption that w, v are to be 
continuous, then in any particular problem we must keep to one or 
other of these two alternatives and not move between them. For 
example, if we choose the First pair of expressions when z = 6, 
y = 3, we have foc εν tte - ἢ 

v = 38--] --,,(δ-- 6- 1) =1; 
and if we then proceeded to take the srconp for the near values 
x = 6-21, y = 3, we should have 


tarde ane Rtn = ]-- 1.1 -- -- ΟἹ, 
v= 3-—1+,/(621-64+1)=2411= 8:1. 
The pairs u=2 u=-01 
pao 
=] v= 31 
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clearly reveal a discontinuity from which the pairs 
= 2 το 2. 
c auf ἃ 1 
ὕῦ-ι 1 υ τ 09 
(formed by keeping to the rrrst expressions) would not suffer. 


We can therefore express u, v uniquely as continuous functions 
of x, y by using EITHER the relations 


ὦ Ξ- 11 γ(α --ὰῶν Ὁ 1), 
Biers mde 

ok the relations u=1—/j(x—2y+1), 
ὰ = y—1+/(x—2y+1). 


There is, however, one condition under which this uniqueness of 
expression breaks down. If the values 2, y are chosen so that 


x—2y+1=0, 

we obtain without ambiguity the solution 

w=1, v=y-1l1 
and we can, for values of x, y near to those so chosen, take either 
of the alternative forms of expression while retaining continuity. 
For example, if x = 7, y = 4, then either of the pairs 

u=1+,/(7—8+1) u = 1—,(7-—8+1) 
ἢ = 4--Ἰ -- Ἰ( -- 8-Ὲς 1) ‘i {, = 4-—1+,/(7-—8+1) 

gives | u=1, v=3. 


Moreover, the near values x = 6-9801, y = 3-99 give, for the first 
alternative, 


u = 1+ ./(6-9801 — 7-98 + 1) = 1+./(0-0001) = 1-01, 
{, = 3-99—1—./(6-9801—7-98+1) = 2:99—0-01 = 2-98, 
and, for the second, 
. =1-0-01 =0-99, 
v = 2.99- 00] = ὃ. 
Each of the alternatives (1-01, 2-98) and (0-99,3) maintains con- 
tinuity with the pair (1, 3). 


We may give a geometrical formulation of what we have just 
indicated, taking x, y as the coordinates of a point in a plane 
(fig. 113) in which the line x—2y+1 = 0 is drawn. 


7-2 


92 JACOBIANS 


Consider, say, the solution 
{* = 1- {(ἡ --ῶν- 1), 
υ =y—1—,/(v—2y+1). 
When we have decided on these formulae, the expressions for u, v 
in terms of x, y are unique and continuous. It is implicit, though, 
that x—2y+1 is to be positive, so that the point (x, y) always lies 


Fig. 113. 


‘below’ the line x—2y+1 = 0. If, then, (ὦ, y) moves continuously 
along any curve in the ‘lower’ part of the plane, ὦ and v vary con- 
tinuously. But it is possible to move over to the alternative formulae 
a = 1—,/(z—2y+1), 
v =y—1+./(~7—2y+1) 
at any point where the path of the point (x, y) meets the line 
x—2y+1= 0. For either choice, the values of wu, v are continuous 
over the ‘lower’ part of the plane, even at the line. 
To summarize: suppose that, at any one point (x,,y,) in the 
‘lower’ part of the (x, y)-plane, the choice of solution 
i = 1+,/(%,—2y,+1), 
υ = y,—1—1(%,—2y, +1) 
_ has been made. The values of wu, v at any other point (x., y,) can be 
obtained by allowing τι, v to vary continuously with changes in 2, y 
such that the point (x, y) moves by a continuous curve (in the ‘lower’ 
part) from (2, y,) to (%2, ya). Ifthe points (5, ¥;), (2, ya) are joined by 
a curve not meeting the line x—2y+1 = 0, the values at (72, y2) are 
U= 1 (%_—2ygt1), ᾧ τὸ ψ5--ὶ -- ἰ(ας -- 2.52. 1). 
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If, however, the points are joined by a curve meeting the line, a 
change can be made if desired as the point (x, y), moving along the 
continuous curve, meets the line. Once the decision has been made, 
the formulae to be used for wv, v remain unique and must be kept. 


The only ambiguity in the choice of functions lies at the line 
x—2y+1=0, and it is worthy of explicit remark that u, v are 
genuinely two-valued functions of x, y there. If the point (2%, y;) 
lies on that line, then each of the pairs of values 

ἐμὲ 1+ νας -- 2... 1) a ΤΣ 1— V(%_—2y_+ 1) 

v = ¥g—1—/(%,—2y.+ 1) V = Yg—14+ γί, --υνς Ὁ 1) 
may be obtained at (x2, ψ4) while retaining continuity as the point 
(x,y) moves along a continuous curve from (2, y;) to (X2, Ye). 


We now turn our attention to the differentials of the given 


functions a=ur?+Qv, y=uty, 


namely, dx = 2udu+2dv, dy =du+dv. 
Solving these for du, dv by the usual process of elimination, we 
have the relations o(i 1) du = de = By, 
2(u—1)dv = —dx+ 2udy, 
leading to the partial differential coefficients 
Bo ϑι,. καριῦπον κ 

2(u—1)’ oy w—l’ 
FF tins egtielrs Wk de Stes 
Ox 2(ι-- 1) oy u—-l 

There is, however, one case where this analysis breaks down. 
When 


e= lI, 
the relations for dw, dv give simply the equations 
O= dzx—2dy, 


0 = —dx+ 2dy 
; , : Ou Ou Ov ov 
which do not involve du, dv at all, and the coefficients Bn’ By? Bx’ dy 
do not exist. Moreover, this relation u = 1 is exactly equivalent to 


the earlier relation ΜΕ ΣΟῚ 
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In other words, the differential coefficients ἐν ὅν, if ὅν do not exist 
at precisely those points (namely, those for which x—2y+1 = 0) 
where the expression of u, v in terms of x, y is not single-valued. 

It may be worth while to comment that the variables x, y are 
related at the ‘doubtful’ points for which wu = 1, by the expression 


z—2y+1=0; 
and that the differentials dz, dy are then connected by the relation 
dz —2dy = 0. 


Finally, let us generalize the functions to exhibit a clearer view 
of what is involved. Write 
c=f(u,v), y= g(u,»), 
where f(u,v), g(u,v) are single-valued functions when the point 
(ω, v) lies within a certain region of the (εν, v)-plane. The differential 
coefficients =, oe =, may, in general, be calculated by an 
appeal to the differentials of the given relations in the form 
«ὦ 5, Ὁ + of 


δ Ὁ, 


dy = δῦ σις Sav; 


for elimination of dv, du in turn gives the equations 


of eg of og og, of 
(nee 2 81) 4 wie te 


(22-4 2) ao - in ὦ 


Bu Dv dv du hed a 


and the partial differential coefficients follow, ἘΧΟῈΡΤ for values 
of ει, v such that 


The left-hand side of this equation is the determinant 
of of 


V2 9] 
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of the coefficients of du, dv in the equations for dx, dy, and is of 
great importance. It is the function whose vanishing in the par- 
ticular case x = u®+2v, y = u+v gave the relation wu = 1 to locate 


the points where ἐπ Ἐν Ἐπ By? τ did not exist. With generalizations 
to larger numbers of ἀν δ this determinant must now be 
studied more closely. 


EXAMPLE I 
Consider similarly the relations 


x=rcos?, y=rsind, 
and show that 
(i) the equations can be solved to give r, Θ uniquely (to within 
multiples of 27 for 0) in terms of x, y except when r = 0, the value of 
@ being then arbitrary; 


S Ox Oy Cudy _ 
(ἢ) or 8 OB er = ° 
ἘΞ Or Or 06 06 
(iii) the partial differential coefficients Ba’ By’? Bn’ ὃν do not 


exist when r = 0. 


2. The Jacobian defined. Taking, first, the case of three 
variables as an illustration, suppose that 


f(u,v,w), g(u,v,w), h(u,v,w) 
are three functions of the variables u, v, το. The determinant 


of 
bu 


Ἢ 

WML 
ML Me Me 

212 2k 


ch ah ah 

Ou Ow 
is called the Jacobian of the three functions, and is often denoted by 
the notation (fg, h) 


~ O(u, v, w)” 
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A similar definition provides us with the Jacobian 
Ai, Bags sj Be) 
O(%1, Ug, «++, Un) 
of the n functions 2,,2%,,...,%, of the πὶ variables w,, uo, ...,%,) in 
the form 
δα, θα) ... δα, 
δι Ug ... δὲ. 
J= Oly My» Ong A 
OU, Og ... Uy 
Ox, Ox Ox 


We shall, in fact, usually deal with the general number n, as the 
restriction to the number 3 affords no real simplification. 


InLustration 1. The variables x, y, z are expressed as functions of 
the variables r, 0, ᾧ by means of the relations 


x=rsinOcos¢?, y=rsin@sing, z=rcosé. 


᾿δία, Y, 2) 
To evaluate the Jacobian δίς, 0,8)" 
By definition 


O(x, y, z) 
O(r,0,9) ὁ 


VY Me VWF 
V2 WMele BWP 
sl? SL VIP 


βίη Θ οοβῴ rcos@cosdé —rsin@sing 
=| sindsing rcos@sing χσβὶπθοοβ ὁ |. 
cos 8 —rsind 0 
Hence, by elementary expansion of the determinant, 


(x, ¥, 2) 
a(r, 0, 4) 


= rsin 0. 
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EXAMPLES II 


δία, y) 

1. Iff τ γοοϑθ, y = rsin#, evaluate ——— a(r, δ)" 
2. Ifx = e“cos¢, y = e“sin ὁ, evaluate (29) 
δίω, 9) 
3. Ifa =rcoshd, y = rsinhJ@, evaluate ee 


4, Ifu=a2%+y?+22, v= ( ἘΣ Ἐ2)3, w= yz+zr+z2y, evaluate 
O(u, υ, w) 
O(a, y,z) 

5. If u=esind cos¢, v=esind sing, w=e'cos0, evaluate 
O(u, v, w) 


O(t, O, 9) ° 


3. The chain rule for Jacobians. ΤῸ prove that, if 11, 2%, ...;%y 
are functions of U1, Ug,...,U,, which are themselves functions of 


ξ,, Es, 9999 ξι then 
δία», Bigs «255 Bq) Ol, Wey .(ὅ., Up) _ O(%,, Xe, .. » Lp) 
O(Uy, Ug, ...» Un) O(Ey, Sa, «++» Sn) Ἢ OE, δ»... aE 
The proof rests on the formula for the product of two determin- 
ants, namely, that 


ἄμ ἄμ ... Ayn δι ὃ, ... On Cyy Cyq ... Cyn 


ie gle Ch δι δ «ὦ δ δὼ δὼ ὧς Com 
where Cy = 5a, 5,;. 
A=1 
[The reader may confine his attention to the case n = 3 if that is 
more familiar. ] 


If we write a. = ott hie Ou, 
Ou,’ 7 06,” 


80 that the two determinants to be multiplied are simply 


O(%y; Legs ...,35.) (hy, Migs 02+» hy) 
(Uy, Ue, em O(E1, $2 erry δ)" 
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Ox, Cu 
then = Pics, ales,’ 
C= ΟΣ Du, δὲ, 
Ox; 


; O(2,, 2 x,,) 
Hence the product determinant is ~~~", as we were to 
P (Es; σας ἧς my i i 


prove. 


ILLUSTRATION 2. ΤῸ verify the chain rule, in the form 


O(u,v) O(z,y)  O(u, v) 
O(a,y) δ(ν,0) O(r, 4)’ 


for the relations εἰ τε αν, v= 2zry, 
Ξε rcos6, εἶ = γ βίῃ θ, 
We have O(u, v) |=! 
O(v,y) | 2y aA 
= 4r? cos 20, 


= 4(x°—y’) 


and O(x, y) cosO πἰηθ 
a(r, 0) | —rsin@ rcosd | 
= f, 
O(u, ; 
Thus Fah ae = 4γ8 cos 20. 
Now u =r cos? 0 +r* sin? 0 
= 73, 
v = 2rcos@.rsin@ 
= r sin 20. 
Thus (uw, v) 2r 
δίγ, Θ) | 2rsin20 2r2cos 20 
= 4γ8ὃ cos 20. 
_ The formula is therefore verified. 


4. The ‘reciprocal’ theorem. 70 prove that 


O( 4, Uo, +++) Un) Ἢ APs, Bas «δ 
Ola e, Wes ashy Me) Oty, δὼ»... Uy) 
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In the preceding paragraph, identify ἔ;, £2, ---» ξη with 
4) Lg) +++) Lye Then 
δία, 22» «-+)Ly) _ Oty, We, «++» Mn) 
O(by, Way +++) Un) (Ay, Hay -++r Vn) 
_ δίαι, ἅ,».: «) Ln) 
"ais Bey »- +) Uy) 


1000 
0100 
bo: dent ,ἴῶ 
0001 
= 1, 


and the result follows. 


ΤΙ τάττον 3. The variables r, 0, ὁ are functions of x, y, z 
according to the relations 


a=rsindcos¢, y=rsindsingd, z=rcosé. 


Το μὰν σας 15.5.0) 
O(a, y, 2) Yy,2) 


The expression of r, 0, ¢ directly in terms of z, y, z is awkward; 
the subsequent differentiations are worse. We therefore use the 
ult (p. 96 
pre te (2,952) _ eng 
a(r, 8, d) “ 


a(r, 0, φ) 1 
O(x, y, 2) ~ sin 0" 


and find at once that 


Remark. The importance of Jacobians centres in many ways 
round two fundamental properties which we are now to investigate. 

First, they give a test to determine whether given functions are 
independent or connected by some ‘functional relation’—the 
phrase being used in the sense that, for example, the functions 

2 =sin?u+cos*v, y= cos*u+sin*® v 
are connected by the ‘functional relation’ 
z+y = 2. 
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Secondly, they give an indication of the ‘magnification’ between 
corresponding (small) figures in, say, the (ει, v)-plane and the (z, y)- 
plane or the (u,v, w)-space and the (x,y, z)-space. The analogy of 
the relation 
Ju = oho 
dx 

between lengths du, dx suggests itself at once. 

The second of these properties will assume importance later 
(p. 132) when we come to study multiple integrals. 


5. Dependent functions. ΤῸ prove that, if there is a functional 


relation connecting the n functions 21,%9,...,%» of the n variables 


Uy, Ug, +++,Uy, then the Jacobian 


ae Hs. Vins :0vsy Be) 


(Uy, Uy, «-+) Un) 
has zero value everywhere. ᾽ 


Suppose that the functions 2,,2%»,...,%, are connected by the 
relation f (2, %q, ...,%,) = 0. 
Differentiate with respect to w,, Ue, ...,u, respectively. Then 
Of Ox, Of 0%, Of oz, 
Be, Ou, * Bx, du, +" δα, du,” 
Of ὃν, | Of Oi, Of Ox, 


Gx, Bulg δα, Oty“ * Bex, Duly” 
Of Cx, δῇ Ox, Of Ox, 
be, ou, Oona 


Now it is a familiar theorem in elimination that, if there are n 
Renan, Qf, +dgf+...+a,£, = 0, 
6, £1 ὐ,ξ, Ὁ vee +0, δι = 0, 


€, 8, + Cobgt... + nb, =0 
among the n variables £,,&,...,£,, then 
Og Bog 000.) Re 
δε. ee sits 
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provided that £,,&5,...,g, are not themselves all zero. Applying 
this to the n equations obtained above, so as to eliminate 


Cae cs 
ΝΣ" ἢ τη μων 
we have the required relation (interchanging rows and columns in 


the definition) Ox, Oay ... ὅὃχ, 
Ou, Ou, ... Ory 
Ot, Oe ss OF 
Js| —4 = _"* |= 0, 
Cuz OUg 2 
OX, Oe 


provided that ef ᾿ ει sakes ἐς are not all zero; and the proviso is 
excluded by noting that, if it held, the given function f(x,, 2», ...,Z,) 
would be independent of each of the functions ,, » ...,%,—a case 


of little interest. 


We state without proof the converse result: 

If n functions of n variables have zero Jacobian everywhere, then 
there is a relation connecting them. 

InLusTRATION 4. 70 examine whether the three functions 


atyt+z, w+y*+2%, yzter+ay 
are related. 
The Jacobian of the functions is 


1 1 1 
Jm| 22. Qy 2 |, 
Yt2 τὰ τὰν 
and it is easy to prove that J = 0, so that the functions are related. 


(In fact, 
(v+y +z)? — ἰδ + y* + 2*) — 2(yz+2e+ zy) =0.) 


EXAMPLES III 
1. Verify that τα ᾽ = 0 for the functions 


x=sin?u+cos*v, y = cos*u+sin*v. 
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2. Determine whether the functions 


e+y+23, xyz, x+y+z2 
are related. 


3. Determine whether the functions 


αϑ y3+28—B3xyz, ewty+z, x+y?+2 
are related. 


4. Determine whether there are values of the constant a for 
which the functions 


e+yi+z22+axyz, στε, 2*+y?+22*-yz—zr—xzy 
are related. 


6. Ratio of areas under transformation. (Approximate 
theory.) The work now to be given is of fundamental importance, 
and exhibits one of the most characteristic features of Jacobian 
theory. Detailed discussion is difficult at the present stage, but 
the principles should be followed closely. 

Suppose that the relations 


ἴω = f(x,y), shan σία, ν) 
define u, v as single-valued functions* of the variables x, y. Take 
(x,y) as coordinates in one plane and (u,v) as coordinates in 
another. 

To each point P in the (2, y)-plane corresponds a point P’ in the 
(u,v)-plane, and to any figure in the (z,y)-plane corresponds a 
figure in the (u,v)-plane. Our purpose is to prove that, if A is the 
area of a small figure round the point (x,,y,), and if A’ is the area of 
the corresponding small figure round (u,,v,), then, approximately, 


Ἄς O(U, O(Uy, V3) 
the sign being adjusted to make the right-hand side positive. 
To do this we break the area A up into a number of small triangles; 
the figure A’ is then broken correspondingly into a number of small 
figures which are also triangles (with sides approximately straight). 


* It is understood that the relations can be solved to give 2, y as single-valued 
functions of u, v, possibly for restricted ranges of values of the variables; for 
example, the relations u=2z*, v=y* can be solved uniquely in the form z=./u, 
y=./v for points in the positive quadrants. See also the Note at the end of the 
next paragraph (p. 106). 
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Once the theorem is established for a small triangle, the general 
theorem will follow by addition (since all signs are adjusted to be 
positive), so we confine our attention to triangles. 
Let us take, then, a small triangle with vertices P, (x,,¥;), 
Ps (2; Y2), Py (%s, Ys), where 
=%th, y=yth, 
ὥς ε Φιρκι ἢ, yg=Yytk'’. 
It is known that the area of the triangle P,P, P, is A, where 


% yy 1 
2A=+/% y 1), 
ὥς Ys 1 
the sign being adjusted to make A positive. Thus 
vy 1 1 


24 --΄Ψἰ] αιτῇ y,+k 1 
ath’ ytk 1 
A & 
=+ ᾽ 
kh’ 6’ 
after reduction by subtracting the first row from each of the others. 
In the (wu, v)-plane we have a triangle P; P; P3, the coordinates of 
whose vertices are given by the relations 
Uy = f(x, Y1); y= g(x, Y); 
Us =f(%2,Y2) =f(tith,yit+h), Veg= σίας, Ye) = σία, +h, y, +h), 
Us = f(%3,Y3) Ξε [5 +h’, y, +k’), 3 = ρίας, Ys) = G(X, +h’, y, +k’). 


Expanding by Taylor’s theorem, we obtain the approximate 
relations 


δ 8 ὃ ὃ 
+h——+k—— Ve = 1,t+h—+k—, 
“δέ. hla 7% a tay OY; 
of 2, af ,o9 4, ag 
=4,+h’——+kh’—, v=9,+h'’—+hk' —. 
κύνα Say κ᾽ pA ποτ Ry 


Hence, by a reduction similar to that adopted for A, the area of the 
triangle P’Q’R’ is A’, where, approximately, 
f+. og . ὃ 
h—=—+k—— h—+k— 
dm," Oy,’ dn, OY, 
Ε] 
a me: BAK 


Oy,’ Oat Oy, 
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and so, using the well-known theorem on the product of two deter- 
minants (or verifying by direct calculation), we have the formula 


f ὃ 

h Ox 

2A’ = + . x ag! ἣ 
h’ κ' of og 

Oy, OY; 


if o(f,9) 
538 O(a, γι) 


Hence A’'=+ 4.36) 
O(24, ψι) 


Remembering that 
%=f(%1,%), υχ Ξε σία, ψ.), 
we obtain the formula for a small area in the vicinity of a point 
(%,¥,), in the form pe phe ates,, δὴ 
δ(α,, 41)’ : 
the sign being adjusted to make the right-hand side positive. 


7. Ratio of volumes under transformation. (Approximate 
theory.) The formula of the preceding paragraph may be extended 
readily to three variables. Suppose that the relations 


ἴω = f(x,y, 2), Ψψ σία, ψΨ, 2), ι-Ξ h(x, y, 2) 
define u, v, w as single-valued functions of the variables z, y, z. Take 
(x,y,z) as coordinates in one three-dimensional space and (u, v, ὦ) 
as coordinates in another. We prove that, if V is the volume of a 
small figure round the point (x1, y,,2,), and if V’ is the volume of the 
corresponding small figure round (u,, v1, W,), then, approximately, 
O(u,, V1; 10.) 
V’ = .---- -- ΤΥ, 
* δία, ψι, δι) 
the sign being adjusted to make the right-hand side positive. 
Breaking the volume into small tetrahedra analogous to the 
triangles of § 6, we consider a particular one, of vertices (2, ¥;, 24), 
(2, Yo, 24)» (Xs, Ya. 253)» (Xa, Ya» 24), Where 
ἄς παι, Yeo=AWwtd Ζ, ΞΖ. ΤΥ, 
ἄς τ +p’, Yy=YWitT, ὅς τ 2. ἘΤ΄, 
᾿ς Ξε Δ .ἘΡ", ψετπε νι 4", Ζ, τ 2. "75. 
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The volume J is given by the formula 


%y YW 


6V =+ 


εξὶ 

δὰ 
δ εδ (δ 
μι μ͵λ μα ὑπὸ 


5} σον 
»" φ' ν" 
on subtracting the first row from each of the others in turn. 


Aftertransformation we obtainatetrahedronof vertices (ω,,υ, ει), 


Where y, = f(x, ψι» 2), = θα,» ψι» 21), Wy = Wey, ψι»2ι), 
Ug τὸ [(ς Ἐ}, γι Ἑ 4,21. 7), Y= 9(%,4+7,y¥,+9,%, +7), 
τῦς = h(x, +p,y¥,+9,% +7); 
etc., and so, using Taylor’s theorem, we have the approximate 
relations of of Ὁ 
“ΠΈΣΕ +a τῆς, 
Ug = Uy Pay q ὃν, Oz, 


oh h 
W, = mips tas tre, etc, 


Hence, approximately, the value of 6V’ is 


OP τς δ ey. og oh oh oh 
wate ae oe, re) Se: tome, Bch per 0 ae 
᾿ OX, ᾿ OY; Oz, ᾿ rq ay, ᾿ +4 oy,” Cay 
» Of oF . ie. et ee. ae ae δα 


re, Pata tr, ρ»"- +9" +9" —— 
δε." Τα, Τὸ Oe’. ” Be, ἢ by,” ὅς, 


o(f, σ, Ὁ) 
“ΟΜ τος 4. 
Ξ “75 Yx,%)’ 


8 ΜΙ 
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͵ οἱ » Vy, ι) 0 0 C 
that Vi = + Von A) Τὴ ὦ if if 
so tha = (Ya %) State without proof the formula for () im terms of ( ) ; ΕἸ ). : 


ov 


Note. The formulae ν (=) , and show that 


(a), 


Ζ 
Aa Ale) γος pO) “ af. 9) _ af, σ) 2a) 
O(a, Ya) O(a» Yas %) δία, ) (u,v) δία, y)" 
imply that the Jacobians Hence, or otherwise, show that 


O(Uy,%;)  O(Uy, Vy, Wy) 

O(%,Y1)’ (4, ψ., 2) 
do not vanish. We have seen in §1 (p. 90), that this vanishing is 
associated with points where the ‘inverse’ of the given trans- 
formation is not determined uniquely; but we do not propose to 
go more fully into the phenomenon. 


5). (e),~ ae 
" ῳ.--(ἢ, a), 


4, Show that, if τὸ, v are functions οὗ x, y such that x, y are 
expressible uniquely as functions of uw, v, then 
Ou dx Ov dx Ou dy dv dy 


aon deo “ GeOu δὲδν 
REVISION EXAMPLES XII with other similar equations. 
University Level If Gu Ou 
δ(λ,μ). - (u,v) _| ὃς ὃν 
1. Evaluate δέ, 9) in terms of A, μι, where a(x, ψ) ὃυ ov , 
2 y? at y? Ox dy 
Ἅπας πα ἐτ--Ὁ ὅν l, — +——_ = 1. 
at+A’b+A atm b+pu prove that et ba 
2. Ifa, b,c, ἃ are constants such that ad—bc+ 0, and if | wee ; si Ὁ; 
z=e+y®, w=ab+by, y=ch+dy, a7 es) 5 (735) 
show that 2 Oe 2 οὗ 5. When v is eliminated between the equations 
Δ ϑυῦμι Minit δῦ [εξ a y=f(x,v), 2=9(x,0), 
Oz 8% Oz δὲς [ἰδίω,Ψ) the equation z= f(x,y) 
Gyo dy? δηδξ dy? is obtained. Prove that 
3. Six related variables f, g, u, v, x, y are such that any two of a = ἘΞ τ 
u, v, x, y can be taken as independent variables and the others are : 
differentiable functions of these two. The partial derivative of f Verify this result when . 
with respect to x when y is constant is denoted by (ef/¢x),, and so y = xcosu—asinv, 
on, and 2 2 2=xsinv+acosy, 
(2) (2) a being a constant. 
O(f,g)_| θυ), \ey/s : , ᾿ ς 
O(a, y) ag ag 6. Each of the variables 2, y, z is a function of the variables 
=) (2). u, v. A partial derivative such as τς is formed on the assumption 
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that v is kept constant. The symbol (=) implies that x is expressed 


ν 
as a function of εὐ, y and y is kept constant during differentiation. CHAPTER XVII 


The symbol 2) is defined by the relation MULTIPLE INTEGRALS 


O(u, v) 
ox cy The ideas already explained for the integrals of functions of a single 
O(a, y) _ du ow variable may now be extended to ‘multiple integrals’ of functions of 
O(u,v) | δὲ ey ; several. The aim of this chapter is to give a clear picture of what is 
ὃν ov involved, but it must be recognized that the treatment is introduce 
δι alae, y) [ey tory. Once the general processes are understood, a detailed investi- 
Prove the relations (=) = cera [Ξ. gation may be found in a text-book of analysis; until the processes 


are understood, the necessity for the details can hardly be grasped. 
We begin by the extension to functions of Two variables. 


(~) _ @y,2) / (ez) 
Ox), O(u,v)} (u,v) 


7. If, y are defined as functions of z by the equations 
f(x, y,z)=0, g(x,y,2z) = 0, 


1. Double integrais. Let 


f(x,y) 
be a function of the two variables x, y, defined for values of x, y in 


Sasi ΒΗ dx o(f,g) [&(f.9) the region R of the (x, y)-plane (fig. 114) bounded by a simple closed 
dz δίψ,:)} A(x,y)’ curve C (without ‘ crossings’). 
where of ag y 
O(f.g)_| ὃν ὃν 
δίψ,) | af ag | 
dz 02 
Find in this way when the defining equations are 


22+ 2z—x—-y+2=0, 
xsyt+y?—x+y—-1=0. 


8. If x, y, z are functions of uw, v with continuous first-order 
partial derivatives, show that 


2) ay 4 O22) gy 4 9) ἀν -- 0, 


(uw, v) ma thu, v) *3(u, 0) v) 
where cy oz Fig. 114. 
δίψ, z) _ Ou du 
(u,v) ὃν |" Divide the region R, in any manner, into a large number of small 
a» Ἂς elements, of which a typical member has area w,;, and let M;;, mj; 
Oe be the greatest and least values of f(x, y) in w,;. Form the two sums 
Find δ, in terms of τ, v if 84 = UM; ὡς» 


a=uit+v, y=utv’, z= err, 84 = LM; W;4;, 
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summed for all the elements w,;. We call S;;, 8), the wpper and lower 
sums respectively for this particular subdivision. 

Now suppose, analogously to integration for a single variable, 
that the number of elements ὡς; is increased indefinitely, the size 
of each being diminished indefinitely. In favourable cases, S;; tends 
to a limiting value S which is independent of the method of sub- 
division, and s;, tends similarly to a limiting value 8. 

The function f(x, y) is said to be integrable over the region R if 


S=s. 


The common limit is then called the pouBLE INTEGRAL of f(z, y) 
over R. 

In what follows we shall assume the integrability without further 
discussion, restricting our choice of functions to that end. 

The actual evaluation, which is our main purpose, is made simpler 
by noticing that, if (€;,7;) is any point in the element of area w,;, 
then, by definition of ᾿ς,» πῃ 


Mi; >f (Ei 95) > Mays 
so that UM ;,0;; > Xf (F;, 0;) Weg > VIN Wy. 


On proceeding to the limit, the two outside sums assume the same 
value, to which 
Σί(ξ,, 3) Ox 


must therefore tend also. Hence the greatest and least values M,;, mj; 
may be replaced in the definition by the value of f(x,y) at any point 
in the element of area w,;. 


InLusTRATION 1. Z'o find the moment of inertia of a uniform 
rectangular lamina, of sides 2a, 2b and density p, about a line through 
the centre parallel to the sides of length 2b. 

Take axes, as shown in the diagram (fig. 115), through the 
centre and parallel to the sides. Divide the rectangle into a large 
number of small elements of which a typical member has area w,,, 
and let P(x;,y;) be a point inside w,;. The moment of inertia about 
the axis Oy is defined to be the limit (if it exists) of the sum 


PLT Wy 


taken over all the elements as their number increases indefinitely 
while the size of each decreases indefinitely. 

We have said nothing so far about the method of subdivision, 
but obviously an orderly system must be used if the calculation is 
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to be at all manageable. The obvious way here is to divide by lines 
parallel to the axes (fig. 116). We therefore begin the calculation by 
considering the elements in a certain strip parallel to Oy (denoted 
by darker lines in the diagram), later completing the summation 
by adding the results obtained in all such strips. 


¥ 


Fig. 115, 


Fig. 116. 


Suppose, then, that the strip parallel to Oy has a fixed width 
which we may call dx,, and that a typical element in it has height 
dy;; then ὡς; = δα, δύ. Suppose, too, that a typical point inside this 
element has x-coordinate ας, taken to be the same for each element 
of the strip. The required sum is 


ρ Σ Σ aj OX; OY;. 


To sum the elements in the strip, we make the summation with 
respect to 7, keeping x, and dx; constant. By the definition of simple 
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integration with respect to the variable y, this summation leads in 
the limit to the expression 

b 


ed ajox,| dy 
i -ὸ 


ὃ 
Ξ-ρΣ αἵ oe υ] 
ἐ —b 
= 2bp > xz Ox;. 
i 


To sum over all the strips, we now let δῷ; become indefinitely 
small and proceed to the limit. This gives 


οὐρὶ atda 
—a 


=f 


= $a°bp. 
If M is the mass of the lamina, so that M = 4pab, then the moment 
of inertia about Oy is Ma? 


ILLUSTRATION 2. ΤῸ find the moment of inertia of a uniform lamina 

of density p, bounded by the ellipse 
oy? 

att pa) 

about the axis Oy. 

The moment of inertia is, by definition, 

PLZ 4y, 
as before (p. 110). 

For the method of subdivision, we again proceed by drawing 
lines parallel to the axes, and then summing ‘up’ a typical strip 
between the lines 

t= Lis “= XL, + δα, 
as indicated in the diagram (fig. 117). The line x = x; meets the 
ellipse in two points whose coordinates are 


bo JCB) fo 8/2) 


The strip is divided into a number of rectangular elements by 
the lines y = const., save that there are ‘end-effects’ at the curved 
boundary. It is easy to arrange the subdivision in such a way that 
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there is only one irregular element at each end of the strip—for 
example, the arc of the ellipse may be divided into a number of 
equal sections and the lines parallel to the axes drawn through the 
points of division. 

The sum ρΣαζω,, 


may now be considered under two headings: 

(i) the sum for all those elements w,; which are complete 
rectangles; 

(ii) the sum for the remaining elements, at the boundary of 
the ellipse. 


Fig. 117. 


We first prove (what the beginner may prefer to regard as obvious) 
that the sum (ii) tends to zero as the subdivision approaches its 
limit. 

The value of 2? is certainly less than a*, and the area w,, of an 
element is certainly less than its ‘height’ dy; times the length ds, 
of the element of arc bounding it (fig. 118). Thus the sum is less 


= Σραϑ ds, dy 5 


Hence, if y (which will tend to zero) is the greatest ‘height’ of any 
rectangle in the whole subdivision, so that dy;<7, the sum is 


less than ραδηΣδα,, 
or pa*nP, 
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where P is the perimeter of the ellipse. In the limit, as 7 tends to 
zero, this sum also tends to zero, and so the sum (ii) is zero in the 
limit. 

Since the ‘end-effect’ is zero, we may neglect it at once, and 
regard the strip of elements as running from precisely the value 


-» /(1-2) to τὸ /(1-%). Summing up the strip, we then 
have the sum +by(1—2¢'/a") 
praz dx, Ϊ y 
—by(1—z*/a") 


= 2pEba3 | (ι- 2) δα; 


ὃς; 


Fig. 118. 


Summing now for all the strips, we obtain the limit 


20)" a {{1-- ὃς de, 


If x = asin @, this is 
ἐπ 
2pba* sin? @ cos @ cos 0d0 
-π 


ἐπ 
= 4pa%s|sin?2040 
—j}7 


= teat" (1-008 46) 40 


= ζπραϑῦ. 
Since the mass of the lamina is 
M = πραῦ, 
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the value of the moment of inertia is 
1Μα3. 

2. Notation. The method just illustrated for evaluating 
double integrals carries within itself the suggestion for the notation 
to be adopted. For Cartesian coordinates, where the subdivision 
is by rectangles, the integral is the limit of the expression 

Uf (wis Ys) 0x, SY; 


(evaluating the function, for convenience, at the corner (x;,y,) of 
the rectangle), and is written in the form 


| f(x,y) dudy. 
We have just proved that the value of 


| | δρῶν 


over the rectangle | «| <a, | y| «ὦ is 


30°, 

and that its value over the interior of the ellipse b®x? + a®y? = a*b? is 
ina*b, 

The value of _ | dady 


over the region R bounded by a simple closed curve is equal to the 
area of ἢ. 

The preceding illustrations show how, in practice, a double 
integral is usually evaluated by considering it as a repeated integral 
with respect to the two variables in succession (in either order). 


3. Triple integrals. The work in space follows by natural 
extension from that in a plane. Let 


f(x,y, 2) 
be a function of the three variables x, y, z, defined for values of 
x, y, z in the region R enclosed by a simple closed surface F. 
Divide the region R, in any manner, into a large number of small 
elements, of which a typical member has volume τῳ» and let V;;;,, 
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m;;, be the greatest and least values of f(x, y,z) in τῳ. Form the 
wpper and lower sums Sop DM Tue: 


8... = LM gE isn: 

summed for all the elements το. If S;;, tends to a limiting value § 
and 8;,,, to a limiting value s (independently of the method of sub- 
division) as the number of elements το increases indefinitely, the 
size of each decreasing indefinitely, and if, further, S = s, then we 
say that the function f(x, y,z) is integrable over the region R. The 
common limit is then called the TRIPLE INTEGRAL of f(x, y,z) over R. 

As in the case of double integrals, we may for actual evaluation 
take the value of the function at a point (ξ;, ἡ,» ¢) inside το» and 
proceed to calculate the limiting value of the sum 


uf (ξ;» 7,» Cx) ΤΙ" 

When the coordinates are Cartesian, the method of subdivision 
is naturally into ‘boxes’ of volume δα; dy; dz,, and the limiting value 
of the sux Uf (Xs Y5, Ξε) δα OY; δῶ, 
is denoted by the symbol 


Ϊ f(x,y, 2) dxdydz. 


The method of evaluation usually consists in considering the 
triple integral as a repeated integral with respect to the variables in 
succession, taken in whatever order appears likely to be most 
convenient. 


InLustrRaTIon 3. 70 find the volume enclosed by the ellipsoid 

whose equation is ¥ vets : 
that 

Divide the ellipsoid nhs small th ’ by planes parallel to the 

axes of coordinates. A typical box, situated at the point (x,, y;, z;,) 

has volume δα, dy;6z, (with small corrections, which are negligible 


= 1, 


for sufficiently fine subdivision, at the surface itself*) and the total — 


volume is the limit of the sum 
Σδα, dy; d2,,. 


To begin with, keep x,, y,; fixed, and integrate with respect to z 
up the ‘tube’ of cross-section δα, dy;. This gives, for the volume of 
the tube, the formula 

Ldx, dy;| dz, 


* Compare the more detailed treatment on p. 113 for an analogous problem. 
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between appropriate limits, or 
Σ δα, δὲ, [2]. 


Now, when x = 2;, y = y;, the value of z varies along the tube from 
(approximately) 


to το /( -ξ.-"Ὁ, 


The totality of these tubes may be located by their sections in 
the (2, y)-plane z = 0, which they cut in small rectangles, covering 
the ellipse εν 

+3 
in that plane. In order to add all the tubes, we first keep, say, 
a, constant and sum for all values of y within the ellipse, giving 


2czer,[ (1-2-5) ay 


between appropriate limits of integration, namely, 


yall 1-4/3) 


To effect this integration, write 


{{8-- 


=] 


so that the integral is 


the substitution y = asin? 
ἐπ 
then gives i acos@.acos6d0 
- ἐπ 
.» 5’ 5. 
ὃ 2° 
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Hence the volume is 2cXdx,—. x δ3{1-- “ἢ, 
26 α 
or, in the limit, 


. x 2a? 
mbe| (1 -=) dx = mbe (20-5 a 


= $nabe. 


EXAMPLES I 


[It is most important that the reader should, at this stage, 
develop the faculty of, as it were, following the integration mentally 
up the filaments and then seeing the filaments move to cover the 
whole figure. The Illustrations in § 4 show how the calculations are 
effected in practice, discarding the language of summation and 
proceeding straight to integration; the examples now given should 
be written out in detail like the above models until the picture of 
the process is completely clear. The later work should then be easy 
to follow.] 


1. Find the moment of inertia, about the y-axis, of a uniform 
lamina of density p bounded by the straight lines 

y=l+a7, y=l-2z, y=0. 

2. Find the mass of a lamina, lying in the positive quadrant, 
bounded by the axes and the circle 2?+y?= 1, given that the 
density at any point is k(1+2*). [Consider Xk(1 +27) dx, dy;.] 

3. Follow the method given in Illustration 3 (p. 116) to prove 
that the volume of a sphere of radius a is ὅπαϑ. 


4, Find the mass ofa sphere of unit radius, given that the density 
at distance r from the centre is &(1 +,r?). 


[Consider ΣΚ( + a3 + y} + 20) da, dy, dz,.] 


4. The evaluation in practice of multiple integrals. We 
now give some typical examples to show how multiple integrals 
are evaluated in normal practice. The picture is that of summation 
- of elements; the language is that of definite integration, 


I={{ aydedy 


over the area in the first quadrant bounded by the straight line x = 1, 
the circle x? + y? = 8, and the parabola y = 32". 


(i) To evaluate 
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For a given value of x, the value of y ‘up’ the line # = const. 
varies from 42" to -Ὁ γ(8 -- αΞ) (fig. 119). Hence 


I - [rae] 


- Al xdx(8—x* -- 1χ ἢ) 


=5] (Se —a?— 425) de. 


Fig. 119. 


Now the range of values for x is (as the diagram indicates) from 
x = 1 to the value of x when the circle meets the parabola, namely, 
x = 3, Hence 


ζω 5 {_(@2—29—429)de 


=a] aot γος got | = 1103-90 


= 213, 


(ii) Τὸ evaluate I={{ xevdedy 


over the rectangle bounded by the lines x= 1, x= 2, y= 1, y=k, 
where k>1 (fig. 120). 
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shows at once that it may be used to evaluate the volume, ‘above’ 

the region ἢ, contained between the surface z = f(2,y) and the 
ὺ 

plane z = 0. (Compare the integral | f(x) ἀξ as the area, ‘above’ 
a 

the interval (a,6), contained between the curve y = f(x) and the 

line y = 0.) In fact, Μ᾿ ως; is the volume of a cylinder, standing on 

the area w,;, of height equal to or greater than f(z,, Y;); and m,;0,; 

is the volume of a cylinder, standing on the area w,,, of height equal 


to or less than f(x,,y;). The whole volume, obtained by summing 
these cylinders and proceeding to the limit, is thus 


[| se.) azdy, 


ItLustratTion 4. ΤῸ find the volume of the ‘box’ whose base is 
the plane z = 0, whose sides are the planes x = —1, x = 1, y = —2, 
y = 2 and whose ‘top’ is the surface given by the equation 


x? + 5y?+2z = 80. 
The volume is | (80 —a* — By?) dady 
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Keeping ᾧ constant, first integrate with respect to y, from y = 1 
toy = k. Since 


' 
| xe*¥dy (x constant) 
1 


= le] = (εἴ —e*), 
ν 


over the rectangle x = 1, y = +2. 
Integrating first with respect to x, we have 


| | 80x δον ὯΝ ἡ ἂν 


5 --1 


- [Γ΄ (00--α- 5y4)—(—80+ 4+ 5y%)}dy 


Fig. 120. 


we have t= { (e—e)dz 


ee 3 
=|; -ὡ 
1 


1 
- Ὁ. ὅν ed 
τ΄ τὰ e* +e 


5. The double integral as a volume. We have seen (p. 13) 
how the function f(x,y) is ‘represented’ by the surface whose 


equation is z= f(x,y), 


where z is the ‘height’ of the surface ‘above’ the point (x, y) in the 
plane z = 0. The definition of the double integral 


| f(@, y) dudy 


by means ofthesums 2M,;0,;, Lm; 


= [ (1504-1042) dy 
[sw 
= (235 — 80) — (258 + 80) 


= 584. 


6. ‘Elements of area.’ Double integrals. In the definition 
of a double integral as a limit of summation, we dealt with the 
argument in terms of Cartesian coordinates, so that the element 
of area w,; was naturally chosen as a rectangle dx, dy;. The next 
problem is to consider what form of element should be taken when 
polar coordinates are used instead. 


9 M ΠῚ 
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Suppose that the function to be integrated, over a region R 
bounded by a simple closed curve C, is expressed in terms of polar 
coordinates in the form f(r,0) 


Dividing the region into elements of area ὡς; and taking the value 
of the function at a point τ,» 0, inside w,;, we reach as before the 


summation Ef (rg, θη) Wy. 


Fig. 121. 


The value of the integral is the limit of this sum as the number of 
elements is increased indefinitely, the size of each decreasing 
indefinitely; it is assumed, as usual, that this limit isindependent of 
the method of subdivision and of proceeding to the limit. What 
we have to do now is to find the most convenient shape for the 
element w,; and then to obtain a corresponding expression for 
its area, 

The curves r constant are concentric circles whose centres are 
at the pole O, and the curves @ constant are straight lines through O. 
We are therefore led to consider an element, such as that shaded 
in the diagram (fig. 121), bounded by ares of circles of radii r,, r;+ dr; 
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and by segments of straight lines inclined at angles 0;, 0,+60, to 
the initial line Ox. Subdivision into elements of this type will 
obviously be convenient for work in polar coordinates, provided 
that the expression for the area assumes a reasonably simple form. 

Now the area of a sector (fig. 122) of angle 60; and radius 7, + dr; 


pa pcan fino Hrg+ dr.) 80,. 


(It bears the ratio 60;: 27 to the area of a complete circle of radius 
γι δὴ.) Hence the area of the element, being the difference of two 
sectors, of radii r;+ δὴ,» r, respectively, is 
A(r; + dr,)? - r?} 30; =f; or; 60; + 3(dr,)? 06;. 
The summation required for the integration is thus 
Σίν, 9;) 7,67, 80; 
ἜΑΣΙ, 95) (0r;)? δθ,, 

where we have temporarily used the ‘corner’ (r;,0;) of the element 
ὡς, a8 the point at which the function is 


evaluated. x bth 
Consider the second summation 
U=$2rf(r;, 0;) (dr,)? 00;. Fig. 122. 
We prove* that it tends to zero with dr,, 60;. 


It is assumed that the given function is bounded throughout the 
region F (fig. 121); suppose that it is numerically less than a certain 


constant K. Then, U <4K3X(8r,)200,, 
numerically. 

Suppose next that, when the subdivision is made, all the elements 
dr; are less than a certain number p, which itself is to shrink to zero 
for the limit. Then, numerically, 


U <4 KpxXor,00,, 

where we have succeeded in making the summation linear in dr;. 

Suppose further that the region # is bounded in extent, so that 
all points of it lie within a certain distance D of the pole O (fig, 123). 
For summation along a fixed sector δϑ,, the value of Xér,; certainly 
cannot exceed 2D, the diameter of a circle containing the whole 
region as Hence U < KpD=00,, 
numerically. 

* The proof of this point may be omitted at a first reading. 
9-2 
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Finally, the whole region R is covered as the radius vector 0; 
swings from 0 = 0 to 0 = 27 (or less), so that, by summation of 06;, 


U <2n7KpD, 
numerically. 
. Since K, D are bounded, the right-hand side tends to zero with 
p, and so U tends to zero also. 
Returning to the summation 
of (r;, 0;) "; or; 60; 


+42f(7;, θὴ (dr;)? 06;, 


----- ,-- 


| 
| 
| 
ἢ 
| 
| 
| 
| 
i 


Fig. 123. 


we are now able to ignore the second expression, since it vanishes 
in the limit. Hence the value of the integral is the limit of the sum 


Sf(re, θη) τεδειδθ, 
namely, | f(r, 0)rdrdé. 


This is therefore the formula for the integration of the function 
(γ, 0) over the region R, 


In actual practice, the formula just obtained is usually reached 
by directing the approximation to the element of area itself rather 
than to the second summation which we have just considered. Thus 
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the elements are conceived as approximate rectangles of which two 

adjacent sides are of length 7,60;, dr;, giving an approximate area 
7,00, δ, 

in agreement with our formula, As a mnemonic this treatment is 

useful, but it seems more satisfactory to trace the fate of the total 

‘correction term’ and to see that it does really vanish. 


At any rate, the formula for the double integral of a function 
f(r,@) over a region R assumes a value which is the limit of the 


summation Ef (r θη) "87, 80,, 
written in the form | Ϊ ιν. δ)νάναδ, 
ILLUSTRATION 5. 70 find the moment of inertia of a uniform lamina 
in the form of a cardioid r = a(1+cos8@), 
aboutaline through the origin perpendicular 


to its plane. 
The moment of inertia is, by definition, 


ρΣη ἴδω,» 
or, in the limit, 
p| | ()rarao, 


taken over the interior of the cardioid Fig. 124. 
(fig. 124), where p is the density. 

Keeping @ constant, integrate first with respect to r ΘΙ r=0 
to r = a(1+cos@). The result is 


p Ϊ ἀ017γ3] = fatp| (1+ ο08 60) dé. 


Now the whole cardioid is traced as 0 moves from zero, through 
7, to 277. Thus we have the expression 


ie'p) "(1 44.0080 + 6.00820 +4.008* 6 + cost 6) 49, 
Qn Qn 
eae | cos θά = | cos*d0 = 0, 
0 0 


27 27 
| 606° Od0 = a, | cos*Od0 = $n. 
0 0 
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Hence the moment of inertia is 

I = jatp(27 + 67 + $77) = ξὅ παρ. 
Also the mass is, similarly, 


Μ = pf | rdrdo = pf dota 


= jatp| (1 +0080)*d0 = Yo%p(2n +n) 


= ὅπαϑρ. 


Hence I = 35 Ma’, 


Note. It may be helpful to use this Illustration to clear up a 
point which sometimes confuses the beginner. The sum 


pur’ Or; 60; 


is really being interpreted in two ways: 
(i) we derive it first in the form 
pxur}.r,dr,00; 
as the integral of the function pr? for elements of area r;dr,60; over 
the interior of the cardioid r = a(1+ cos @); 

(ii) we treat it next as the sum of terms pr? for elements of area 
6r,60;, as if the variables r;, 0; were RECTANGULAR CARTESIAN 
COORDINATES. (It was in this way, though not explicitly stated, that 
we reconciled the sum pir3er,00, 


with the formula purty 
leading to the integral p Ϊ γ3ἀν4θ.) 


The region of integration with 0, r as Cartesian coordinates is that 
between 0 = 0 and @ = 27 lying ‘under’ the curve r = a(1+cos@) 
—just like the area between z= 0 and x = 27 lying under the 
curve y =a(1+cosz). This region is indicated in the diagram 


(fig. 125). 
In other words, we may say that the algebra involved in the 
summation p&r3dr,80, 


for 0<r<a(1+cos@), 0<0< 27 has two ‘geometrical’ interpreta- 
tions: 
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(i) in polar coordinates r, 0, the integral of pr® over the interior 
of the cardioid r = a(1+cos@), in accordance with our original 
derivation; 

(ii) in Cartesian coordinates 0, r (compare 2, y), the integral of 
pr* (compare py*) over the area bounded by the straight lines 
r = 0, 8 = 0, 0 = 27 and the ‘cosine’ curve r = α(1] - 6080). (Com- 
pare the area bounded by y= 0, x= 0, 2 = 27 and the curve 
y = α(Ἱ] - 6008 2).) 

᾿ 


Fig. 125. 


7. ‘Elements of volume.’ Triple integrals. 

(i) CYLINDRICAL cooRDINaTES (p. 2). When a function to be 
integrated throughout a given volume is expressed in terms of 
cylindrical coordinates in the form 


f (ρ, ῴ, 2), 
the summation Uf (Pis Djs Ze) τας 


is easily effected by taking the elements of volume as ‘slices’, 
z = constant, standing on elements of area like those used (p. 122) 
in a plane for polar coordinates. By immediate extension we 


obtain the formula 
ΣΙ (Pi, Pj; %) Pi Op; δῴ, ὃζ,, 


leading to the limit Ϊ Ϊ Ϊ Νριᾧ, Ἰράραφαν. 


As an example, we obtain a formula which we shall use almost 
immediately afterwards. (Several other methods of calculation 
are available.) 
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InuustTrRaTiIon 6. ΤῸ prove that the volume of a ‘sector’ of a sphere 
of radius a bounded by a cone of vertical angle a, with its vertex at the 
centre of the sphere, is 27a2(1— cosa). 


Take the axis of the cone as the z-axis and the centre O of the 
sphere as origin for a system of cylindrical coordinates p, ᾧ, z. The 
volume is to be calculated from the formula 


Fig. 126. 


Keeping p, ¢ fixed integrate with respect to z up a ‘filament’ 
such as QP in the diagram (fig. 126). If this line meets the plane 
z = Oin NV, the length ON is p. Then, since P is on the sphere, 

NP = (a?—p?)t, 
and, from the triangle ONQ, in which ZOQN =a, 


NQ = ροοΐα. 
Hence the volume, being 


[fowl], 


is {| pdpd¢ {(a* — p*)t — p cot c}. 
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The double integral is to be taken over the circle of radius asin « 
in the plane z = 0 (shown dotted in the diagram), which is the pro- 
jection on that plane of the circle in which the cone meets the sphere. 
Keeping p fixed, integrate with respect to ¢ from 0 to 27, giving 


asina 
V = 27 Ϊ p{(a*—p*)t—pcota}dp. 
0 
By elementary integration, we have 
V = 2n[—}(a*—p2)} — 4p cot Jenin 
= 2na*{(—4}cos*a—4sin* « cot ~)+ 4} 
= ὅπαϑ(] -- οο88 ~—sin® & cos a) 
= ξπαϑ(] -- ὁ08 α). 
(ii) SPHERICAL POLAR COORDINATES (p. 2). The evaluation of 
th 
ie UP (ris 5s Pre) Tage 


for spherical polar coordinates is more awk- 
ward. The surfaces r constant are spheres, the 
surfaces @ constant are cones, the surfaces ¢ 
constant are planes, and it is natural to use 
them to define the elements 7,;,. We shall 
therefore proceed as follows.* 

(a) Calculate the volume of the part of a 
sphere of radius r between two cones of 
vertical angles 0, 8+ 60. | 

(b) Calculate the difference of the volumes 
obtained in this way from two spheres of radii 
r, r+ Or; the figure may be visualized as the 
‘ring’ obtained by rotating the area shaded in the diagram (fig. 127) 
about the axis Oz. 

(c) Calculate the volume of the element of this ring contained 
between the planes ¢, 6+ δῴ. 

We take these calculations in turn. 

(a) Using the formula obtained in the preceding Illustration 6, 
we have, for the volume between the two cones, 


27r3{[1 — cos (0 + 68)|—[1 — cos 07} 
= 2mr*{cos 0 — cos (0 + δ0)}. 


O 


Fig. 127. 


* But the timorous reader may proceed straight to the ‘mnemonic’ method 
given on p. 132. ᾽ 
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But, by the mean-value theorem (Vol. I, p. 61), we have the 
ΗΝ cos (0 +80) = cos -- εἴη θ' δθ, 
where @’ lies between 9 and +60. Hence the volume is 

ξπγϑ sin θ' 30. 

(Ὁ) Since 0, 60 are now being kept constant, the volume of the 

‘ring’ contained between spheres of radii r, r + dr is 
an{(r + dr)? — 7} sin 6’ 60 
= 2πίγ3 dr + r(dr)? + 4(dr)3} sin & 68. 

(c) The fraction of the ring cut off between the two planes is 

6¢/27, and so, in all, we have 
Tip = {72 Or +r(dr)? + 4(dr)*} sin 0’ 60 δῴ. 


Inserting suffixes, we may express the summation for the 
integral in the form 1+U+V, 


where I= ΣΙ; 0;, P;) ror, sin 0; 00; δῴ,, 
U=2Xf(r;, θ,, Py) 7(Or;)? sin 0; 00; ODns 
V = 42h (1, 9;, x) (dr,)8 sin 0; 00; Ody. 
The calculation, so far, is exact. 
We prove that U (and similarly, by implication, V) tends to 
zero in the limiting process: 
It is assumed that the given function is bounded throughout the 
volume of integration; suppose that it is numerically less than a 
certain constant K. Then 


U < K=r,(dr,)? sin 6; 60,50, 
numerically. 


Suppose next that, when the subdivision is made, all the elements — 


dr; are less than a certain number p, which itself is to shrink to zero 
for the limit. Then 


U < Kp&r, dr, sin 0; 60,8, 
numerically, where the summation now involves ér;, linearly. 
Suppose further that the volume of integration is bounded in 
extent, so that all points of it lie within a certain distance D of the 
origin O. For summation along a fixed filament near the line at the 


position 0;, ¢,, the value of &r,dr; is certainly less than DXér, 
numerically, or, at the worst, 2D*—since the sum Lér, along the 
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filament cannot exceed the length 2D of the diameter of the sphere 
containing the whole volume of integration. Hence 


Ὁ <2KpD*x sin 0; 30,5, 
numerically. 
Moreover, sin θ᾽ is numerically not greater than unity. Hence 
U <2KpD*X00;5¢,, 
numerically. 


Keeping δῴ, constant, the sum £60; cannot exceed 7, as the radius 
vector θ constant swings from 0 = 0 to 6 = 7. Thus 


U < 27KpD*X6¢,. 


Finally, the whole volume is covered as the plane ¢, rotates 
from ¢ = 0 to ¢ = 27 (or less), so that 


U <2n7KpD*2n, 


U «4π Κρ 
numerically. 
Since K, D are bounded, the right-hand side tends to zero with 
p, so that U tends to zero. 
Tt follows at once that V tends to zero also. 
The summation giving the integral is therefore obtained from 
the expression 


I= >& f(r,9;, $4) 17 6r, sin 0; 00; δῴᾳ, 
(limit) 


or, replacing 0; by 9; in f(r;, 4;, $,)—as is legitimate since (p. 116) 
any point in the element of volume τος may be chosen at which to 
evaluate the function—we have 


I = Ἔ ftp θ᾽, φμ) r; Or, sin 6500; δῴ,. 
(limit) 


This limiting sum is, by definition, the integral denoted by the 
notation 
r= [} f(r,0, $)r2sin 0 dr do dg. 


rsin Odrd0dd 


for the ‘element of volume’ is usually recaptured by using as a 
mnemonic the figure shown in the diagram (fig. 128), which illu- 
strates the surface of a sphere of radius r on which the two ‘circles 
of latitude’ 0, 9+ 60 and the two ‘circles of longitude’ ¢, 6+ δῴ are 
drawn. These delineate the shaded area, which may be regarded 


The formula 
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as a small ‘rectangle’ of sides 7 60 (down the longitude) and r sin 06¢ 
(across the latitude). The element of area on the sphere is thus 


sin θ δῦ δῴ, 
and multiplication by the ‘height’ ér gives the volume 
r?sin 0dr δθ δῴ 


of the element, 


8. The elements by Jacobians; change of variables. We 
proved in the preceding chapter (§$ 5, 6) that, for a transformation 


u=f(z,y), v= g(x, y) 
between two planes in which 2, y and τ, v are taken as Cartesian 


coordinates, a small area A in the (2, y)-plane is transformed into 
a small area A’ in the (w, v)-plane, where 


O(u, v) 
oh? ws + A ’ : 
᾿ aie O(x, y) 
with the similar relation 
ne O(u, v, w) 
ne δία, y, 2) 


for volumes. 
Taking first the case of a plane, consider the evaluation of the 


integral 
| F(x, y) dxdy 
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over a region R. This is effected by means of the summation 


LF (E;, 3) Wy. 
Under the transformation 


U = f(x,y); v= g(x,y), 
the function F(z, y) will become a certain function G(u,v) of u, v; 
and evaluation at the point (£;,7;) will correspond to evaluation at 
a point (A;, “;) of the (ει, v)-plane. The mesh of lines cutting out the 
elements w,,; of the region FR in the (x, y)-plane will give rise to a mesh 
of lines cutting out the elements Q,,; of a region S in the (ει, v)-plane. 
By definition, the integral of the function G(w, v) over the region S 
is just XLG(A;,4,;)Q,;, but this ts Nor equal to the summation 
XP (E;,9;) 0,4 over R. For, though each G(A;,,4;) is equal to the 
corresponding F(£;,7,;), the element Q,, is the multiple + ee) 


O(a, y) 
(evaluated at the point in question) of the corresponding element 


ὡμ- To make the sums equal, we must therefore reduce Q,, in the 


O(u, v) (x, y) 
ratio +1 δίς, y)’ or (p. 98) + 5,0)" Thus 
Os: 73) 


EF (Es) Oy = ΣΘίλρ μὴ | ay μὴ 


so that, in the limit, the relation for the double integrals is 


[freien ff cme 


Finally, although we have given this explanation in terms of 
rectangular Cartesian coordinates, the summations themselves 
are independent of their geometrical interpretation, provided that 
the correct range of the variables τι, v is taken. Compare the 
explanation given for the Illustration 5 on p. 126. 

In particular, if u, v are the polar coordinates 7, 0 given by the 
relation 


D5; 


dudv. 


t= r cos 8, ¥Y¥= rsin8@, 


then O(a, y) " 
o(r, θ) 


οοϑδῦ —rsin@ 


sin@ γοοβθ ee 


so that the formula of transformation is 


| Ϊ ΨΩ χαν « Ϊ Θίν, θ) ἀνθ 


over appropriate ranges of values for r, 0. 
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(Compare p. 124 for the ‘element of area’, rdrd@.) 
Similar considerations apply to volumes. In particular, for 
cylindrical coordinates p, ¢, z, where 


x=pcosd, y=psing, z=2, 


we have ᾿ οοδϑτῴ —psind 0 
(%,Y,2) _ . = 
δ(ρ, 9,2) sin ᾧ amt ᾿ =P» 


and the formula of transformation is 


[| Pe.2dvdyae = [{{| στρ, $,2)pdpagae 
and for spherical polar coordinates r, 0, ὁ, where 


x=rsindcosd, y=rsin@sing, z=rcos8, 


we have 
βίη 9 οοϑῴῤ reosOcosd —rsind sind 
(x,y,z) | .. : ‘ ‘ 
δ», 7, Φ = βίη θ βίῃ ᾧ rcos@sing τ βὶη 9 οοβ ᾧ 


cos 0 —rsin@ 0 
sin@cos¢ cos@cosd —sing 
= ?r*sin@| sin@sing cos@singd cos¢ 
cos 0 —sin 0 0 
= r* sin θίοοβ 8 cos @ - βίη θ sin 0} 


on expanding in terms of the last row. Hence 


δία, ¥, 2) 
o(r, θ, 4) = rein 0, 


a Ϊ Ϊ Ϊ Fie,y,s)dedyde-m | | | G(r, 0, $) r?.sin Odr dO ἀφ. 


(Compare p. 131 for the ‘element of volume’, r?sin Odrd0d¢.) 
The formulae 


{ F(x, y) dady ={[ cu») δ v) 
{I F(x, y,z) dadydz -{{{σω, w) 


it O(%,¥)| 4 
d(u, v) 


O(a, y, 2) 
δίω, υ, w) U, v, w) 


dudv, 


dudvdw 
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enable us to make a transformation from variables x, y, z to new 
variables u, v, w. They are analogous to the formula 


F(x) dx =| Fw) = du 


for a single variable. 
The Illustration which follows shows how such transformations 
may be used to simplify the evaluation of multiple integrals. 


ILLUSTRATION 7. 70 show that the triple integral 


2 
ΤΕΣ 
28 
taken throughout the volume common to the ellipsoid S10 ptat 1 
and the part of the cone a*x? + f*y? = 2*, for which z is pre has 
the value πίαβ. 
We can convert the ellipsoid into a sphere by the substitution 


c=au, y=bv, 2= cw. 


Since a 0 0 
Ot % 2) _)9 & 0|=abe 
O(u, v, w) , 
0 0 ὁ 
we require the integral 


I 24 92 
fi καρ VAM τυ + w*) abedudvdw 


throughout the volume common to the sphere u? + v? + w? = 1 and 
the part of the cone aa*u* + £7b*v? = οξιυξ, for which w is positive. 
Transform now to spherical polars by the substitution 


w=rsindcos¢?, v=rsin@singd, w=rcosd, 
We have 


SIL | (Feast) drdOdg = τ | | | sect tan Odrdodd. 


Keeping @, ¢ constant, integrate with respect to r from r = 0 to 


r = 1, giving 
τὰ Ϊ sec? tan θαθαφ. 


Now for points on the bounding cone, we have 


αϑαϑ sin? 0 cos? + fb? sin? 0 sin? d = c2 cos? 9, 
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Keeping ¢ constant, integrate with respect to 0 from 6 = 0 to 
6 = A, where 


tant = ΤΣ ΤΣ, ΣΤ 
“ἐμ Bm 


05 
=o a*a* cos* ὁ + £76" sin? d ἀφ, 


The cone is traced out by the rotation of ¢ from 0 to 27. Since 
(as will appear at once) we are to divide numerator and denominator 
by cos?¢ to complete the integration, we avoid crossing the 


‘infinite’ values of tan ¢ which will arise by using symmetry in the 


four quadrants and Ce 


. 49 οὗ 
Xai, τ SISTERS 19, 


or, putting tan ᾧ =f, 


a0 Ϊ ο αϑα3- βεὺτιβ 


2a dt 
| πᾷ 
alll 


val 


REVISION EXAMPLES XIII 
University Level 
1, Show that the polar equation of the circle of unit radius which 
passes through the origin and has the initial line as diameter is 


r = 2cos0. 
Integrate the function 1 —(1/r) over the area of this circle. 
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2. A thin lamina covers the positive quadrant of the circle 
x+y? = a® and has density pay at the point (x,y). Find its centre 
of mass. 


3. Integrate sin (x+y) over the square 0 «ὦ «π, 0<y<m7; and 
xe” over the rectangle 0<a<a, 0 «ν «ὃ. 
4, Evaluate | xydxdy 


extended over (i) the triangle of sides x = 0, y= 0, x+y = 3a; 
(ii) the triangle of sides y=0, y=a+a, x+y = 3a; (iii) the 
quadrilateral of sides a=0, y=0, y=a+a, x+y= 3a and 
vertices (0,0), (0, α), (a, 2a), (3a, 0). 
5. Show that the value of 
Ϊ | dxdy dz 
vy (xzy)* 

is four times the volume in the (u,v,w) space which corresponds 
to V, where z=uvww, y=uvw, z= uw’. 


[Joven 


over the positive quadrant of the ellipse x?/a* + y?/b? = 1. 

7. Calculate the volume of the solid bounded by the surface 
z= xy and the planes x = 0, y= 0, a+y =1,z=0. 

8. Evaluate the double integral 


| (x+y)*dady 


taken over the area bounded by the curves zy = 1, αξῃ = 1, and 
the line x = 2. 


9. A plate has a plane triangular base and its thickness at any 
point is proportional to the sum of the perpendicular distances of 
the point (measured in the plane of the base) from the three sides 
of the triangle. Prove that its mean thickness is equal to the 
thickness at the centroid of the triangle. 


10. Evaluate the integral 
| (x+y+a)dxdy 


taken over the circular area 2? + y?<a?. 


10 M Ir 


6. Evaluate 
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11. Evaluate the integral 


[fs 


2 42 
where - = = +55 and the region of integration is the elliptic area 


al a 
mt pS! 


12. If D is the interior of the closed curve whose equation, in 
ordinary polar coordinates, is r = 2a(1+ cos @), evaluate the double 


ss | Ϊ » (a? + y*)# ἀχάν, 


where x = rcos@, y = rsin@. 


13. A bounded function f(x,y) is defined in a certain simple 
closed region R of the x, y plane. Explain briefly what is meant 
by the double integral 


Ἐπ | | Jeuddedy. 


The variables x, y are expressed in terms of two new independent 
variables τὸ, v so that f(x,y) = F(u,v) and R becomes a region S of 
the u, v plane. State, without proof, the expression for J as a double 
integral over 9. 

If Ris the region bounded by the four parabolas y? = ax, y? = Az, 
x* = by, α = By, where 0<a< A, 0<b<B, find new variables wu, v 
such that § is a rectangle, and hence show that the area of R is 


}(A—a)(B—b). 


14. A solid sphere is bounded by the surface x?+y?+2? = a?, 
and a cylindrical hole is drilled through the sphere, the boundary of 
the hole being part of the surface x*+y?— bx = 0, where 0<b<a. 
Express the volume V of material removed as a double integral 
taken over a suitable region of the zy-plane, and derive the formula 


V = al (a? —r?)t γάγαθ, 
where this double integral is taken over the region defined by 


0<0< hn, 0<r<bcosf. 
Find V when ὦ = a. 
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15. Show that 


[J ewer+yydedy = 5 [ἰώ αι, 


where the double integral is taken over the area 
20, y>0, w+y?<1. 
Evaluate the integral 


Ϊ ny(b?x? + a2y?)tdady 


over the area x>0, y>0, bx? + a®y? < αϑ3, 


16. Ifa, y, z are the perpendicular distances of a point from the 
sides of an equilateral triangle of side 2a, prove that the mean value, 
with respect to area, of xyz taken over the interior of the triangle 
is a5 /3. 

17, A uniform solid of density p is in the form of an anchor 
ring, generated by rotating a circle of radius a about a line q in its 
plane and distant / (>a) from the centre of the circle. Show that 
the moment of inertia of the ring about q is 27%pa*l(l? + 3a*). 


18, Evaluate the integral 
| Ϊ (ax* + by? + cz*) dadydz, 


taken throughout the spherical volume x? + y?+2z2<1, 


19. Prove that 
(i) the area of a plane sector extending from the origin to a 
curve 2 = f(u), y = φ() is 
co UY 
1 
5 | de ay |éus 


»"--- 


du du 


(ii) the volume of a cone extending from the origin to a surface 
a= f(u,v), y = d(u,v),2= (u,v) is 


Zs ¥ 8 

rrr) oe ὃψ & 

All Cu ou du |dudv. 
Ox Oy 0% 
ov dv ov 
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20. Express the integral 


ftav{{” seve 
in polar coordinates. . tia 


Evaluate the integral when 
fey=Z—% 
? e+ y? Ξ 
21. Give an account of triple integrals sufficient to enable you 
to find the moment of inertia of the ellipsoid 
ἃ 98 
about the axis y = z = 0. 


22. Prove that the volume in the positive octant which lies 
inside the surface y?+2? = 4axz and outside the surface y? = 4k*ax 
(where &< 1) and between the planes x = 0, x = cis 

ac*{cos-1 k —k ./(1—k?)}. 

23. If γι, 7ᾳ are the distances of the point (z, y) from the foci 
(+c,0) of the ellipse 2?/a?+y?/b? = 1, show, by means of the 
transformation x+iy = ccosh (u+iv) or otherwise, that 


Ϊ Ϊ [ 87) ἀνάν anh 
5 κ΄ 
the integral being taken over the interior of the ellipse. 


24. By transforming to polar coordinates, show that the 
integral of the function (x? + y?)?/(ay)?, taken over the area common 
to the circles 2? + y? = aw, 2*+y = by (where a> 0, ὃ.» 0), is αὖ. 

Verify the result by transforming to variables w, v given by the 
relations wx = vy = 2*7+y". 

25. A uniform plane lamina, of density p per unit area, is formed 
by the area common to the ellipse x?/a? + y?/b® = 1 and the parabola 
2ay* = 32x. Find the moment of inertia of the lamina about the 
axis x = 0. 


26. Show that the value of 
| εἰν--αὐψ τα) da αν 


over the triangle bounded by the three lines x = 0,y = 0,x+y=1 
is equal to (e? — 1)/4e. 


REVISION EXAMPLES XIII 141 
27. Prove that the quadric 1+z = a*x?+ by? meets the sphere 
a +y*+2? = 1 in a curve the polar coordinates (r,0,¢) of any 
point of which satisfy the equations 
(a?—1)+(b?—1)tan?¢ 
a*+b* tan? d 3 
and prove that the volume common to the sphere and the cone 


joining the origin to this curve is 277/3ab. It is assumed that a, ὃ are 
each greater than unity. 


28. Prove that the moment of inertia of a uniform lamina in the 
form of a limagon r = a+6cos@ (where a>b) about a line through 
the origin perpendicular to its plane is 


(8a4 + 24α35 + 364) M 
8(2α5- 8) ᾿’ 
where ΤΠ is the mass of the lamina. 


29. Show that the volume common to the sphere x? + y?+ 2? = r2 
and the ellipsoid 2a* sin® α + 2y? cos® α +-2? = γ3, where 0 «α <4z, is 


$13 (47 + 2a cosec 2a). 


r=1, 0050 
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CHAPTER XVIII 
THE SKETCHING OF CURVES 


It is useful to be able to indicate by a rough sketch the general 
shape of a curve of given equation 


f(@,y) = 0, 


and the aim of the present chapter is to show, by means of particular 
examples, how to set about the problem. Details of the theory must 
be sought in a text-book on plane curves. 

We saw in Volume I (p. 59) what features are most helpful 
when the equation can be reduced to the simple form 


y = (2), 


where g(x) is a function of x only. In the more general cases, with 
which we are now concerned, symmetry, gradient, concavity and 
so on—especially, perhaps, symmetry—form valuable guides, but 
more detailed analysis becomes necessary. 

The following examples illustrate curves which the reader should 
be able to draw freely with little trouble, and they will be used 
without further explanation when required later. They are all of 
the form y = Ax", where x is a rational number, and so can be 
sketched by methods already given. 


EXAMPLES I 
Sketch the curves given by the equations: 


l. y=2%. 2, y= 2, 
3. y®? = 2%, 4. yma, 
5. yx? = 2, 6. ¥ = κα", 
1. y= 2°, 8, ἢ = 2, 
9. =z, 10. 9° = 24. 
11, y= 2%. 12, y® = 2°, 


We confine our attention to curves for which the function 
f(x,y) is a polynomial in x, y. Such a curve is said to be ALGEBRAIC. 
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1. Towards a general method. Though our interest is centred 
in curves whose equation f(x,y) = 0 cannot be solved for y in terms 
of x, the general method may be illustrated by an example in which, 
exceptionally, the solution can be effected. This particular curve 
has what is called a ‘singular point’ at the origin, and we shall have 
to show later how such points may be located more generally. 
Apart from this, our treatment follows a standard process involving 
three stages: 
(i) examination of the curve near the origin (or, if necessary, 

at other singular points); 

(ii) examination of the curve at ‘distant’ points; 

(iii) use of symmetry, special points (for example, points where 
the curve crosses the axes), ideas of continuity, and so on, to 
complete the curve. 


ILLUSTRATION 1. 70 sketch the curve 
ψ = x + 423, 
The solution for y in terms of x is immediate, giving the two 
PONE elite, ἀν - αἰ νυ 


The curve is symmetrical about the x-axis; the positive sign corre- 
sponds to points ‘above’ the axis and the negative to points 
‘below’. It will therefore be sufficient to consider the former, 


namely, y= a(1 ΕΞ 4a)t, 


(i) Small values of x. When z is small, an approximation to the 
value of y may be obtained by means of the binomial theorem. Thus 


y=x(1+2x+...) 
2+ 2e2+..., 
The first approximation y=2 


shows that the curve lies very close to the line y = x, which is a 
tangent at the origin. The next approximation 


Y=x+ 2x? 


then shows that the value of y for the curve near the origin exceeds 
the value of y for the tangent. Hence the curve lies above the 
tangent, as the diagram (fig. 129) indicates. 
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Fig. 132. 


ee ee ee ~ 


eg οὧὧἱἱν “ὦ ὦ “ ——, ὦ ἀπ σα 28 se 06 


146 THE SKETCHING OF CURVES 
(ii) Large values of x. When z is large we again use the binomial 
theorem, but this time to express y in descending powers of x. Thus 


1\% 
= 273 vs 
ψ 514 Ὁ) 


= 2214+...) 


Qo +daot+.... 
The first approximation y=2z4, 
or y= 22 γα, 
shows that, for large values, x is necessarily positive, and the curve 
behaves very like the simpler curve ψ = 423 shown (for the 
‘upper’ half) in the diagram (fig. 130). 

(iii) Intermediate values. We have now, using symmetry, reached 
the stage indicated in fig. 131, and our problem is to ‘join up’ the 
ares to complete the whole curve. The shape is made clear by noting 
that the curve meets the x-axis where x2(1+ 4x) = 0, that is, in the 
origin and in the point ( -- ξ, 0); and the y-axis where ψ = 0, that is, 
in the origin only. We therefore obtain the form given in the 
diagram (fig. 132). 

[It must be admitted that the final steps involve an appeal to 
a certain amount of intuition, and it is an intuition which the 
reader should develop. There is no other reasonable way of joining 
the ares, given that the axes may not be crossed except at the point 


(—4, 0).] 
EXAMPLES II 
Sketch the curves given by the equations: 
1, y* = α-- δ, 2. y® = χ3-- 4.χ3, 
3. α τὸ y?+ 4y3, 4, y= y?—4y/4, 
5. y?—2ay+z° = 0, 6. ψ5- 4χυ -- 4χ = 0. 


2. The method for small values οὗ χ. We now seek a method* 


οὗ ‘successive approximation’ to be used when simpler devices, 


such as the binomial theorem of the preceding illustration, are not 


* We attempt to give a general method applicable to the many and varied 
complications. This one method, when correctly grasped, can be applied very 
widely, and seems to require a minimum amount of specialized information. 
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available. We give a number of examples to show the characteristic 
features which arise most frequently. In exposition, the first 
example will be given in detail, but after that less will be said in 
explanation as the process becomes more familiar. 

The basic property, on which the method depends, will be seen 
by returning for a moment to the preceding Illustration 1 (i) 
(p. 143), where we obtained the successive approximations (for 
small 2) 

Y=x, 
y=x + 2.3, 


The first approximation is linear, corresponding to the tangent 
y = x; the next approximation is obtained by adding to the first 
a term of higher degree. This is what we seek to do in the examples 
which follow, where each approximation is obtained from its 
predecessor by bringing in further powers of x. 


ILLUSTRATION 2. Το sketch the curve 
2x—y+x*?+ 3xy—y> = 0 
for small values of x. 
Near the origin the terms of lowest degree are most important, 
and we have the first approximation 


2x —y=0 
or , y= 22. 
The next approximation is therefore of the form 
y= 2x+Azx" (n>1), 


where A, ἢ are constants to be determined, and where ἢ is greater 
than 1 since the term Az” has to be small compared with 2a. Sub- 
stituting this approximation in the equation of the curve, we have 


the relation 
— Ax" +x? + 3x(2a+ Ax”) — (2. -- Ax”) =O, 


We are concerned with small values of x, so that terms of lowest 
degree are most important. In particular, in each of the expressions 
in brackets, namely, 2x+ Aa", where n>1, the term Az” is neg- 
ligible compared with 2x and may therefore be omitted; as a point 
of technique, we always imply such omissions by dots, thus: 


— Ax” +27 + 3a(2a+...)—(27+...)8=0, 
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We then see that the final bracket, with terms in x*, may be neglected 
entirely in the presence of earlier terms in x®. We therefore have 


th on 
e position — Ax” + 722+ ...=0, 


All the other terms have been omitted by comparison with these; 
indeed, the comparison has been with the term 7x", which is of 
lower degree than anything left out. 

The essential point of the process now arises, namely, that, by 
choosing the particular values A = 7, n = 2, we reach a relation 


Fig. 133. 


in which the left-hand side has no terms involving z to the power 
2 or less, so that it approximates very closely to the (zero) value 
on the right. In other words, the values A = 7, n = 2 yield the 
approximation y= 20+ Ta? 


fitting the given equation of the curve to a high degree of accuracy. 

It follows that the curve at the origin lies very close to the line 
y = 2x, which is the tangent there; also, since 72? is positive, the 
value of y for the curve exceeds the value of y for the tangent. 
Hence the curve lies ‘above’ the tangent (fig. 133). 
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ILLUSTRATION 3. (Double point.) To sketch the curve 
2.5 — 3ay+y?+a°y+y>? +24 = 0 


for small values of x. 


The terms of lowest degree give the approximation 
22° — 3xy + y?=0 
or (y—2) (y—2a)=0. 
There are therefore two possibilities: the curve may behave very 


like the straight line y = x, or very like the straight line y = 2z. 
The curve has a DOUBLE POINT (fig. 134) at the origin, with these 


lines as the tangents. 
We take the two cases in turn: 
(i) When y=, 


the next approximation is 
y=2%+Ax™ (n>1). 
Experience shows that, when the terms of lowest degree factorize, 


the approximation should be inserted in the factorized form of 
equation. Thus, the equation is 


(2 -- ψ) (α -- 2) ὁ αῦν τ νὴ +2x* = 0, 
and the corresponding relation is 
(w— Ax) (— Ax”) +2°(2 + Ax”) + (2+ Aa”)? + a4=0., 


Omitting terms involving powers of x which (remembering that 
n> 1) are seen as negligible, we obtain the relation 


(w—...)(—Aa™) + 27(a+...)+(e+...)?+...=0, 
or — Ag"! + 243 +... 0. 
For closest approximation we cancel the terms of lowest degree, 
toning A=2, n=2. 
Hence the second approximation to the value of y is 
y=2+ 22%, 
showing that the curve lies ‘above’ the tangent. 


(ii) When y= 22, 
the next approximation is 
y=2a4+ Ax. 
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(It is convenient to use the same letters A, n, though, of course, 
they now have fresh values.) 

Substituting and omitting obviously negligible terms, we have 
the relation 

(—Aza”) (—a—...)+29(2a+...)+(2e4+...)®?+...20, 
or Ax) + 10a3+...20. 
For closest approximation, we have 
A=-10, n=2, 

giving y= 22 — 102, 
so that the curve lies ‘below’ the tangent. 

The shape of the curve near the origin in indicated in fig. 134. 


Fig. 134, 


ILLUSTRATION 4. (Cusp.) Το sketch the curve 
4a? — ἀφ + y? + 16x*y + 4y3 +24 = 0 
for small values of x. 

(The characteristic features of this equation are that the terms 
of lowest degree are quadratic, forming the square of a linear 
function which is not a factor of the cubic terms; that is, 27 —y is 
not a factor of 16x*y + 4y°.) 

The terms of lowest degree give the approximation 


4a? — day + y?=0 
or (y—2x)?=0, 
so that y= 22. 
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The next approximation is 

y=22+Ax"™ (n>1). 

Substituting and omitting obviously negligible terms, we have 
A*x?” + 1629(2a 4+...) +4(2a+...)8+...0, 

or A279” + 6493 + ...0, 
For closest approximation, we choose A, ” so that 

A?+64=0, 2n=3, 
or A=48,/(-1), n=. 


Fig. 135. 


At first sight the presence of γί -- 1) seems alarming, but it is 
essential to the argument. For the approximation is 

y= 20+ 8 γ{-- 1) αἕ, 
or y= 2x + 8x ./(—2). 

It follows that, for small values, 2 is necessarily negative; and 
also that the curve lies on either side of the tangent y = 2x, one 
‘branch’ being ‘above’ it and the other ‘below’, as in the diagram 
(fig. 135). 

At such a point, the curve is said to have a cusp. 

Inuustration 5. (T'acnode.) Το sketch the curve 

4.3 — day + y* + 3. — xy + 2ay?—y3 — 12a4+y4 = 0 
for small values of x. 


152 THE SKETCHING OF CURVES TACNODE 153 


(The characteristic features of this equation are that the terms 
of lowest degree are quadratic, forming the square of a linear func- 
tion which is also a factor of the cubic terms; that is, 2x—y is a 
factor of 2a3 — ay + 2ay?— ἡ = (2. —y) (x? + y?).) 

Bearing in mind the remark on p. 149, we use the ‘factorized’ 
form of equation 

(20 —y)? + (2e—y) (αὐ + ν᾽) — 1224 Ὁ γ᾽ = 0. 
As in the preceding example, we consider the approximation 
y=22+Aa™ (n>1). 
Substituting and omitting obviously negligible terms, we have 
A%y2” — Ax {x? + (2a+...)*}-—12at+ (2e+...)#=0, 

or A*yin — § Ag+? + 4at+ ... 0. 
We have retained all terms which might conceivably be wanted; 
the others are seen at once to be negligible by comparison with 
their immediate neighbours. . 

It is now clear that the value ἢ = 2 involves all three terms, so 

that we have, for closest approximation, 
A?—5A+4=0, n=2. 
Thus A=lor4 n=2. 
We therefore have the two approximations 
y=2e+2%, 
y= 2x + 42, 
The curve thus consists of two touching branches, each lying 
‘above’ the tangent on either side of the origin (fig. 136). 
At such a point, the curve is said to have a TACNODE. 
We conclude this paragraph with an example, not involving any 


essentially new principle, to demonstrate the slight modification 
necessary when the tangent at the origin is the y-axis: 


Fig. 136. 


Fig. 137. 
InLustTRATION 6. 70 sketch the curve 


at+et—yt+ys = 0 (Note that πὶ is greater than 1, not 0; the first approximation was 


x=0.y.) 
On substituting and omitting obviously negligible terms, we have 
A*yn +... —yt+...0, 


or Aty2n — yt + ...=0, 


Ir M III 


for small values of x. 

The approximation z?= 0 suggests that we interchange the roles 
of x, y and attempt to express x in increasing powers of y. We thus 
make the approximation 

zzAy” (n>1). 
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For closest approximation, we choose A, so that 
At*=1, 2n=4, 
or A=+1, n=2. 
Hence we have the two approximations 
aay’, az—y’, 
leading to the form (a tacnode) shown in the diagram (fig. 137). 
GENERALISATION. A point P on a curve given by an equation of 
degree k in x, y is called a MULTIPLE POINT, of multiplicity p, if an 
arbitrary line through P meets it in only k—» points other than P. 
It follows at once that, if P is the origin, the equation contains 
no terms of degree less than p and at least one term of degree p. 
(Compare Illustrations 3, 4, 5, where k = 4, p = 2.) 


EXAMPLES III 


Use the method of this section to obtain the approximations near — 


the origin for the curves given in Examples II (p. 146). 


3. The method for large values of *. Very similar methods 
may also be applied to study the form of a curve at large distances 
from the origin. Referring again to the preliminary Illustration 1 (ii) 
(p. 146), we recall that we were able there to express y in decreasing 
powers of x; we had the successive approximations 

yx=2x8, yx 2ah+1t. 

In the more general case, we obtain our first approximation 
from the terms of highest degree in the equation f(x,y) = 0. The 
subsequent improvements are obtained just as in the case of small 
values of x. The Illustrations which follow exemplify a number of 
typical curves. 

(Note. The calculations will of themselves justify the use of the 
terms of highest degree for the first approximation; when 2, y are 
both large these terms have obvious importance. But it is possible 
for, say, x to remain finite while y becomes large, as in the familiar 
example of the rectangular hyperbola zy = 1. We hold this case 
back for some time (pp. 172-3).) 


ILLUSTRATION 7. 70 sketch the curve 
208 — αν — 2ay? + y® — 2.3" — Qay + 2x+ ὃν = 0 
for large values of x. , 
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When z is very large, only the terms of highest degree are sig- 
nificant, so we have the first approximation 


2x3 — αἷρ — Qay? + y?=0 
or (w—y) (2a—y) (w+y)=0. 


We therefore consider in turn three approximations to the value 
of y, given by y=a, y=2x, ψ- --α. For convenience, we shall 
complete the calculations for all three cases before commenting 
on the geometrical interpretations. 

(i) For y=, the next approximation is 


y=x+ B2?, 
where Ὁ is less than 1. Substituting this value in the equation 
(x—y) (2a —y) (w+y) — 2 -- 2xy + 25 -ἰ ὃν = 0, 
where the terms of highest degree are expressed in factorized form, 
we have — Ba? (x — Bx?) (2+ Ba?) 
— 2.3 — 2χί(α -ἰ Bx?) 
+ 25 + 8(x+ Bu?)=0, 


We are now interested in the larger powers of x, so that we may 
omit terms of lesser degree. Remembering that p< 1, we thus have 


— Bu?(x—...)(2u+...)—2u*— 2u(v+...)+...20 
or — 2Burt?—4o2+...=0, 


each omission (indicated by dots) being justified on comparison 
with neighbouring terms. 

For closest approximation we choose B, p so that these terms 
vanish, giving Bu—2, 9 =; 


Hence the approximation is 
y=x—2, 


(ii) For y= 2z, the next approximation is 
y= 24+ BueP (p<l) 


(for fresh values of B, p), so that, on substituting and keeping terms 
of highest degree, 


(—x—...)(— Ba?) (8a+...) — 2.3 — 2a(2a+...)+...20 
or 3BuPt? — 62?+...<0, 


so that B=2, p=0. 
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Hence the approximation is 
y= 2x + 2. 


(iii) For y= —z, the next approximation is 


y>=—2+BaP (p<l), 
so that, as before, 


(2a—...) (Ba—...) (Ba?) — 2. -- 2a(—a+...)+...20 
or 6 Bart? +0 .a7+...20. 
Strictly speaking, B is essentially non-zero, so we ought to con- 
tinue the approximation by retaining further terms. But com- 
parison with the two preceding approximations shows us that the 
vanishing of the coefficient of x is, in a sense, ‘accidental’. We are 


about to examine all three cases further, so, for uniformity of 
treatment, we pause here and keep the approximation in the form 


y= —2. 

We have therefore reached the three approximations 

y=x—2, 
y= 2x + 2, 
y= --. 

Experience shows that the next approximations are best ob- 
tained by replacing the terms of highest degree (x — y) (2a —y) (w+y) 
in the given equation by the corresponding product 

(%—y — 2) (2e—y+2) (x+y) 
and making the consequential adjustments. Now we have, by 
direct multiplication, the identity 

(w—y — 2) (2 --τῖ| +2) (u+y) 

= (ὦ -- ψ) (2x—y) (w+y) — 2a" — xy — 4α — ἄν, 
so that the equation of the curve is expressible in the form 
(v—y—2) (2a—y+2) (x+y) + 6x+12y = 0. 

Consider the three approximations in turn: 

(i) Let y¥=x—2+ χα, 
where q is now less than 0. On substituting and keeping only terms 
of highest degree, we have 


(—Cx2) (2+...) (2a—...) + 6x +12(x—...)=0, 
or — 2Caxt? + 182+...=0, 
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so that, for closest approximation, 


C=9, g=-l. 

Hence γα 265, 

(ii) Let y=2x+2+Cx% (¢<0). 
Then (—2-—...)(—C2®) (8a+...)+6a+12(2e+...)=0, 
or 3C x22 + 302+...0, 
so that C=-10, g=-1. 
Hence γαϑω 2- Ὁ, 

(iii) Let y= —x+Cxt (g<0), 


where we know from the earlier work that g must be less than zero 
although, in this particular case, the constant term is absent. 
As before, 

(2a—...)(3a—...) (Cx%) + 6x+12(—2+...)=0, 


or 6Catt? — 6a +...=0, 
so that C=z1, g=-—l1. 
Hence er reo 

© 


To summarize, we have obtained the three approximations 


9 
YAt—2+-, 


The geometrical significance may be appreciated most readily 
by referring to the sketch (not drawn to scale) in the diagram 
(fig. 138). 

When z is large, y is very nearly equal to one or other of the 
expressions x — 2, 2x-+ 2, —a, and so the curve (for large values of 2) 
lies very near to the three straight lines 


y= x-2, 
y = 2x+2, 


y= -- α. 
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Straight lines which the curve approaches in this way are called 
ASYMPTOTES of the curve. 

The further approximations then show that, if x is large and 
positive, the curve lies ‘above’ the asymptote y = x—2, ‘below’ 
the asymptote y = 25: - 2, and ‘above’ the asymptote y = —a; if 
x is large and negative, the curve lies ‘below’ y = x—2, ‘above’ 
y = 2x+2, and ‘below’ y = —x. These features should be checked 


by reference to the diagram. 


Fig. 138. 


InLustTRATION 8. (Parallel asymptotes.) Το sketch the curve 
x + x*y — xy? —y* + bay + by? —19%—lly+6=0 


for large values of x. 
The terms of highest degree give the approximation 


αϑ + xy — xy? — y= 
or (x—y) (w+y)*?=0. 
We therefore have two approximations 
y=, 
yee ie 


of which the second may be expected to present fresh features; we 
begin with the normal case y=2. , 
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Taking the equation of the curve in the form 
(~—y) (v+y)* + bay + 6y?—19%—1ly+6 = 0, 
we make the approximation 


y=x+ BaP (p<1), 
so that 
(—Ba?) (2a+...)?+ 6a(a+...)+6(a+...)?—...0, 
or —4Bypt2 + 1227 + οὔθ, 
Hence B=3, p=0O, 
and the approximation is y=2+3, 
Suppose next that 
y= —2+BxP (p<1). 


We notice that the terms of degree one less than the highest 
(that is, the quadratic terms) also contain x+y as a factor—a 
point that should always be investigated when repeated factors 
occur—and so we write the given equation in the form 


(w—y) (2+y)* + 6y(%+y)—19%—1ly+6 = 0. 


Application of the routine process of substitution needs care. 
At first sight, we have 


(2. --...) (Bx?) + 6(—a2+...) (Bu?) —...=0, 
or 2B%z2P+1 — 6Barti+ ...=0, 


so that p = 0. But a glance at the equation then shows that the 
terms —19x—1ly are also involved, and we therefore incorporate 
them in the relation 


(2. —...) (B®) + 6(—a +...) (Ba?) —19"—11(—a+...)+...=0, 


where the terms now omitted are all negligible by comparison with 
their immediate neighbours. Thus 


25%x°P+1 — § Berti — 8r+...=0, 


and we take p=9, 
285 -- 68 -- 8 = 0, 
or B=-1, +4, 
giving the two approximations 
y=—«-l, 
y= —x+4, 
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We have obtained in all an asymptote 
y=2+3 
and two parallel asymptotes 
y= —2-1; 
ψτ- --,χχἹ 4. 
For the next approximation, we begin as before (p. 156), with the 


identity (obtained by direct multiplication of the three expressions 
derived from the approximations already obtained) 


(ν --α, --  33ξ (ψ -Ἐαι - 1) (ν Ἐπ -- 4) 
ΞΞ (υ -- “) (y+2x)*—6xy -- 6y? + 135 - 5y + 12, 
so that the equation of the curve is 
(y—x—3) (ν Ἐα - 1) (ψ Ἐπ -- 4)-- θ(α - ψ) -- 18 = 0. 
(Note that we have again exhibited the factor x+y in the term of 


highest remaining degree.) 
The approximation 
y=2+3+Cxt (¢<0) 


gives Cx 2x +...) (2+...) +6(2e+...)—...=0, 
or 4C xt? 4+ 124+ ...0, 
so that C=-3, q=-l, 


and the approximation is 
3 
id tha 
The approximation 
y= --,ο--Ἰ - τ χα (ᾳ «0) 
gives (—2a+...)(Cxv%)(—5+...)+6(—1+...)—18=0, 
or 10Cxt+1-— 24+ ...=0, 
so that C=32, q=-1, 
and the approximation is 
Ξε: gi ee 
δα" 
The approximation 
y= —2+44+Cx% (¢<0) 


gives (—2x—...) (5+...) (χα) +6(4+...)— 180, 
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or —100xtt1+6+...=0, 
so that C=2%, qg=-1l, 
and the approximation is 


3 
y= —2+4+ Ἐπ" 
We have thus obtained three asymptotes, of which two are 
parallel; the further approximations show that (i) when 2 is large 
and positive the curve lies ‘below’ y = «+3 and ‘above’y = -- ᾽ἡ: -- ἰ 
and y = —a+4; (ii) when z is large and negative the curve lies 
‘above’ y = 2+3 and ‘below’ y = —a—1 and y = —x+4. These 


Ny 
Ζ 


Fig. 139. 


properties should be checked by reference to the diagram (fig. 139) 

which exhibits the main features of the curve but is not drawn to 
scale. (Note as ‘guides’ that the curve meets the y-axis in the 
points (0,1), (0,2), (0,3), that it crosses the asymptote y = x+3 at 
the point (0, 3), and that it does not meet the parallel asymptotes.) 
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EXAMPLES IV 
Use this method to find the asymptotes of the following curves: 
1, 2a?—-3ay+y?+3x—2y—1=0. 
x? —y?+ 2e—2y+1=0. 
(x+y) (ὦ - ψ) (2a —y) + 42? — Qay + 2y?- 2y+2+7 = 0. 
(v+y)(e—y)?+2*—ay—d3y = 0. 
. (w—y) (e+y)? + 40? + 4ay+1 = 0. 


InLustRATION 9. (Parabolic asymptotes; also, completing the 
sketch.) Το sketch the curve 


xv — ay —ay?+y>—2x*+ 3xy = 0, 
The terms of highest degree are 
(y+) (ν -- α), 


so that, for large values of x, we have to consider the two approxi- 
mations y= -- ἃ, 


oR bo 


Y=. 
The first of these follows standard pattern, and we may deal with 
it quickly. The equation is 

(y+x) (y—x)?—2* + 3ay = 0, 


so that the approximation 
y= —2+BuP (p<1) 
gives (BaP) (—2x+...)?—2*+ 3a(—x+...)=0, 
or 4ByrPt2— 4y2+ ...0, 
so that po=0, Bel, 
and there is therefore an asymptote 
y=-—2z+1. 


The given equation, in the form 
| (y+%—1)(y—2)?+y(x+y) = 0, 
leads from the approximation 


y= —2+1+Cx% (¢<0) 
to the relation 


(C2) (—2a+...)?+(—a+...) (1+...) =0, 
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or 4Cxtt? a+ ...20, 
so that we take q=-l1, C=}, 
ΕΝ i 
giving y= —x+1 Lor 


Thus the curve lies ‘above’ the asymptote when 2 is positive and 
‘below’ when z is negative. 


The approximation Y=2 
introduces new features, so we treat it with more detail. Writing 
y=x+ BaP, 
Oo x 
Fig. 140. 


we have the relation 
(22+...) (Ba?)?—2*+ 3a(x+...)=0, 


or 2Br%y*Pt1 + 2724... =0, 

so that p=%4, Β-- τὶ (¢=¥(-1)). 

Thus y=xt./(—2), 

which shows that x must be negative on this branch when numeric- 
ally large. 


The curve this time does not approach any straight line. On the 
other hand, it behaves like the two branches 
y=2a+/(-2), y=2—(—2), 
which, when considered together, appear as the two ‘arms’ of the 
parabola (fig. 140) (y—x)?+a = 0. 
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Before commenting further, we proceed to the next approxima- 
tion, writing the equation of the curve in the form 


(y +2) {(y—«)? +2} — 22" + 2xy = 0, 
or (y+%) (ν --᾿α᾽ εὖ να) (y—x—t /x) + 2a(y—z) = 0, 
The approximation 
πα τ εὺ απ χα (q<}) 


gives (2a +...) (2ὲ γα -" ...) (Cx%) + 2a(t Ja+...)=0 
or 4iCattt + Qiak +... <0, 

so that q=0, C=-}. 

Hence y=x+t α-- ἢ, 


with a similar result when i is replaced by —i. Thus a very close 
approximation to the curve for large values of x is found in the 


δος, (y—2+4)?+2 τ 0, 


which is called the ΑΒΥΜΡΤΟΤΙΟ PARABOLA of the curve. 

The earlier parabola, (y—«)*+2 = 0, gives a good idea of the 
general appearance of the curve, and is usually found sufficient in 
practice. When more accurate ‘placing’ is required, the asymptotic 
parabola must be used. 

This completes the investigation of the curve for large values 
of x. In order to sketch the whole curve, we now look at the small 
values. The terms of lowest degree give x~0 and y=ja, so we 
examine each of them in turn. 

For x0, we proceed in ascending powers of y and consider 
(compare p. 152) the approximation 


eAy”™ (n>1). 
Substituting in the equation (arranged now in terms of ascending 
degree) “(3y— 2x) + 2° — ay —ay?+y3 = 0, 
we have Ay" (3y—...)+...+¢y3=0, 


retaining terms of lowest degree in y. Thus 
3Ay"t1+ 3+ ...=0, 
so that, for best approximation, we have 
m=2, A=-}, 
Hence = — ἐγξ, 
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showing that, near the origin, the curve lies to the left of the 


tangent x = 0. 
Similarly, if y=}x+Anu" (n>1), 
then 2(3Aa")+29—a(4e4+...)—x(4o+...)?+(fe+...2=0 
or 3Aanrt + 38a3+...0, 
so that n=2, A=-H1, 
Thus y= }r—167?, 


showing that, near the origin, the curve lies ‘below’ the tangent 
y = $e. 


Fig. 141. 


We have therefore reached the position given in the diagram 
(fig. 141). 

For the rest of the curve we must revert to more general con- 
siderations. Unfortunately we receive no help from symmetry, but 
we are able to calculate the coordinates of the points where the 
curve crosses the axes and the linear asymptote. 

When x = 0, we have y = Ὁ only; when y = 0, we havez = 0,0, 1; 
when x+y = 1, we have x = 1, y = 0. Hence the curve does not 
meet the y-axis except at the origin, and it cuts the z-axis and 
the asymptote where they cross at the point (1,0). Thus we obtain 
the form in the diagram (fig. 142), which, however, is not drawn 
to scale. 
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VN 


Fig. 142. 


EXAMPLES V 

Find the asymptotic parabola of the curves: 

1, 2(2y—2) (y+2)*?+ 92y = 0. 

2. y® = 4χϑ(χ -" υ). 

3. x(y+3x)? = y?+9. 
ILLUSTRATION 10. (Loop and inflexion.) To sketch the curve 

x—yi+a7%(a+y) = 0. 
For large values of x, we have 
(x —y) (αὐ + αν + ay? + ay® + y*)=0, 
and the only real approximation is 
Y=. 


(The presence of terms not resoluble into real linear factors may 
indicate the presence of a loop. But this is not necessarily the case.) 


Writing y=x+BaP (p<1), 

we have (— Bax?) (5at+...)+27(2a+...)=0, 
or — 5Burt4+ 2a3 + ...=0, 

so that B=%, p=-1. 
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We therefore have the approximation 


2 
YLT ~,» 


showing that the line y = xis an asymptote, and also that the curve 
lies ‘above’ it when is positive and ‘below’ it when z is negative. 

When z is small, we have the two approximations 2*=0, x +y=0, 
which we consider in turn. 


If zzAy” (n>1), 
then —y> + Ar*y*"(y+...)=0, 
or — y+ Atymti +... 0, 
so that n=2, A=+1. 
We therefore have the two approximations 

zany", 
= —y?, 


indicating the presence of a tacnode (p. 153). _ 
If =-—2+Anz" (n>1), 


then x'—(—x#+...)§+27(Aa")=0, 
or 245 + Agrt2+ ...=0, 
so that ; n=3, 4-π η - 2. 


We therefore have the approximation 
y= —x— 22%, 


Thus the curve lies ‘below’ the tangent y = —x when z is positive 
and ‘above’ when =z is negative. 

When the curve crosses the tangent in this way, the curve is 
said to have an INFLEXION at the point. 

This particular curve has both a tacnode (tangent x = 0) and an 
inflexion (tangent x+y = 0) at the origin. 

We have so far reached the state indicated in the diagram 
(fig. 143), and our problem is to combine these arcs into a single 
curve. 

We note that, since replacing x, y by —2, —y does not affect the 
equation, the curve is symmetrical about the origin. Also it meets 
the axes in the origin only, the asymptote x = y in the origin only, 
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and the line x+y = 0 in the origin only. Thus the curve does not 
cross any of the lines shown in the diagram, and so its form is that 
indicated in fig. 144, The arcs within the area bounded by the lines 
« = 0,x+y = 0 join up by means of loops. 


4. To locate the multiple points of the curve f(x, y)=0. 
Suppose that the point (ἢ, /), initially unknown, is a multiple point 
of the curve Ἠξὼ αὐ 


Referred to axes through that point parallel to the given axes, the 
equation assumes the form 


f(z’ +h,y' +k) = 0, 
where f(h, k) is, of course, zero. Expanding by Taylor’s theorem 
(p. 57), we have 
MP LD a 
f(h, k)+ (- mt 3 δ f(h, k) +... = 0, 
Of(h,k) , δ (ἃ, kb), , Pf(h, k) 
oh” * Ok 4 tot ane 


By hypothesis, the new origin is a multiple point of the curve and 
so the coefficients of x’, y’ both vanish. Thus 


Of(h,k) _ 0 Of(h, k) _ 0 
a ee 
In other words, the multiple points, if any, are those whose coordinates 
satisfy simulineously the three equations 
f (x, y) = 0, 
of (x,y) Of (x,y) _ 
“"- πὰ ΞΞ 0 a οὶ,.,»-- 0. 


or v’2+...=0.. 


ILLUSTRATION 11. 70 find the multiple points of the curve 

y® + 2 + 2ay — by? — 22+ 14y—11 = 0. 
The multiple points are those for which, on differentiating partially 
with respect to z, y in turn, 

22+ 2y—2 = 0, 
3y? + 2a—12y+14 = 0. 
Substituting for 27 in the second equation, we have 
3y?— 14y+16 = 0, 

or (y—2)(3y—8) = 0. 
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Fig. 143. 


oO AZ 


Fig. 144. 
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Thus y=2 or 8, 
so that, from the equation 27+ 20 -- 2 = 0, 

ἃ τὸ --1 or --. 


Thus possible multiple points are (--1, 2), (—%, 8). But it is easy 
to check that the solution z = —1, y = 2 satisfies the given equation 
whereas the solution x = —3, y = $ does not. Hence there is a 
multiple point at the point (—1, 2). 


5. Some rules for finding the asymptotes in simple cases. 
The method which we have given enables us to find the asymptotes 
of a curve f(x, y) = 0 with reasonable ease, and also (what is often 
very important) to estimate how the curve lies in relation to them. 
There are, however, one or two rules available in simple cases. 

The equation of the curve is 


1, ϑ = Un(%, Y) + Una(%,y) +... = 0, 
where u,(z,y) is a polynomial homogeneous of degree ὦ in 2, y. 
We begin the process of finding the asymptotes by factorizing the 
polynomial w,,(z, y) in the form 
U, (2, y) = (y — ky x) (y— kya)... (y—k, 2) 
and we restrict ourselves to the case when the constants k,, ko, ..., ky 
are real and distinct. With these restrictions, we can formulate the 
following rules: 
(i) Tum “PARTIAL FRACTIONS’ RULE FOR ΑΒΥΜΡΤΟΤΈΕΒ. If 
Un_1(2, Y) [Up (x, ψ), 
when expressed in partial fractions, assumes the form 
Un—1(%, y) αι ας an 
me - τς - 9, Ὁ 
ὠ,(α,}) y-ke y—kyx y—k,x 

(where it is assumed also that u,_,(x,y), u,(x,y) do not have common 
factors), then the equations of the asymptotes are 

y—kz+a,=0 (¢=1,2,...,n). 

The given equation, after division by w,,(z, y) is 
14 nal, Y)  Un—2(%, ψ) πολ δ 
U, (2, y) U,(%, y) 


Xn Un_o(%, υ ) 
ι ρα k,% U,(2, ¥ ) 


or Ppa ek 
she? 


+...=0, 
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Substitute the approximation 
y=k,x+BuP (p<1) 
in this form of the equation. Then 
1 +345 + SESE EE ἘΠ... Ἢ ΞΕ Ee 
Un_o(X, kx + Ba?) rome 
U,(x, kx + Bu?) ; 
Since p<1, and since none of k,—ky,...,4,—k,, is zero, the sig- 
nificant terms for large values of x are 


m4 πῇ 
1+ pi, +--0, 


so that, for closest approximation, we take 
Pp = 0, B = — ἃ... 
Hence yk, 2z-—a, 
eading to the asymptote 
y—k,x+a, = 0; 
and similarly for kg, ks, ...,k,. 
(ii) THE “TERMS OF HIGHEST DEGREE’ RULE FOR ASYMPTOTES. 
If f(x,y) can be expressed in the form 
(y—k,x+c,) (y — kav +.) een (y—k,%+ Oy) +Un_o(%, y) = 0, 
where, after the product of the n linear forms y — ka + ας, the rest of the 
equation consists of terms of degree n—2 at most, then, provided that 
ky, kg, ...,k, are all different, the lines 
y—kx+a, =0, 
y = kex + Ae = 0, 
y—k,x+a, =0 
are asymptotes of the curve f(x,y) = 0. 
Consider the approximation 
y=k,x+Bar (p<). 
On substituting, we have 
(Bu? +a) {(ky—ky)at...}... (6 —k,) Ot. Ἐν 9(@, kya +...)0, 
or 
(Ba? + αι) {(ky — he)... (hy — hey) 27 +. FA Vy o(@, Kyat...) =0. 


12-2 
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The terms of highest degree are (with p< 1) 

(k, —ky) ... (ky —k,) Bu? +" + (hk, — ky) ... (ky —&,) ae" 1 +... =0, 
and so (since k,, ka, ...,k, are distinct) we have for closest approxi- 
mation the relations 
p=0, B=-m%, 
so that the approximation is 


y=kx—o%, 
leading to the asymptote 
y—kx+a, = 0; 


and similarly for kg, ks, ..., ky. 
ILLUSTRATION 12. Τὸ examine the asymptotes of the curve 
y® — Gay? + Lla®*y — by? +27 + γ᾽ -- 2. Ὁ 8 = 0. 
Wehave τρία, y) = (y— 2) (y— 22) (y— 32), 
Un—i(%, ψ) =O? +y*, 


Un—1(%,Y) _ x+y? 


iain u,,(x,y) - ψ --αὴ (y—2z) (y— δε) 


Hence the asymptotes are 
y— x+1=0, 
y—2x—5 = 0, 
y—32+5=0. 
Moreover, we have the identity 
(y—x+1)(y—2x—5) (y—3x+5) 
= y* — 6xy* + Lla®y — 6y3 + 2? + y? + 802 — 25y — 25, 
so that the equation of the curve may be expressed in the form 
(y—x+1) (ν -- 2. -- δ) (y—3x+5) —32x+ 25y+ 28 = 0, 


corresponding to the ‘terms of highest degree’ rule just enunciated. 
(iii) ASYMPTOTES PARALLEL TO THE ΑΧΕΒ. If the equation of a 
given curve is expressed in terms of descending powers of y in the form 


y""w,(2e) + "M00, 4(tt) +... + Wal) = 0, 
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where w,(x),w,,3(%),...,W,(x) are polynomials in x of degree 
r,r+1,...,n, then, PROVIDED THAT THEY ARE ALL DISTINCT, the 
r straight lines 
w,(x%) = 0 
are asymptotes of the curve; and similarly with x, y interchanged. 
On division by y”~, the equation assumes the form 


ta) + Se 4, Pa) 0, 


which, for large values of y, approaches closely to the form 
w,(x) = 0. 
ILLUSTRATION 13. To find the asymptotes of the curve 
xy? — 3xy? + a+ 2y2?+a—y = 0. 
Since U(x, y) = αΞν3, 


there are no asymptotes except those parallel to the axes. The 
coefficient of the highest power of y (namely y?) is 


x* — 3a + 2, 
and the coefficient of the highest power of z (namely, z?) is 
y* +1. 

The latter expression does not give rise to real asymptotes, and so 
the asymptotes are 
x—l=0, 

z—2= 0. 

6. The radius of curvature at a multiple point. The method 
given earlier (Vol. II, p. 114) for calculating curvature can be 
adapted to meet the case where the curve has a multiple point. 


For we have shown how to obtain, for each branch, an approxi- 
mation to the equation of the curve in the form 


y=mxz+Axz” (n>1), 
Hence we have, for this branch, the approximate relations 


dy oy 
qa mt Aa” A 


FY n(n —1) Ag, 


174 THE SKETCHING OF CURVES 


At the origin, = is zero if n>2 and meaningless (‘infinite’) if 


n < 2. We may therefore confine ourselves to the normal case n = 2, 
when 
dy\ _ dy\ 
aaah: ire 
The curvature x is given by the formula 
dy 
da? 
K= Γ jdt ; 
t+ (@) 


dx 
ΡΒ Ὁ 
(1 την)" 


ILLUSTRATION 14, Τὸ find the radu of curvature of the branches 
at the origin of the curves 


(i) 2a*—3ay+y*?+a°y+y>+a4 = 0, | 
(ii) x?+at-—yt+y5 = 0. 
(i) We obtained (p. 149) the approximations 
Y= 2+ 22%, 
y= 2. -- 102". 
For the first, y’=1+42, y”=4, so that 
4 
Κ τε απ = ,2: 
and, for the second, y’ =2—20z, y” = -- 20, so that 


— 20 
plan ir ie ie 


so that 


(ii) We obtained (p. 152) the approximations 
easy", 
= —y?, 


Interchanging the roles of x, y we have 
dx 
dy ΞΞ Ὁ 2y ᾽ 


d*x 
dy" 
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so that 


7. The sketching of curves from their polar equations. 
We now outline very briefly a few of the considerations to be kept 
in mind when sketching a curve of given polar equation F(r,@) = 0, 
where 7 may be positive or negative and where 0 may take any 
positive or negative value (not necessarily restricted to an interval 
of 27). We confine our attention to examples in which the equation 
may be solved to give r as a single-valued function of @ in the form 


r= f(0). 


It seems harder to give a systematic 
treatment for polar coordinates than for 
Cartesians, and there is much scope for 
ingenuity. A routine first step may well 
be to construct a table of values of r for 
particular values of 9; for example, if 
f(@)involves the trigonometric functions, 
it is useful to put 9 = 0,47, ὑπ, ἐπ, ἐπ and 
so on. 

The formulae which we obtained in 
Volume II, Chapter X, may be used, 
and one or two of the results are re- 
peated here for convenience. 

Let P be the point (r, @) of the curve. The direction of the tangent 
at P is determined by the angle ¢ ‘behind’ the radius OP; if a radius 
vector centred on P is imagined to rotate about P from the position 
PO until it is first in line with the tangent at P, the angle of rotation 
is denoted by ¢, and it is proved that 


Fig. 145. 


dé 
tan d = ἔπ» 


This relation defines ¢ uniquely, since it lies between the values 
0,7. 
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For the curvature x, we obtained the formula 


dr\? dr dr\?\% 
wi, ede αἰ cea en ee τοῦ 
ἐπί () τα + (ao) |» 
the square root implied in the denominator being positive. Now 
(Vol. 11, p. 114) the sign of x is the same as the sign of Ὕ when x is the 
parameter defining the curve; and (Vol. 1, p. 54) the concavity is 
‘upwards’ or ‘downwards’ according as is positive or negative. 
In the present work, the parameter is 0, and the relationship 
connecting the concavity with the sign of the above formula for « is 
more complex. It remains true, however, that the curve has an 
2 
inflexion where x = 0, that is, where ΣΡ +o = 0. 
Note that the expression on the left assumes a simpler form on 
replacing r by its reciprocal τ = 1/r. Then 
ΕΝ ΟὙΙ τος 
dd ωδάθ’ ἀθ ωδάθδ᾽ ὡδλαᾶθ) ’ 
so that the expression is, on substituting and simplifying, 


du 
(τ + u) / us, (1) 
In particular, the points of inflexion satisfy the equation 
+u=0 
9 
provided that the expression (1) changes sign for such a value of 0. 


Singularities 

The location of singularities is, of course, very helpful. The origin 
will certainly be a multiple point if the equation f(@) = 0 has more 
than one (significantly) distinct solution. Suppose, more generally, 
that there is a singularity at the point (p,a), where p = f(«). 
Since r is a single-valued function of 6, this can happen only if 
EITHER the value of r corresponding to a value a+27, a+ 4n, 
a +67, ... of θ is also p, on if the value of r corresponding to a value 
atin, at37,a+57,... of 15 —p. For example, if f(9) is a trigo- 
nometrical function of 0, it is often (but not always) sufficient to 
consider the points, if any, expressed in the alternative forms 
(p,a), (—p,a+7); and these are found by solving the equation 

f(a) +f(a+m) = 0. 
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Thus, for the curve r = sin? + cos 20, 


we should consider the equation 


{sin θ + cos 26} + {sin (8 + 77) + cos (20 + 277)} = 0 


or cos 20 = 0, 
iving ie Re 
gl 7 oe 4° 4? 4? “ee 


and singularities at the two points 


ae ee ἘΠ σεν 
Ae’ ar 1 a a ν᾽ Π)Ξ:(-- Ὥ.- 
Asymptotes 

The identification of the (linear) asymptotes is a fairly simple 
process, but the justification of the method seems less easy. For 


convenience, we write 1/r = u, 1/f(@) = g(@), so that the equation 
of the curve assumes the form 


u=9(9). 
For large values of r, its reciprocal w is small, and so we begin by 


finding the solutions, if any, of the equation g(?) = 0. Let 0=a 


be such a solution, where 
g(x) = 0. 


Take an associated system of Cartesian coordinates in which the 
positive x,y axes are the radii 0 = 0, ἐπ respectively. Let ὶ be a 
value of 0 near to a, and let @(5, #) be the corresponding point of 
the curve (fig. 146), so that 


1 
3 = 9(6). 
The line through Q in the direction ἃ is given by the equation 


xsina—ycosa = ὦ βὶη (α -- β) 
_ εἰπία-- β) 
9() 

Now suppose that β approaches the value a, and suppose that, 
as it does so, : 

ial SER) νὰ 

fora g (2 ) 
Then the line takes up the limiting position given by the equation 


xsina—ycosa =h, 
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The distance of the point Q (bcos f, ὃ βίῃ f) from this line, being 
+ (bcos £ ΤΡ cos -- ἢ) 
ἘΞ ἢ mat ἡ, 
g(P) 
tends to zero as β tends to a. Thus the line whose equation is 
lim 518 sin (a — β) 


xsina—ycosa = 
9 yes ΜΠ ΗΝ 


Fig. 146. 


is an asymptote of the curve, being approached ever more closely 


by the point Q of the curve as f tends to a. 


In accordance with the usual rule for the limiting value of the 
ratio of two functions each of which tends to zero, we have the 


relation sim Si (4-8) dite~ ἘΝ: cos (a—f) 
pow J9 : δῆς. σ΄ (ῥῚ) β-α 
mune 
g'(x) 
Thus, in normal cases, the equation of the asymptote is 
—1 


xsina—ycosa = ——. 
ε 9. (α) 
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InLustRaTIon 15. ΤῸ find the asymptotes of the curve 
r(cos @—sin 0) = 1+cos0+sin@. 
cos 0 —sin 0 
1+cosé+sin@° 
The asymptotes are obtained from the equation 
cos @—sin @ = 0, 


If 1/r = u, then u= 


so that θ = jn, ἐπ, .... 
If 6 = jz, then, since 
(0) = —sin@—cos@ _—_(cos#—sin6@)?* 
9 Ὁ) = (+ οοϑθ βίη 6) (1+cos0+sind)?’ 
we have 9 (dn) = —2/(1+/2), 
and so the equation of the asymptote is 
LI aA RS, 
J2 5 2 ’ 
or x—y=/2+1. 
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The value 0 = $7 gives x—y = —./2+1, and other values of 0 


give repetitions of these two asymptotes. 
ILLUSTRATION 16. 70 find the asymptotes of the curve 
r /{cos 0(cos @—sin #)} = 1. 
If l/r=u, then w= ,/{cos0(cos 6 —sin 6)}. 
The asymptotes are obtained from the equations 
cos@=0, cos@—sind = 0, 


so that @ = ἐπ, $2, ...; θ.- ja, ἔπ, .... 
If 9 = ἀπ, then 


ps ee 
join vos βίοοϑ β-- sin β}} ~jin. loon πὶ “i 


Similarly for 0 = ὅπ, Thus one asymptote is the line 


x= 0. 
θ = jn, then 


—. (cos f—sin f) 
lim sin (}7 — β) = lim V2 = 
prin νίοο8 f(cos 8—sin f)} ὁ pin /{cos B(cos # — sin 4) } 


=}, ("soar |~° 
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Similarly for 6 = ὅπ. Thus the other asymptote is the line 
z—y=0. 
(We could not use the formula involving g’(«), since g’(47), σ΄ (477) 


are both ‘infinite’, The answer would have been correct, but the 


justification might not.) 


Note. It may be remarked that, in general, if the equation of 


the curve occurs in the form 


rG(9) = F(0), 
then g(@) = ray 


G0) _ (0) F'(8) 
Τὴ) {Γ0}} " 
Hence the asymptotes occur when 6 = a, where 


σ΄ (ΘᾺ = 


G(a) = 
G'(a 
so that 9΄ (α) =, Fay 
and the equation for the asymptote is 
xsina—ycosa = τ ° 


This is the form often given, and it can be useful; but it needs care, 
Thus the curve ἢ 


"= [cos 0(cos 9: sin 8}} 
just discussed leads us, with 
F(@)=1, G(@) = ψίοοβ 6(cos 6 —sin 6)}, 
into theoretical difficulties. On the other hand, with 
F(@) = /{cos @(cos 0—sin 6)}, G(0) = cos 6(cos6—sin 6), 


we avoid these difficulties, at the expense, perhaps, of some mental 
confusion about the apparently arbitrary selections possible for 
F(@) and G(@). 


Finally, it should be remembered that the curve 


r= f(0) 


may have a circular asymptote if f(@) tends to a finite limit as 6 tends 
to infinity. Thus, if , 
᾿ {ἴα [{6) = 
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then the curve (of ‘spiral’ shape) approaches asymptotically the 
circle 


r= a. 
For example if r = tanh@, 
e? —¢~# 
then, since lim tanh 9 = in F368 =], 
θ--» ὦ eo +e~9 + eo 
the asymptotic circle is fei, 


REVISION EXAMPLES XIV 


‘Scholarship’ Level 
1, The equation of a curve is 


xy? — 22+? = 0, 
(i) Find the equations of the tangents at the origin. 
(ii) Find the equations of the real asymptotes. 
(iii) Show that the numerical value of y is always less than 
that of x. 
(iv) Show that the numerical value of y is always less than unity 
(v) Sketch the curve. 


2. The equation of a curve is 
(a? +y?)? = 2*(1+y?). 
(i) Discuss the symmetry of the curve with respect to the 
origin and the axes. 
(ii) Determine the coordinates of the points at which the 
tangents are parallel to the z-axis. 
(iii) Prove that the curve has no real asymptotes. 
(iv) Find the points of intersection of the curve with the axes, 
and show that x is always numerically not greater than 1. 
(v) Sketch the curve. 
3. Find the asymptotes of the curve 
αϑ -- αὐ -τὰῷ = 0, 
and prove that the origin is a point of inflexion. 
4. Find the tangents at the origin andthe asymptotes tothe curve 
w(x* — y*) = γία + 2y). 
Prove that each asymptote meets the curve at a finite point 
which lies on the straight line 


185 - ὃν = 6. 
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5. Trace the curve 

x + xy? +2%—y*? = 0, 

Draw the tangents to the curve at the origin and find the equation 
of the real asymptote. 

6. Find the asymptote of the curve 

a—y® = ay(x+y). 

There are two other tangents to this curve which are parallel 

to the asymptote. Find their equations and the coordinates of 


their points of contact. 
Trace the curve. 


7. Sketch roughly the curve 


6y? = x(x —y) (w—3y); 
find its asymptotes and indicate them on your sketch. 
Show that each asymptote meets the curve in a single finite 
point and that the three points so obtained are collinear. 


8. Find the asymptotes of the curve 
(z+y—1)% = a+y%, 
and prove that they meet the curve only at infinity. 


Prove also that there is no point on the curve for values of x 
between A and 1, where A is the real root of the equation 


3A8 + 3A2—3A+1 = 0. 
Give a sketch showing the general form of the curve. 


9. Find the equation of the straight line which is asymptotic to 
Prove also the following facts and give a sketch of the curve: 
(i) the origin is a cusp; 
(ii) no part of the curve lies between x = 0 and x = 4; 
(iii) the curve consists of two infinite branches, one lying in the 
first quadrant and the other in the second and third quadrants. 
10. Trace the curve 
(x? —y?)?—4y? +y = 0. 
11. Find the asymptotes of the curve 
(y—«)? (y +2) —(y—2) (ἐν -- α) + 2y—a = 0. 
Investigate on which side of each asymptote the curve lies, and 
trace the curve. 
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12. Find the asymptotes of the curve 
(x+y) (~—2y) (w—y)* + ϑαγία -- ὐ) Ἐῶ γῆ = 0. 


Determine on which side the curve approaches each of the ends 
of the parallel asymptotes. 


13. Find the asymptotes of the curve 
e+ 222° +a = ry*?—Qry+y. 


Show that there are values which y—-2 never takes, and sketch 
the curve. 


14, Find the asymptotes of the curve 
228 — δχϑυ — xy? + 6y® — 4.3 — 4ary + 15y?— 62+ 9y+1 = 0. 


Give a rough sketch of the curve, indicating clearly on which 
side of the asymptotes the curve lies. 


15. Give a rough sketch of the curve 
x — xy? — 24° + 3ry + νυ -- ὃν = 0, 
and indicate clearly on which sides of the asymptotes the various 
branches of the curve lie. 


16. Trace the curve 
Qay? + 2(.3 — 2+ 2) y—(x?-—5a+2) = 0. 
Prove that at no finite real point of the curve is the tangent 


parallel to the z-axis. 


17. Determine the asymptotes of the curve 


(y— 1)? (y?— 4a?) = 3ay. 
Investigate on which sides of the asymptotes corresponding 
branches of the curve lie, and trace the curve, 


18. Trace the curve 
y?(x® + 23) = 3(.3 — 4a + 2.5}, 
19. Show that «—y = 3 is an asymptote of 
(x—y+1)(e—y—2) (x+y) = 82--ὶ, 
find the other asymptotes, and sketch the curve, 
20. Find the asymptotes of the curve 


a(a+y) = «2+ 4y, 
and trace the curve. 
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21. Trace the curve 
αὖ —a*z*y —b®zy?+y5 = 0, 
What modification occurs when a? = b2?2 


22. Trace the curve 
(x? — y?)? a3a = 0, 


23. Find the asymptotes of the curve 
at + 32%y + 2x?y? + Qary + 3x+y = 0; 


determine on which side or sides the curve approaches each 
asymptote, and where it cuts the asymptotes. 


24. Trace the curve 
y® — 3*y + 4a(x*—y*) = 0, 


finding its asymptotes and the curvatures at the points where 
x = 0 meets it. 


25. Find the equation of the tangent at the point (1,2) to the 
curve given by ay(et+y) = 2244241, 
and determine the point at which it intersects the curve. Find the 
asymptotes and trace the curve. 


26. The lines whose equations are x = y, x =—y, x = 2y are 
the asymptotes of a cubic curve which touches the axis of x at 
the origin and which passes through the point (0,6). What is the 
equation of the curve? 


27. Sketch the curve 2z?—y? = 7, 


Find (i) the position of the centre of curvature of either branch of 
the curve at the origin, and (ii) the area of the loop. 


28. Sketch the curve zy = 23+), 


and find (i) the radii of curvature at the origin of coordinates, 
(ii) the area of the loop. 
29. Investigate the curvature at the point (1, 2) of the curve 
(y—2)? = x(a—1)2, 
30. Trace the curve 
4(x? + 2y? — 2ay)? = χϑίαϑ- Qy?), 
and find the radii of curvature of the two branches at the origin. 
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31. Find the asymptotes of the curve 

ary Ἕν" = x? — 4y?, 
and trace it. 

Find the cubic which has x+y = 1 as an asymptote and touches 
both axes at the origin, the radii of curvature there being 1 and 2 
units in length. 

32, Find the asymptotes of the curve 

are 
ai—as 


o = 
and find the radius of curvature at the origin. 
Sketch the curve. 
33. Find the curvature at the origin of each of the branches of 
the curve δ: ἀρ ay 4g? = dey — Dy", 
and trace the curve. 
34, Find the radii of curvature at the origin of the curve 
22? + xy — 3y? + 223 —a*y + 4y4 = 0. 
35. Evaluate [. V{(b—x)/(x—a)}dz (a<b) 
by means of the substitution x = asin? @ +6 cos? @, or otherwise. 
Make a rough drawing of the curve 
a? + 3xy?— 3a(x*—y?) = 0, 
and show that the area of its loop is 3a?. 
36. Find the coordinates of the node of the curve 
(et+yt+l)y+(x+y+1}%+y' = 0, 
and the area of the loop at the node. 
37. Make a rough sketch of the curve 
x(a? + y?— 3a?) = 2α8, 


and prove that the area between the curve and its asymptote is 
37a", 


38. Sketch roughly the curve 
y*(a +22) = a(a*—2%), 
and find the area of one of its loops. 
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39. Trace the curve 


(2? +y*)? = l6axy?, 
and find the areas of its loops. 
Prove that the smallest circle that will completely circumscribe 
the curve has radius 3 /3a. 


40. Trace the curve y? = χϑία -- αὶ 
for a= 1,0,—1. Find the area enclosed by the loop in the case 
a=-l. 
41. Trace the curve 
16a*y? = b®2*(a — 22), 


where a,b are positive, and find the area enclosed by the loop. 
If 16a? = 3b?, show that the perimeter of the loop is 40. 


2/2 
42. Trace the curve y?= = = a - 
and find the area of the loop. 
43. Trace the curve 
a + y% —3axy = 0, 


and find the area of the loop. 


44, Give a rough sketch of the curve 
+ y5 = 5ax*y?, 
Find the area of the loop of the curve. 


45. Sketch the curve whose equation in polar coordinates is 
r = 1-Έ 008 20. 
Prove that the length of the curve corresponding to 0< 6 < 27 is 
8 + log (2+ 4/3). 
46. Sketch the curve whose polar equation is r? = a*(1+3 cos 0) 
and find the area it encloses. 


47. Trace the curve r = a(2cos@—1), find the areas of its loops, 
and show that their sum is 37a?. 


48. Trace the curve r = a(cos@+cos 20), and show that the 
curve crosses itself at the points where r = 4a /2, 0 = +4n. 

Prove that the area of that portion of the largest loop that is not 
common to the other loops is $a? /2. 


RE VISION EXAMPLES XIV 187 


49. Trace the curve r = 2+3cos 20, and find the areas of the 
loops. 


50. Trace the curve r = a(sin 6 —cos 20), and find the area of the 
loop which passes through the point (2a, }7). 


51. Sketch the curve 
γ(1 -- 2 0080) = 3acos 20, 
and find the equations of its asymptotes. 
52. Trace the curve 
rcos@+acos 30 = 0. 


Show that the area of the loop is a*(2— 47), and that the area 
enclosed between the curve and its asymptote is a*(2 + 47). 


53. Determine the asymptotes of the curve rcos3@ =a, and 
sketch the curve. 


54. Sketch the curve r(cos@+sin 0) = asin 20, and find the area 
of the loop of the curve. 


55. Trace the curve rcos@ = asin 30, and prove that the area 
of a loop is 4a*(9 /3— 4m), 
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Examples IX: 
. 4090 . . 69 . . 11-19 . 
ANSWERS TO EXAMPLES ἦν OZ eee /o 
CHAPTER XIV 
Examples X: 
Examples I: 1, 2y?, ἀψὺ, 3.5, 4y, 4a. 
1, Continuous. 2,3,4. Discontinuous, 2. Gary*, 12a*%y8, 12a%y?, 24ry®, 362%y?, 
3. e*siny, e* cosy, —e*siny, e* cosy, —e*siny. 
Examples 11: 4. 0, y, —ay-2, 0, —y-. 
1. aay”, Saty*. 2. ys πέρ ἐν 5. 2siny, 2acosy, —z*siny, 2cosy, —2zsiny. 
3. sINYy, TCOSY. 4. —ysinay, —xsin zy. 6. 0, —e%, ze7,0,e%. 7. 0, —y-*, 2ay-%, 0, ὃν" 3. 
5. e* cos 2y, — 265 sin 2y. 6. ὁδὶ ὃν 1 ox 
7. A;0 8.) SS ἘΞ ΞΞΞΙ, Ὁ, “a ag» © 
1,0. (1+a2)? (1+2%)’ ? {τ 
9. hy/ Ju, Jz. (2a? + 8% + 4) ye®, (2a? + 41) ες. 
10. 3a%ytz-*, 4a8y8z-2, — 2a3ytz-3, 11, 425, 3y?, 22. 
12. e*siny cosz, e* cosy cosz, —e*siny sinz. | Examples ΧΙ: 
13. 008 ἃ; siny sinz, sinx cosy sinz, sin# siny cosz, 9a 4a2 ox a 
14. 3(1+2)% ev, 2(1+2)% eM, y(1+2)% eM, i. ψ’ yp 2. ~ By? By 
ἐπ Xz xy ae Darke 
ac eS μα 1+ yz?’ 1+ yz: 3. yey 
2eyt+y? by(x+y) (αὔτ αν ἐν) 
Examples IV: 4, Day + at 2 (a+ 2y)® ὺ 


1. 2y, 2%, 2rsin 20, 2r cos 20. 
2. 32", — ὃν, 3(2° — y?)/r, — 8xy(x%+y)/r. 


Examples VI: REVISION EXAMPLES X 
1. 16035, 2. (cos?#—sint) esint, 
3.0. 4. 219(2 + 345) 2, 8, v(x, y) =e7(xsiny+ycosy), F(z) =ze*. 
26, Suu Fu 
Examples ΚΠ: S29y—Ly Ix 
1. cos 20 sin? ¢ cos ¢, sin@ cos @ sin ¢(2 cos? ὁ — sin? 4). 36. $bsinC, bsinC, jacot A. 
2. 2(a+y+z) —2xy*2?, ete. 44. (i) 2, (1) 1. 
3. e*+V{cos (x —y) —sin (x —y)}, e*+”{cos (ὦ -- ψ) +sin (x—y)}. 45. u=ax+b, v=acy+d, w=cz. 46. Asinh»v. 
4. 6u(u?-+-08)Flog (u?— 08) + MOE EE 49, ὁ | + (cob B+ 2 cot 4)α]. 
6v(u® + v®)2log (u®— v2) — 2v(u? + υϑ)8 50. p(y) = (ay +6) (ὁψ Ὁ 41 —(a'y - δ (cy+d). 


w—ye Take y=eJ?~*4v, with similar form for X. 
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REVISION EXAMPLES XI 
1, Max. } at (2,1); min. — at (-- 8, —}). 
3. 8, 4, 8. 4, c*/(a?+b?+1), 
5. Max. at (A, A)if A? > a*; min. at (A, A) if A2< a?; no stationary 
value at A? = a?, 
6. 4/3 abc. 
7. Roots of X/*(br?—1) (cr?—1) = 0. 
8. Min δα at points such as (— 4a, 2a, 2a). 
[(a, 0,0), (0, α, 0), (0, 0, α) are not genuine points. ] 
11. syk, th. 
12. Greatest is (La??P-2))—-2/, Least is the least of a, ὦ, c. 
16: 1.2. ¥, 
16. Max. 27 at (2, —1, -- 4); min. 3 at (0, —1, -- 3). 
19. Neither max. nor min. at (—5, 4, 4); min. at (i, I, 1} 


abe (»-- 1)(4-- 1) (»-- 1) 
*°: ποκα τα }* Σ(4-- 1)(ν-- 1) * 
2 y-2 A 


Examples IT: 
l. ν. 2. οἷν 3. fr. 4, 0. δ, e*sin 0, 


Examples III: 
2. No. 3. Yes. 4, α -Ξ -- ἃ. 


REVISION EXAMPLES XII 
A-y 
4 V{—(a+A) (ats) (b+A) (b+4)} 
2(% + 2y +1) (z+1) 
e+y+2 ; 
8. 2ve*t*/(1 + 3.0). 


CHAPTER XVII 
Examples I: 


1. Gp. 2. 3gmk. 4, 327k. 


ANSWERS TO EXAMPLES 191 


REVISION EXAMPLES XIII 


1, π--4. 2. (38a, 3%). 3. 0; (e”—1—ab)/b. 
4, 22q4, Sat, δα, 6. τ αϑδ3, 7. oe 
8, log 2+4%. 10. παϑ. 11. ἐπαῦδο. 

12, 2045, 13. ὦ τὸ ψϑα, υ = αϑ!ν. 


14, 2(3π-- 4)α3|9θ. 1δ. maid, 
90. | “49 | Meet T ACen — 10) ei, 


in 0 
21. 4:7abce(b?+c*). 25. ,a%bp (7 + “1 : 


18. {sm (a + b + c). 


CHAPTER XVIII 


Examples III: 
l. y= +2(1—42). 2. y= +2(1— 22?) 
. ©=+y(1+2y). 4, x=+y(1—2y’). 


3 
5. y = 2x— 42" or y = 42°, 
6. y= —4e—2 ory = 2, 


Examples IV: 

Ι,͵ 2v—y+1=0,2-y+1=0. 
x—y=0,r+y+2=0. 
x+y+4=0,2-y+2=0, 2. --ἶὖ -- ὃ τ 0. 
2+yt+4=0, r—y+ $= 0, x-y-1=0. 
. @-yt+2=0,r+y=0,x+y4+2=0.° 


om 9 


Examples V: 
1, 2(y+a—})*+3a = 0. 2. (22+})? =y. 
ὃ. (y+3x+3)? = 9x, 


REVISION EXAMPLES XIV 
1. (i) zty=90, (i) yt 1 =0. 


2. (ii) (0,0), (+=, +=), (iv) (0,0), (+1,0). 
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. £=0,r+y=0. 


4, Tangents y = 0, ++ 2y = 0; 


. 5/2, 385,4/18. 


Asymptotes x+2 = 0,7+y-—}4=0,x-y-8=0. 


. «=1, 
. &©-y= fa; x—y—a=0 at (0, --α and z—-y+ja=0 


at (—§a, $c). 


. ©=2,x2-y+3=0, x-—3y-1=0. 

. *=ly=l1rz+y=0. 9. x-y+1=0. 

. Asymptotes are rt+y+1=0. 

. ©-y+1=0,x-y+}4=0,r+y-$=0. 

- e+y—3 = 0, r—-2y+2=0, r-y—2=0, r—-y+}4=0. 
. ©=0,r-y+2=0,r+y=0; 


y—x is not between — 2,0 or 2, 4. 


. *+y+1=0, x—2y—3 = 0, 27+ 3y = 0. 
. Asymptotes arex—1 = 0,4+y—2=0,2-y+1=0, 


Asymptotes are x = 0, y= }4,x+y-4=0. 


. ¥-1l=0,yt2x=0, 19. w—-y+2=0,r+y=0. 


z+y=0,7+2=0. 


. Asymptote is x+y =0. If a? = b?, the line x+y = 0 is 


part of the locus. 


. Asymptotes x+y = 0. 
. «= 0, cutting at (0,0); x+y = 0, cutting at (0,0), (1, -- 1), 


«+ 2y = 0, cutting at (0,0), (3, —4). 


. Asymptotes y=4a, 4+3,/32—-3y+4a=0; curvature 


—% 2a each branch at origin and a— at (0, 4a). 


. 6r+y = 8, (ὦ, —5); asymptotesz = ly = l,a+y+2 = 0. 
. (z?—y?) (w—2y) = 2b*y. 
. (2, —2), (-- 2, — 3); 5%. 


4; 2, 29. 2,-3.2. 80. 2a, a. 


. £+4=0,y-1=0,r%+y-3=0; 


day = 223 + 9x°y + Bry? + γϑ, 


. z=+a,y=0; ja. 
3. ὁ for y=0; — 5/5 for y= 35. The asymptote is 


«+y =—%, and asymptotic parabola (y—2x+4)? = 2. 
35. 4m(b—a). 


36. 
38. 


SERBS 


49, 
50. 
53. 
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(-- 1,0); 4. 37. The asymptote is ὦ = 0. 

(4a —1)a*. 39. 47a’. 

. 41. dab. 

3/3. 43. 3a. 

δα, 46. a*{7+2./2—sin— (% /2)}. 
. aX(7—$./3), a2(27 +3 V3). 


Ai tan (/5)+ ὃ V5; ar —44¢ tan (/5) — ὃ 5. 
(ἐπ - ἧς /3) a’. 51. α“ ,ϑεν-α 3. 
ϑα-ατεῦ, Ἐν γ8-τ- ξα. 54. 1α3(π-- 3). 
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INDEX 


Algebraic curve, 142 
Asymptotes, 158, 170 

circular, 180 

parabolic, 162 

parallel, 158 

in polars, 177 
Asymptotic parabola, 164 


Chain rule, 25 
for Jacobians, 97 
Change of variable, 41 
by operators, 43, 46 
Coefficients, partial, 10 
of higher order, 38 
a 052 


spherical polar, 2, 129 
Curve, algebraic, 142 
multiple point, radius of curvature, 
173 
tangent, 50 
See also Sketching of curves 
Cusp, 150 


Differential, 28, 31 
exact, 51 
Differentiation, partial 
introductory example, 8 
notation, 10 


‘Element of area,’ 121 
‘Element of volume,’ 127 


Geometrical interpretation, 4 
of differentiation, Cartesian, 13 
of differentiation, polar, 16 
dx 300 
of Pa δε “" 18 


Homogeneous functions, 47 


Implicit functions, 30, 49 

Increment 
for several variables, 11 
small, 36 

Inflexion, 167 

Integrals, multiple, 109 
double, 109 


evaluation in practice, 118 
as repeated integrals, 116 
triple, 115 

Integrating factor, 54 


Jacobian 
chain rule, 97 
definition, 95 
of dependent functions, 100 
introductory example, 90 
for ratio of areas, 102, 133 
for ratio of volumes, 104, 134 
‘reciprocal’ theorem, 98 


to distinguish, 77 
general conditions, 74 
undetermined multipliers, 79 
Moment of inertia, 110 
of cardioid, 125 
Multiple point of curve, 154 
to locate, 168; in polars, 176 
radius of curvature, 173 


Operator, 43, 46 
Parabola, asymptotic, 164 


Sketching of curves 
asymptotes, 158, 170 
circular, 180 
parabolic, 162 
in polars, 177 
cusp, 150 
general method, 143 
inflexion, 167 
large values of z, 154 
location of multiple points, 168; in 
polars, 176 
from polar equation, 175 
small values of x, 146 
tacnode, 151, 154 


Tangent, 50 
Taylor’s theorem, 57 


Volume 
in cylindrical coordinates, 127 
of sector of sphere, 128 
in spherical polar coordinates, 129 
See also Integrals, multiple (triple) 


μὴ 
ee 
Ε- 
τς 
Ο 
ΌΔ 
Ζ 
° 
—_ 
Ζ 
— 
-- 
° 
~” 
n 
aa 
m 
Ξ- 
[58 
Σ 
fe) 
Ὁ 


Besise 
ane 


ἫΓ 


3s 


ion 


scope 
Te Tinset Ditacdiousd Sermlecions 


such is beyond the 


HE 
Hite 


standard. No 
wise to omit 
bg eee nlc 
Journal of the A.M.A. 


and lava 


11} 


ae ἢ 


commended... .The series 


racively presented and We look forward to the 


Journal of Education 


ja 


He 
late 


ae 


Bas z 


BY ΕΚΑ. MAXWELL 


THE METHODS OF PLANE PROJECTIVE 
GEOMETRY, BASED ON THE USE OF GENERAL 
HOMOGENEOUS COORDINATES 
(available in both cloth and paperbound (Students’) editions) 


GENERAL HOMOGENEOUS COORDINATES 
IN SPACE OF THREE DIMENSIONS 


ADVANCED ALGEBRA 
IN TWO VOLUMES 


AN ANALYTICAL CALCULUS FOR 


SCHOOL AND UNIVERSITY 
IN FOUR VOLUMES 


FALLACIES IN MATHEMATICS 
(available in both cloth and paperbound editions) 


ΠῚ TOA snynoyed yeonAjeuy Uy : TTA MX VW 


Published by 


CAMBRIDGE UNIVERSITY PRESS 


Bentley House, 200 Euston Road, London, N.W. 1 
American Branch: 32 East 57th Street, New York, N.Y.10022 


517 


PRINTED IN GREAT BRITAIN 


Ὁ 521 05698 5 


